This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Bowtie2 wrong alignment score in local mode with ignore-quals?

Used bowtie2 2.2.6. I have the following 2 alignments, with exactly the same CIGAR. The only difference is that there is a mismatch in the second alignment. My scoring options are: --rdg 2,1 --rfg 2,1 --score-min L,18,0 --np 1 --ma 1 --mp 5,1.

For the first alignment AS = 21 (M) - 2 (XO:i:1) - 15 (XG:i:15) + 22 (M) = 26, which matches.

For the second alignment AS = 21 (M) - 2 (XO:i:1) - 15 (XG:i:15) + 21 (M) - 5 (XM:i:1) = 20, which is diffrent from AS.i:24 in the alignment. Is there something I'm missing?

Alignment 1: tag_344761|1 0 MI0000167_mmu-mir-142a 4 255 6S21M15D22M * 0 0 CTTTCTCATAAAGTAGAAAGCACTACTTGTAGTGTTTCCTACTTTATGG IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII AS:i:26 XN:i:0 XM:i:0 XO:i:1 XG:i:15 NM:i:15 MD:Z:21^AACAGCACTGGAGGG22 YT:Z:UU Alignment 2: tag_457953|1 0 MI0000167_mmu-mir-142a 4 255 6S21M15D22M * 0 0 GTCTCTCATAAAGTAGAAAGCACTACTTGTAGTGTTNCCTACTTTATGG IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII AS:i:24 XN:i:0 XM:i:1 XO:i:1 XG:i:15 NM:i:16 MD:Z:21^AACAGCACTGGAGGG9T12 YT:Z:UU

alignment sam bowtie

0 answers

No answers yet.

Log in to answer this question.