SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   BWA: specifying SAM/BAM file header fields before read alignment? (http://seqanswers.com/forums/showthread.php?t=5078)

nora 05-11-2010 06:59 AM

BWA: specifying SAM/BAM file header fields before read alignment?
 
Dear all,

I have found the modification of SAM/BAM file headers quite cumbersome and time-consuming, so I wonder if it is planned to include an option into BWA to specify header fields before starting the alignment, or if there is a supereasy, superfast way of adding header lines to a SAM or BAM file without overwriting the beginning of the alignment?

Currently I am splitting the output SAM file into header and alignment, modifying the header, and concatenating them back together. This needs about 15 minutes for a 40 million reads SAM-file, which seems a bit much knowing that I only want to add a couple of alignment info lines.

Many thanks!

drio 05-11-2010 10:11 AM

I don't think there are plans to incorporate that feature anytime soon. Feel free to email the mailing list.

Dump your alignments in BAM format and process the header with picard. You should reduce the processing time.

nora 05-12-2010 09:33 AM

Thanks for your advice!

gerald2545 12-04-2010 09:11 PM

Hello
did you try "samtools reheader"?
Never used, but seems to do what you want

Gerald


All times are GMT -8. The time now is 03:54 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.