Question: Finding Snps From A Bam File
1
gravatar for Varun Gupta
7.1 years ago by
Varun Gupta1.1k
United States
Varun Gupta1.1k wrote:

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 • 4.1k views
ADD COMMENTlink modified 7.1 years ago by Pierre Lindenbaum129k • written 7.1 years ago by Varun Gupta1.1k
1
gravatar for Gabriel R.
7.1 years ago by
Gabriel R.2.7k
Danmarks Tekniske Universitet
Gabriel R.2.7k wrote:

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

ADD COMMENTlink written 7.1 years ago by Gabriel R.2.7k

You pasted it before me. You won. :-)

ADD REPLYlink written 7.1 years ago by Ashutosh Pandey12k
1
gravatar for Ashutosh Pandey
7.1 years ago by
Philadelphia
Ashutosh Pandey12k wrote:

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?

ADD COMMENTlink written 7.1 years ago by Ashutosh Pandey12k
1
gravatar for Pierre Lindenbaum
7.1 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum129k wrote:

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 ?)

ADD COMMENTlink written 7.1 years ago by Pierre Lindenbaum129k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1518 users visited in the last hour