This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Allele Frequency reported by Samtools and IGV Do Not Agree

Hi all,

I use the mpileup command to create the pileup file with this command :

samtools mpileup -v -u -f file.fasta mybam.bam -l my bed.bed > myvcf.vcf

The problem is that When I visualize a region using IGV, the allele frequency percentage reported by IGV does not agree with that reported by the mpileup command (the QS values) in any given position. Has anyone else faced this problem? Is there a solution?

Thank you,

frequency igv allele snp-calling
• 1,757 views • link
updated by shiyang_bio • shiyang_bio • written by OmBibo • OmBibo

1 answer

It's possibly a result of the filter in samtools. By default samtools discards some reads with low quality while IGV utilized all reads. You may check the parameters for samtools and include all reads.

• 0 views • link
written by shiyang_bio • shiyang_bio

Thanks Shiyang for your reply. For samtools, I used --min-BQ 0 to disable filtering of low-quality bases. But no changes, it shows the same values as outputs. Please if you have any idea about other parameters or filters to consider for samtools? Your help is really appreciated.

• 0 views • link
written by OmBibo • OmBibo

Log in to answer this question.

More posts like this