SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   CIGAR field of SAM file from Bowtie alignment (http://seqanswers.com/forums/showthread.php?t=19846)

qwsqe 05-07-2012 12:37 AM

CIGAR field of SAM file from Bowtie alignment
 
In a SAM file generated from Bowtie, I found one entry has '*' in the CIGAR field. What does this mean, bad alignment?

Thanks

arvid 05-07-2012 12:44 AM

It means "not available" in the SAM spec. Typically, this happens when the read isn't aligned (RNAME - 3rd field - is also '*' and 0x4 is set in the FLAG field).

qwsqe 05-07-2012 05:31 AM

Thanks, arvid.
To what extend, Bowtie sets the FLAG to 0x4 for a read, and make the CIGAR field "not abailable". What are the measure and threshold bowtie uses for setting this?

Thanks

arvid 05-07-2012 05:39 AM

Well, using the settings (default or your specific choices), Bowtie found no reportable alignment for that read.
It is not possible to make a general statement about this behaviour here as it depends entirely on the settings you provided bowtie with and the nature of your reads and reference. Check the bowtie manual, typically you should look at alignment and reporting modes - these and their sub-options control alignment sensitivity and reporting filters.

qwsqe 05-08-2012 05:28 AM

Hi Arvid,

Thank you very much for the help.


All times are GMT -8. The time now is 12:14 PM.

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