This is a test version of Biostars. For the public version, visit https://www.biostars.org.
filtered.bcf file much larger than raw.bcf file

After using mPileup (Sam Tools v. 1.6) I get files that are in the neighborhood of 12 GB. Next, I use Bcftools filter (bcftools v. 1.6) to filter out variants with read depth less than 2 or greater than 100. The filtered files are around 85 GB! Why are the filtered files so much larger? After all I am filtering out sequence right? Here's the bcftools filter script

bcftools filter --threads 20 raw.bcf -e 'DP<2 || DP>100' > filtered.bcf

Any thoughts would be appreciated. Thanks

sequencing alignment
• 2,125 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by PKSWEE11 • PKSWEE11

1 answer

your generating a VCF file (plain text), not a binary (small) BCF.

see options:

-o, --output <file>           write output to a file [standard output]
-O, --output-type <b|u|z|v>   b: compressed BCF, u: uncompressed BCF, z: compressed VCF, v: uncompressed VCF [v]
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thanks! That was very helpful

• 0 views • link
written by PKSWEE11 • PKSWEE11

Log in to answer this question.

More posts like this