SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BFAST and Variant Calling nexgengirl Bioinformatics 8 01-02-2013 06:03 AM
variant calling with SAMtools krueml Bioinformatics 2 09-28-2012 06:36 AM
samtools mpileup did not call variant lchang1101 Bioinformatics 2 09-24-2012 08:22 AM
Problems using GATK for variant calling pooja_singh RNA Sequencing 1 07-07-2012 11:23 AM
variant calling kjaja Bioinformatics 1 11-04-2011 08:16 AM

Reply
 
Thread Tools
Old 03-31-2013, 09: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
Reply

Tags
bioinformatics, mpileup, samtools, variant calling

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:35 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO