Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 264,681 A→G 34.8% P193P (CCA→CCG thiL → thiamin‑monophosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE264,6810AG34.8% 43.0 / 36.5 66P193P (CCA→CCGthiLthiamin‑monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base A (21/22);  new base G (11/12);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TGCAGGTTGCCGATGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGCCATCGATCTCT  >  minE/264626‑264742
                                                       |                                                             
tGCAGGTTGCCGATGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTAttt                                                  <  1:2121360/69‑1 (MQ=255)
tGCAGGTTGCCGATGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTAttt                                                  <  1:651493/69‑1 (MQ=255)
tGCAGGTTGCCGATGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTAttt                                                  <  1:3116952/69‑1 (MQ=255)
tGCAGGTTGCCGATGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTAttt                                                  <  1:3223385/69‑1 (MQ=255)
           gATGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggca                                        >  1:2117721/1‑68 (MQ=255)
           gATGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggca                                        >  1:778256/1‑68 (MQ=255)
           gATGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggca                                        >  1:235038/1‑68 (MQ=255)
           gATGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggca                                        >  1:176353/1‑68 (MQ=255)
           gATGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggc                                         >  1:1075671/1‑67 (MQ=255)
            aTGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggcagg                                      <  1:1417450/69‑1 (MQ=255)
            aTGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggcagg                                      <  1:1453255/69‑1 (MQ=255)
            aTGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggcagg                                      <  1:478853/69‑1 (MQ=255)
            aTGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggcagg                                      <  1:1625237/69‑1 (MQ=255)
            aTGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggcagg                                      <  1:2585835/69‑1 (MQ=255)
            aTGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggcagg                                      <  1:507498/69‑1 (MQ=255)
            aTGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggcagg                                      <  1:1878919/69‑1 (MQ=255)
            aTGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggcagg                                      <  1:1185909/69‑1 (MQ=255)
            aTGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggcagg                                      <  1:1117959/69‑1 (MQ=255)
              gCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggcaggca                                    >  1:2534668/1‑69 (MQ=255)
                     aTGCTGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGcgcg                             <  1:3122812/69‑1 (MQ=255)
                     aTGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGcgcg                             <  1:2558106/69‑1 (MQ=255)
                     aTGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGcgcg                             <  1:2500467/69‑1 (MQ=255)
                     aTGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGcgcg                             <  1:956933/69‑1 (MQ=255)
                     aTGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGcgcg                             <  1:1565049/69‑1 (MQ=255)
                         ggACTACTTGATCAAACGTCATCTCCGTCCGTCGCCGCGTATTTTACAgggg                                          >  1:831445/1‑52 (MQ=255)
                         ggACTACTTGATCAAACGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCt                         >  1:2018596/1‑69 (MQ=255)
                         ggACTACTTGATCAAACGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCt                         >  1:822534/1‑69 (MQ=255)
                         ggACTACTTGATCAAACGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCt                         >  1:697626/1‑69 (MQ=255)
                         ggACTACTTGATCAAACGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCt                         >  1:1781843/1‑69 (MQ=255)
                         ggACTACTTGATCAAACGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCt                         >  1:422017/1‑69 (MQ=255)
                         ggACTACTTGATCAAACGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCt                         >  1:3049892/1‑69 (MQ=255)
                         ggACTACTTGATCAAACGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCt                         >  1:2614291/1‑69 (MQ=255)
                         ggACTACTTGATCAAACGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCt                         >  1:2618405/1‑69 (MQ=255)
                         ggACTACTTGATCAAACGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCt                         >  1:143247/1‑69 (MQ=255)
                         ggACTACTTGATCAAACGTCATATCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGa                            >  1:3066101/1‑66 (MQ=255)
                               ctTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAgg                                            >  1:2000304/1‑44 (MQ=255)
                               ctTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCaaa                   >  1:849850/1‑69 (MQ=255)
                               ctTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCaaa                   >  1:2789525/1‑69 (MQ=255)
                                 tGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGggcagg                                      >  1:2900218/1‑48 (MQ=255)
                                 tGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCaaa                   >  1:1483407/1‑67 (MQ=255)
                                 tGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCaaa                   >  1:3332862/1‑67 (MQ=255)
                                 tGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCaaa                   >  1:42071/1‑67 (MQ=255)
                                 tGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCaaa                   >  1:1297452/1‑67 (MQ=255)
                                 tGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCaaa                   >  1:1034550/1‑67 (MQ=255)
                                 tGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCaaa                   >  1:2243192/1‑67 (MQ=255)
                                      aaaCGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGct            <  1:2738004/69‑2 (MQ=255)
                                      aaaCGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGct            <  1:197743/69‑2 (MQ=255)
                                      aaaCGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGct            <  1:2405977/69‑2 (MQ=255)
                                      aaaCGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGct            <  1:1080055/69‑2 (MQ=255)
                                      aaaCGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGct            <  1:2461575/69‑2 (MQ=255)
                                      aaaCGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGct            <  1:559322/69‑2 (MQ=255)
                                      aaaCGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGct            <  1:1357874/69‑2 (MQ=255)
                                      aaaCGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGct            <  1:1743962/69‑2 (MQ=255)
                                      aaaCGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGct            <  1:389836/69‑2 (MQ=255)
                                      aaaCGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGct            <  1:1642911/69‑2 (MQ=255)
                                      aaaCGTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGct            <  1:2713268/69‑2 (MQ=255)
                                         cGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGCCAt          >  1:2653703/1‑68 (MQ=255)
                                         cGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGCCATc         >  1:2437295/1‑69 (MQ=255)
                                         cGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGCCATc         >  1:2708732/1‑69 (MQ=255)
                                         cGTCATCTCCGTCCATCGCCGCGTATTTTAAAGGGGCAGGCACTgcgc                              >  1:949600/1‑48 (MQ=255)
                                          gTCATCTCCGTCCGTCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGct            <  1:2917084/65‑2 (MQ=255)
                                          gTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGCCATCg        <  1:1729436/69‑1 (MQ=255)
                                          gTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGCCATCg        <  1:1347707/69‑1 (MQ=255)
                                          gTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGCCATCg        <  1:2688720/69‑1 (MQ=255)
                                          gTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGCCATCg        <  1:2194457/69‑1 (MQ=255)
                                                tcCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGCCATCTAtctct  >  1:3287369/1‑69 (MQ=255)
                                                       |                                                             
TGCAGGTTGCCGATGCTAAAGATGCGGACTACTTGATCAAACGTCATCTCCGTCCATCGCCGCGTATTTTACAGGGGCAGGCACTGCGCGATCTGGCAAATTCAGCCATCGATCTCT  >  minE/264626‑264742

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: