Dear All, I am not so sure about two Tophat settings. Please help.
1) "Number of mismatches allowed in the initial read mapping" Based on the documantation, my understanding is: the reads are re-aligned to transcriptome/genome if the mismatches in the initial alignment is more than the set number (for example, the default setting is 2). In other words, the re-aligning will continue until the mismatches is equal to or below the set number. Is my understanding correct? If I am right, I have one worry: will Tophat stop re-aligning if the mismatch is below 2 (if I use the default setting). If it is true, the read will not be aligned to where it belongs to (with 0 mismatch). 2) "Number of mismatches allowed in each segment alignment for reads mapped independently" Does this mean that the reads will be cut into segments if the mismatches of alignment is more than the set number? Thanks in advance. Jianguang
___________________________________________________________ The Galaxy User list should be used for the discussion of Galaxy analysis and other features on the public server at usegalaxy.org. Please keep all replies on the list by using "reply all" in your mail client. For discussion of local Galaxy instances and the Galaxy source code, please use the Galaxy Development list: http://lists.bx.psu.edu/listinfo/galaxy-dev To manage your subscriptions to this and other Galaxy lists, please use the interface at: http://lists.bx.psu.edu/