Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 621,170 C→T 41.0% A51V (GCG→GTG)  ycaI → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE621,1700CT41.0% 41.3 / 66.3 95A51V (GCG→GTG) ycaIconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (41/15);  new base T (8/31);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 4.02e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

TGCGTTTCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAA  >  minE/621105‑621235
                                                                 |                                                                 
tGCGTTTCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCt                                                                >  1:2540422/1‑69 (MQ=255)
tGCGTTTCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCt                                                                >  1:2389319/1‑69 (MQ=255)
tGCGTTTCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCt                                                                >  1:1838668/1‑69 (MQ=255)
      tCTGACTCTCTTCTCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          >  1:3022729/1‑69 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          >  1:3048533/1‑69 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          >  1:2662812/1‑69 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          <  1:2600189/69‑1 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          >  1:782987/1‑69 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          <  1:668677/69‑1 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          <  1:1982816/69‑1 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          >  1:1976398/1‑69 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          <  1:1951368/69‑1 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          >  1:916286/1‑69 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          <  1:544071/69‑1 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          >  1:1774373/1‑69 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          <  1:1099592/69‑1 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          <  1:1715458/69‑1 (MQ=255)
      tCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGtt                                                          >  1:995989/1‑69 (MQ=255)
                              ggCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCAtatt                                  <  1:317643/69‑1 (MQ=255)
                                    tATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCag                             >  1:2496062/1‑68 (MQ=255)
                                    tATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCag                             >  1:3123505/1‑68 (MQ=255)
                                    tATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCag                             >  1:519209/1‑68 (MQ=255)
                                    tATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCag                             >  1:3158563/1‑68 (MQ=255)
                                    tATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCag                             >  1:3067077/1‑68 (MQ=255)
                                    tATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCag                             >  1:927663/1‑68 (MQ=255)
                                    tATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCag                             >  1:2757115/1‑68 (MQ=255)
                                    tATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCag                             >  1:1416615/1‑68 (MQ=255)
                                    tATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCa                              >  1:2845236/1‑67 (MQ=255)
                                    tATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATAag                             >  1:2425583/1‑68 (MQ=255)
                                    tATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTCTTTTTCGTTTGGGGCATATTATCag                             >  1:313855/1‑68 (MQ=255)
                                        ccTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCagcaaa                         >  1:801494/1‑68 (MQ=255)
                                        ccTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCagcaaa                         >  1:2802747/1‑68 (MQ=255)
                                        ccTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCagcaaa                         >  1:2908160/1‑68 (MQ=255)
                                        ccTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCagcaaa                         >  1:2078721/1‑68 (MQ=255)
                                        ccTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCagcaaa                         >  1:2963504/1‑68 (MQ=255)
                                        ccTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCagcaaa                         >  1:158448/1‑68 (MQ=255)
                                        ccTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCagcaaa                         >  1:1803480/1‑68 (MQ=255)
                                        ccTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCagcaaa                         >  1:1666281/1‑68 (MQ=255)
                                        ccTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCagcaa                          >  1:2843152/1‑67 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTTACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:408283/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:3218887/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:3180989/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:631276/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:3254395/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:37185/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:3324679/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:669172/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:3417434/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:795180/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:1012951/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:1126306/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:997671/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:1471838/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:1495769/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:95382/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:1715845/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:947785/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:931501/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:1862643/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:1911870/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:1946462/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:826356/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:2051552/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:2139479/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:2209974/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:227752/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:2359183/62‑1 (MQ=255)
                                            ttAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:2448296/62‑1 (MQ=255)
                                             tAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:1149840/61‑1 (MQ=255)
                                             tAAAACCGTCCGTTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAgcagca                           <  1:1141688/61‑1 (MQ=255)
                                              aaaaCCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAAtt                  >  1:899778/1‑69 (MQ=255)
                                              aaaaCCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAAtt                  >  1:998316/1‑69 (MQ=255)
                                              aaaaCCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAAtt                  >  1:1066897/1‑69 (MQ=255)
                                              aaaaCCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAAtt                  >  1:1216772/1‑69 (MQ=255)
                                              aaaaCCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAAtt                  >  1:337101/1‑69 (MQ=255)
                                              aaaaCCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAAtt                  >  1:1805683/1‑69 (MQ=255)
                                              aaaaCCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAAtt                  >  1:1824501/1‑69 (MQ=255)
                                              aaaaCCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAAtt                  >  1:2489609/1‑69 (MQ=255)
                                              aaaaCCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAAt                   >  1:394140/1‑68 (MQ=255)
                                              aaaaCCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAAt                   >  1:2010576/1‑68 (MQ=255)
                                                      ccgtTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAGCagcaaagc                       >  1:1908807/1‑56 (MQ=255)
                                                      ccgtTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAGCAGCAAAGCAAATTTTGTGGGc          >  1:2714586/1‑69 (MQ=255)
                                                      ccgtTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAGCAGCAAAGCAAATTTTGTGGGc          >  1:2675930/1‑69 (MQ=255)
                                                      ccgtTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAGCAGCAAAGCAAATTTTGTGGGc          >  1:2229861/1‑69 (MQ=255)
                                                      ccgtTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAGCAGCAAAGCAAATTTTGTGGGc          >  1:1932597/1‑69 (MQ=255)
                                                      ccgtTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAGCAGCAAAGCAAATTTTGTGGGc          >  1:1918785/1‑69 (MQ=255)
                                                      ccgtTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAGCAGCAAAGCAAATTTTGTGGGc          >  1:335184/1‑69 (MQ=255)
                                                      ccgtTATATCGTGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTAGCAGCAAAGCAAATTTTGTGGGc          >  1:1285982/1‑69 (MQ=255)
                                                             atCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGgaga    <  1:2720329/68‑1 (MQ=255)
                                                             atCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGgaga    <  1:2426488/68‑1 (MQ=255)
                                                             atCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGgaga    <  1:796984/68‑1 (MQ=255)
                                                             atCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGgaga    <  1:553332/68‑1 (MQ=255)
                                                             atCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGgaga    <  1:3104213/68‑1 (MQ=255)
                                                             atCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGgaga    <  1:1225350/68‑1 (MQ=255)
                                                             atCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGgaga    <  1:345790/68‑1 (MQ=255)
                                                               cgcgCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGaaa  <  1:1513147/68‑1 (MQ=255)
                                                                 |                                                                 
TGCGTTTCTGACTCTCTTCGCCTGCGTACTGGCATTTATCCCTGTTAAAACCGTCCGTTATATCGCGCTGACGTTGCTGTTTTTCGTTTGGGGCATATTATCAGCAAAGCAAATTTTGTGGGCAGGAGAAA  >  minE/621105‑621235

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