This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Convert A Soapsnp Output File To Soap, Sam Or Bam Formats

Hello,

I would like to know if there is any tool to convert SOAPsnp to SOAP? Once it is converted to SOAP, I can convert it to SAM using Soap2Sam from reseqtools. Hopefully using reseqtools is not that complicated.

If this is complicated, i wouldn't mind trying the reverse. Somehow converting SAM file to SOAPsnp files.

Background: I ran BWA in a HPC cluster. So I have SAM/BAM files from that test. Now I have to generate SAM/BAM files from the same data by running tests on Hadoop/AWS. I tried using DistMap. But ran into some errors; still trying to fix those errors. Now I am thinking of using Crossbow, which seems to be relatively easy to run on AWS. After that I have to compare both the results in terms of quality, time etc.

Hoping to find some help from Biostar members.

Thanks, Ashwin

sam bam genome mapping paired-end
• 3,358 views • link
updated by Istvan Albert Istvan Albert • written by bashwin.u • bashwin.u

1 answer

The SOAPsnp program' purpose is to assemble a consensus genome based on alignments and a reference genome.

The output of the SOAPsnp program cannot be converted to an alignment format.

• 0 views • link
written by Istvan Albert Istvan Albert

Thanks for replying. Then I will have to think about converting SAM file generated from BWA to Soapsnp to do a comparison of both architectures. I understand I can convert SAM to SOAP by using reseqtools. But I am not too sure how to convert it to SOAPsnp. I think I can do that by sorting the SOAP file using msort and then use SOAPsnp. Is this how I can do?

• 0 views • link
written by bashwin.u • bashwin.u

Log in to answer this question.

More posts like this
  • Sam To Bam
    written by freddy • freddy
    <p>I sequenced a mito genome using paired end method on the mi-seq I mapped all the reads back to reference gnome using <a href='http://bio-bwa.sourceforge.net/'>BWA</a> and …
  • How To Sort Sam File
    written by Chen Sun Chen Sun
    <p>I want to process the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file created by <a href='http://bio-bwa.sourceforge.net/'>BWA</a>, if it is sorted by the sequence coordinate on reference, it will be …
  • Sam To Bam - Loss Of Data Or Just Great Compression?
    written by Clare Clare
    <p>I use picard tools to take a <a href='http://bio-bwa.sourceforge.net/'>BWA</a> alignment <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file and convert it into a sorted <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. Normally this works …
  • How Can Soap2Sam.Pl Be More Fast?
    written by wm wm
    <p>I have some big soap files (&gt;20 G, each) and I need to convert them into <a href='http://samtools.sourceforge.net/SAM1.pdf'>SAM</a> format. </p> <p>I am using the soap2sam.pl …
  • Samtools View Of A Sam File Entire Chromosome
    written by dfernan dfernan
    <p>Hi, </p> <p>I have (what I think) is a valid <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file.</p> <p>When I do:</p> <pre><code>samtools view -S F1i1F.C57.2.Aligned.out.sam chr19_C57 </code></pre> <p>I get the …
  • How To Interpret Crossbow Data
    written by bashwin.u • bashwin.u
    <p>Hello,</p> <p>May I know how to interpret Crossbow output. Is there any chance of building a <a href='http://samtools.sourceforge.net/SAM1.pdf'>SAM</a> or <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file from Crossbow output. …
  • Merge The Sam Files For Dexseq
    written by camelbbs camelbbs
    <p>I have RNAseq data for one sample, with two technical replicates. As I know, I need to merge them to get the gene or exon …
  • Sam To Soap Format Converter
    written by Casey Bergman Casey Bergman
    <p>I am interested in testing <a href='http://soap.genomics.org.cn/soapsnp.html'>SOAPsnp</a> on alignments made from from <a href='http://bio-bwa.sourceforge.net/'>BWA</a> (and other mappers), but SOAPsnp only takes <a href='http://soap.genomics.org.cn/soap1/#Formatofoutput'>SOAPaligner output format</a> …
  • Soap Result Converted To Bam Not Working In Samtool
    written by Sarwar • Sarwar
    <p>I converted the output of soap2 which is map file to <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file using map2sam.pl and then converted this <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> to <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> …
  • Bioscope With Cufflinks
    written by Ramon Vidal • Ramon Vidal
    <p>Does any one here have ever used Bioscope with Cufflinks to do DE analysis? I am trying to use the <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> files generated by …