This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Finding Snps From A Bam File

Hi everyone

I have dna seq dataset for S.cerevisiae genome. I have aligned bam files and the ref seq genome. I would like to know some methods for finding snps from the bam files as input.

I know 2 ways to do it.

  1. Using samtools mpileup and then bcftools and generating a vcf file.

  2. Using GATK unified Genotyper

Any other method if you could tell would be really nice

Hope to hear from you soon

Regards

Varun

snps
• 8,814 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by Varun Gupta Varun Gupta

3 answers

This question has been addressed before : What methods do you use for In/Del/SNP calling?

• 43 views • link
written by Gabriel R. Gabriel R.

There are many tools for finding SNPs but the ones you mentioned are the most used or popular ones. So if I would be in your place I will just use those tools (samtools and GATK). I used them to call SNPs for a particular strain and found 95% of the SNP calls in common.

One of the many posts (search for 'SNP calling' on bistars.org) that discuss the snp calling

What methods do you use for In/Del/SNP calling?

• 42 views • link
written by Ashutosh Pandey Ashutosh Pandey

Ok let's have fun with the semanticmedia wiki interface of http://seqanswers.com

The semantic data of the wki can be downloaded from

http://seqanswers.com/wiki/Special:ViewXML?title=Special%3AViewXML&categories%5BBioinformatics_application%5D=on

the following stylesheet extracts the application having input=BAM and output=VCF


<xsl:stylesheet xmlns:xsl="&lt;a href="http://www.w3.org/1999/XSL/Transform" "="" rel="nofollow">http://www.w3.org/1999/XSL/Transform'
xmlns:str="http://exslt.org/strings"
extension-element-prefixes="str"
version='1.1'
>

<xsl:output method="text"/>
<xsl:template match="/">
<xsl:apply-templates select="/Pages/Category/Page/Template[@Name='Bioinformatics_application']"/>
</xsl:template>

<xsl:template match="Template">
<xsl:variable name="pageName" select="../@Title"/>
<xsl:variable name="input" select="Field[@Name='input format']"/>
<xsl:variable name="output" select="Field[@Name='output format']"/>
<xsl:for-each select="str:tokenize($input,', ')">
<xsl:if test="./text()='BAM'">
<xsl:for-each select="str:tokenize($output,', ')">
<xsl:if test="normalize-space(./text())='VCF'">
<xsl:text>http://seqanswers.com/wiki/</xsl:text>
<xsl:value-of select="$pageName"/><xsl:text>
</xsl:text>
</xsl:if>
</xsl:for-each>
</xsl:if>
</xsl:for-each>
</xsl:template>

</xsl:stylesheet>

result:

$ xsltproc stylesheet.xsl   "http://seqanswers.com/wiki/Special:ViewXML?title=Special%3AViewXML&categories%5BBioinformatics_application%5D=on"

(BTW I don't know well SMW , does anyone know if it's possible to query the RDF data directly from the wiki ? is there a SPARQL point ?)

• 4 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this