SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Alternative Hits in BWA generated .sam (http://seqanswers.com/forums/showthread.php?t=5445)

jtladner 06-09-2010 02:04 PM

Alternative Hits in BWA generated .sam
 
I am working with Illumina sequences that are about 76 bp long. I have been looking through some .sam files generated by samse after aligning to a reference using bwa. I noticed that many of the reads that contain multiple optimal hits do not have XA tags at the end of the line, while many others do.

I have not noticed any consistent difference between reads that have the XA tag and those that don't. Can anyone offer insight into this?

Jason

nilshomer 06-09-2010 02:07 PM

Quote:

Originally Posted by jtladner (Post 19941)
I am working with Illumina sequences that are about 76 bp long. I have been looking through some .sam files generated by samse after aligning to a reference using bwa. I noticed that many of the reads that contain multiple optimal hits do not have XA tags at the end of the line, while many others do.

I have not noticed any consistent difference between reads that have the XA tag and those that don't. Can anyone offer insight into this?

Jason

Did you check the BWA manual?
Code:

samse        bwa samse [-n maxOcc] <in.db.fasta> <in.sai> <in.fq> > <out.sam>

Generate alignments in the SAM format given single-end reads. Repetitive hits will be randomly chosen.

OPTIONS:
-n INT        Maximum number of alignments to output in the XA tag for reads paired properly. If a read has more than INT hits, the XA tag will not be written. [3]



All times are GMT -8. The time now is 04:21 AM.

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