Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 264,912 T→G 45.9% T270T (ACT→ACG thiL → thiamin‑monophosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE264,9120TG45.9% ‑1.3 / 54.4 48T270T (ACT→ACGthiLthiamin‑monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (22/4);  new base G (15/7);  total (37/11)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTG  >  minE/264848‑264977
                                                                |                                                                 
ttGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGc                                                               >  1:2791594/1‑69 (MQ=255)
ttGAACCGGAACAGGCGCTGCGCTGGGCGCTATCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGc                                                               >  1:606041/1‑69 (MQ=255)
              gcgctgcgctGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCg                                                  >  1:2346409/1‑68 (MQ=255)
              gcgctgcgctGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCg                                                  >  1:3099864/1‑68 (MQ=255)
              gcgctgcgctGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCg                                                  >  1:1792627/1‑68 (MQ=255)
              gcgctgcgctGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGt                                                 >  1:1475179/1‑69 (MQ=255)
               cgctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACGGTGCCGGAACTGAACCGTg                                                <  1:228333/69‑1 (MQ=255)
               cgctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACGGTGCCGGAACTGAACCGTg                                                <  1:3018198/69‑1 (MQ=255)
               cgctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACGGTGCCGGAACTGAACCGTg                                                <  1:542218/69‑1 (MQ=255)
               cgctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACGGTGCCGGAACTGAACCGTg                                                <  1:785993/69‑1 (MQ=255)
               cgctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACGGTGCCGGAACTGAACCGTg                                                <  1:844526/69‑1 (MQ=255)
               cgctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACGGTGCCGGAACTGAACCGTg                                                <  1:1173043/69‑1 (MQ=255)
                             ctctctGGCGGTGAAGATTACGAGTTGTGTTTCACGGTGCCGGAACTGAACCGTg                                                <  1:2722240/55‑1 (MQ=255)
                                       gTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGac                        <  1:476271/69‑1 (MQ=255)
                                       gTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGac                        <  1:1229341/69‑1 (MQ=255)
                                       gTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGac                        <  1:1166598/69‑1 (MQ=255)
                                            gATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTg                   >  1:646651/1‑69 (MQ=255)
                                            gATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTg                   >  1:2472374/1‑69 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGct                               >  1:2891109/1‑46 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGct                               >  1:2679181/1‑46 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGct                               >  1:2517998/1‑46 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGct                               >  1:208049/1‑46 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGct                               >  1:1970357/1‑46 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGct                               >  1:1494565/1‑46 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGct                               >  1:1822058/1‑46 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGtt        >  1:207046/1‑69 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGtt        >  1:2061617/1‑69 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGtt        >  1:2591830/1‑69 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGtt        >  1:1972304/1‑69 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGtt        >  1:1376113/1‑69 (MQ=255)
                                                       gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGtt        >  1:1741565/1‑69 (MQ=255)
                                                            tCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:2914778/1‑68 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTg                   >  1:1453306/1‑53 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTAc            >  1:90703/1‑60 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTa      >  1:3010424/1‑66 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:1887161/1‑68 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:3080235/1‑68 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:2128204/1‑68 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:3212225/1‑68 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:2969589/1‑68 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:1679447/1‑68 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:581535/1‑68 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:2938747/1‑68 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:1901244/1‑68 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:1181778/1‑68 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:193466/1‑68 (MQ=255)
                                                            tCACGGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTAcc    >  1:920918/1‑68 (MQ=255)
                                                              aCTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTg  <  1:2825466/68‑1 (MQ=255)
                                                                |                                                                 
TTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTG  >  minE/264848‑264977

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