This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Reference And Alternative Allele Count

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
• 7,561 views • link
updated by Erik Garrison Erik Garrison • written by gabrielle123123 • gabrielle123123

4 answers

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.

• 3 views • link
written by matted matted

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

• 4 views • link
written by Malachi Griffith Malachi Griffith

Thanks, i will look at those option

• 0 views • link
written by gabrielle123123 • gabrielle123123

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!

• 0 views • link
written by Malachi Griffith Malachi Griffith

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.

• 0 views • link
written by Erik Garrison Erik Garrison

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

• 0 views • link
written by Zev.Kronenberg Zev.Kronenberg

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

• 0 views • link
written by Erik Garrison Erik Garrison

Log in to answer this question.

More posts like this