Biostar Beta. Not for public use.
Filter pileup file for completeness
0
Entering edit mode
18 months ago
AP • 90

Hi all,

Is there a straightforward way to filter a pileup file for completeness using samtools (or other programs)? I would like to exclude sites based on the proportion of missing data (e.g. SNP present in 80% of the samples). Basically, I want to do the same than what vcftools --maxmissing does but on a pileup file instead of a VCF file.

Any help would be greatly appreciated! Thanks!

ADD COMMENTlink
0
Entering edit mode
18 months ago
AP • 90

6 months later: simply filter your pileup file using a grep command, removing lines that have a relatively low coverage. For instance: grep -vw '[0-5]' A.mpileup > B.mpileup will remove all lines with coverage below 5.

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1