View Single Post
Old 03-31-2013, 08:53 PM   #1
asebastian
Junior Member
 
Location: USA

Join Date: Jul 2012
Posts: 4
Default Variant Calling with mpileup

Hi all,

I am a new bee to variant calling. I am using samtools mpileup. My criteria includes:
-- minimum total coverage of 5
-- minimum variant-supporting coverage of 2
-- minimum average base quality of q15
-- minimum variant allele frequency of 0.01

This is how I am doing it

samtools mpileup -DS -C 50 -m 2 -F 0.01 -Q 15 -f ref.fa test.bam | bcftools view -bvcg - >test.bcf

bcftools view test.bcf | vcfutils.pl varFilter -d 5 >test.vcf

Is this the way to do it?

Also how do I find "mean position of the variant within the variant-supporting reads" ?

Thanks in advance for all your comments and suggestions.

-A
asebastian is offline   Reply With Quote