This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Bwa: Should Sai/Sam Sequences Processed Match?

Hi,

After running BWA the output includes information such as this:

[bwa_sai2sam_pe_core] print alignments... 0.66 sec
[bwa_sai2sam_pe_core] 197538574 sequences have been processed.
[main] Version: 0.6.1-r104

Should the no. of sequences processed match for the sai and sam files?

For both sai files I have the no. of sequences processed at 147913358

For the sam file it is 197538574

I ran this again this but got the same result. For other sai to sam conversions this number seems to match.

Thanks,

Jom

sam bwa
• 2,864 views • link
updated by Biostar Biostar • written by jomaco jomaco

Unfortunately I specified the wrong read files at the sai > sam stage. These read files contained 197538574 reads instead of 147913358.

Despite my error, I do wonder how it was possible for the conversion to continue with the wrong reads files being specified. What is the purpose for the read files being specified again at the sampe stage? The process seemed to complete successfully despite my error.

Thanks for taking the time to look at my question.

• 0 views • link
written by jomaco jomaco

1 answer

The help for bwa sampe seems to indicate that it has the ability to output multiple hits for paired reads:

-n INT   maximum hits to output for paired reads [3]
-N INT   maximum hits to output for discordant pairs [10]
• 0 views • link
written by Istvan Albert Istvan Albert

If there were more lines than sequences in the sam file (ignoring the header) then your answer would account for that. However, I am not sure how this would explain why BWA has apparently processed more sequences to produce the sam file. If it said 197538574 alignments have been processed then this would make sense. Even so shouldn't the number still be the same at the sai stage?

• 0 views • link
written by jomaco jomaco

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 …
  • Conversion From Sam To Bed With Max Information
    written by roll roll
    <p>I am trying to convert my <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file into bed file in such a way that i will keep some information from the <a …
  • Removing Phix Data From Ngs Data For Analysis
    written by crivenster crivenster
    <p>I have a NGS data containing 95-98% of PHIX data,which i dont need. So after some reading, I used <a href='http://bio-bwa.sourceforge.net/'>BWA</a> to align the reads …
  • Seeing Unexpected Characters (^D,^Q) In The Qual Field Of A Sam File
    written by edwardhust • edwardhust
    <p>i use <a href='http://bio-bwa.sourceforge.net/'>bwa</a> tool to map raw pair end reads to the reference, then output two sai file. then, i use <a href='http://bio-bwa.sourceforge.net/'>bwa</a> sampe …
  • 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 …
  • Do The Sequences In Sam/Bam Header Have To Match The Sequences In Reference Fasta File?
    written by Victor Lyuboslavsky • Victor Lyuboslavsky
    <p>Do the sequences in <a href='http://samtools.sourceforge.net/SAM1.pdf'>SAM</a>/BAM header have to match the number and the order of the sequences in reference FASTA file?</p> <p>For example, if …
  • Bwa Sampe Segmentation Fault
    written by a11msp a11msp
    <p>Hi everyone,</p> <p>I'm running <a href='http://bio-bwa.sourceforge.net/'>bwa</a> in the sampe mode and, after successfully processing &gt;10M reads it fails with a segmentation fault (as follows) on …
  • Predict/Estimate/Find Bwa Parameters From Bam Or Sam File
    written by Houkto • Houkto
    <p>Hi All, </p> <p>I would like to know if there is a way of finding out the parameters of <a href='http://bio-bwa.sourceforge.net/'>BWA</a> aligned <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a>/SAM file. …