This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Missing @Pg Header In Bwa Alignment

Hi All,

I am using BWA Version: 0.7.5a-r405 for alignment. The output sam file has only @SQ lines in the header. The @PG line is missing. For now I added the line manually, but I was wondering if there is a way that the aligner can output the @PG line by itself in the sam file.

Thanks for your help.

bwa
• 4,292 views • link
written by is327 • is327

2 answers

Strange, I'm looking a the souces of bwa-0.7.5a: and book sampe/samse sources have a call to bwa_print_sam_PG()

void bwa_print_sam_PG()
{
    err_printf("@PG\tID:bwa\tPN:bwa\tVN:%s\n", PACKAGE_VERSION);
}

and as far as I can see err_printf writes to stdout...

update : the PG flag was added after vr405:

http://sourceforge.net/p/bio-bwa/code/51/tree//trunk/bwa/main.c?diff=43

+void bwa_print_sam_PG()
+{
+    printf("@PG\tID:bwa\tPN:bwa\tVN:%s\n", PACKAGE_VERSION);
+}
+
• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

Hi, I'm brand new to this community and I'm also brand new to bioinformatics so I apologize if this is considered as resurrecting an old thread but it's exactly the problem I'm having. I'm afraid I don't understand the answer, does it mean that I have an outdated version of bwa installed, which doesn't include the @PG line? Unlike, is327 I don't have the knowledge or experience to manually add the header.

Thanks :)

• 0 views • link
written by smurfben67 • smurfben67

But for some reason it is not written on the sam file.

Thanks.

• 0 views • link
written by is327 • is327

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 …
  • 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 …
  • Clarification For Running Bwa In Parallel‏
    written by Mohamed • Mohamed
    <p>Hi,</p> <p>I tried running <a href='http://bio-bwa.sourceforge.net/'>BWA</a> in parallel by splitting the input FASTQ file into 'n' splits and ran single ended alignment (aln and samse …
  • Cufflinks Or Alternative From Bwa Alignments?
    written by Biostar User • Biostar User
    <p>Can I use Cufflinks starting from reads in <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a>/bam format aligned using <a href='http://bio-bwa.sourceforge.net/'>bwa</a> instead of bowtie? If not, what is the alternative to …
  • Problem Running Pindel With Bam Files
    written by xenophiliuslovegood • xenophiliuslovegood
    <p>Hi, I am running <strong>Pindel version 0.2.4p, March 26 2012</strong>, on a set of <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> files. It chokes on them, printing out:</p> <p><strong>No currentState.Reads …
  • Sam To Bam Conversion Reference Fasta File
    written by Atom Smasher • Atom Smasher
    <p>Hello,</p> <p>I had some data in "eland_export.txt" format and I want to convert it to the <a href='http://samtools.sourceforge.net/SAM1.pdf'>SAM</a> format first and then later to the …
  • Preprocessing The Bam For Gatk
    written by Pierre Lindenbaum Pierre Lindenbaum
    <p>The <strong>GATK</strong> often requires that the <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> contains some extra informations in the header ( Group ID, Group Library, SO:coordinate...).</p> <p>Which commands I should …
  • 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. …