Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,824,410 T→A pseudogene (543/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,824,4100TA100.0% 359.6 / NA 100pseudogene (543/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (0/100);  total (0/100)

GGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGG  >  W3110S.gb/3824360‑3824452
                                                  |                                          
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCTGCAGAGACCGa                              <  1:2302635/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:2475643/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:1988391/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:2097455/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:1823125/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:2239796/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:2324566/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:2467261/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:2472355/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:1569692/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:840629/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:1528191/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:1526187/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:255712/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:354794/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:1298475/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:719821/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:1098970/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:1125899/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:1017396/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:371384/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:707387/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:1129509/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:499360/62‑1 (MQ=255)
ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa                              <  1:386268/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGTAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:1876502/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:2094628/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:807364/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:20918/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:2210415/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:2000684/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:757492/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:734208/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:2337599/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:2445217/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:690385/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:68371/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:638138/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:547968/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:205324/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:2532162/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:360475/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:844713/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:1040896/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:1244355/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:1260859/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:1429590/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:1508215/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:1585419/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:83959/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:1857751/62‑1 (MQ=255)
            ggtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGACGCAGAGACCGACATGAACGTAGt                  <  1:41716/62‑1 (MQ=255)
             gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:1571057/62‑1 (MQ=255)
             gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                  <  1:1574072/62‑1 (MQ=255)
                     atttGCGATCTGGAATACGTTGAAGACTCAGCGGCAGAGACCGACATGAACGTAGt                  <  1:2320809/55‑1 (MQ=255)
                            aTCTGTAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:1065819/65‑1 (MQ=255)
                            aTCTGGAATACTTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:1474692/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGTCATGAACGTAGTGATGACCGAAGACggg  <  1:1841236/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGATATGAACGTAGTGATGACCGAAGACggg  <  1:1126202/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCTAAGACggg  <  1:1865828/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:180726/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:824482/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:729838/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:1052965/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:2254461/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:1129123/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:119054/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:874967/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:53809/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:123527/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:2278339/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:400835/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:813317/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:1385961/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:1848904/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:342832/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:2313002/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:2318739/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:2495596/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:789937/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:1897292/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:175943/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:1688057/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:1841162/65‑1 (MQ=255)
                            aTCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:1820955/65‑1 (MQ=255)
                                ggtaTACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg  <  1:2260502/58‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGTCATGaa                        <  1:1231286/35‑1 (MQ=38)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:98037/35‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:782381/35‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:709122/35‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:31580/35‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:2493904/35‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:2426272/35‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:2414542/35‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:2243523/35‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:1593315/35‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:1568769/35‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:1522694/35‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:1220251/35‑1 (MQ=255)
                                    tACGTTGAAGACTCAGCCGCAGAGACCGACATGaa                        <  1:1093712/35‑1 (MQ=255)
                                                  |                                          
GGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGG  >  W3110S.gb/3824360‑3824452

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑