Reference And Alternative Allele Count
4
2
Entering edit mode
12.2 years ago

Hi, I'm working on polypoid plant. Where using samtools and BCF tools to call SNP. We realize that bcftools have difficulty calling some simplex genotype and i would like doing it on my own. Is there a way to extract the number of reference and alternative allele from a mpileup format?

Thanks,

Gabrielle

reference allele • 6.6k views
ADD COMMENT
1
Entering edit mode

It might be worth checking out other variant calling software. For example I have been using SNVer for pooled and non diploid data. http://snver.sourceforge.net/

ADD REPLY
3
Entering edit mode
12.2 years ago
matted 7.8k

See this post to get most of the way there using samtools mpileup. The other suggestion to consider other SNP callers is good too. I recall that FreeBayes can handle arbitrary ploidy, but I don't have any personal experience with it.

ADD COMMENT
0
Entering edit mode
12.2 years ago

You can use Bio::DB::Sam to get reference and alternative allele counts directly from a BAM file.

ADD COMMENT
0
Entering edit mode
12.2 years ago

Thanks, i will look at those option

ADD COMMENT
1
Entering edit mode

Hi. Welcome to BioStar. For future reference, this comment should have been added as a comment to one of the answers or the original question. Also, the best way to thank people for their answers is by voting up the answers you find useful by pressing those triangle buttons along the left hand side of the page. Good luck and welcome to the community!

ADD REPLY
0
Entering edit mode
11.9 years ago
Erik Garrison ★ 2.4k

Rather than hacking samtools to handle polyploid detection, I would suggest using FreeBayes. It's designed to detect SNPs, indels, and complex variants (such as MNPs and other block substitutions) without restriction of ploidy.

I'm the author, so please email me if you have any questions.

ADD COMMENT
0
Entering edit mode

Can you call multiple bam simultaneously? IE strait to a multi individual vcf? Are you in the Marth lab?

ADD REPLY
0
Entering edit mode

Yes, I'm in the Marth lab. FreeBayes goes directly from (potentially 10s of thousands of) BAMs to a multi-sample VCF.

ADD REPLY

Login before adding your answer.

Traffic: 1543 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6