Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,353,402 G→A 43.9% T37T (ACC→ACT thiM ← hydoxyethylthiazole kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,353,4020GA43.9% 21.7 / 62.6 82T37T (ACC→ACTthiMhydoxyethylthiazole kinase
Reads supporting (aligned to +/- strand):  ref base G (30/16);  new base A (23/13);  total (53/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

CTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGT  >  minE/1353336‑1353468
                                                                  |                                                                  
ctgactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCgg                                                                   >  1:1862017/1‑68 (MQ=255)
ctgactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCgg                                                                   >  1:1231767/1‑68 (MQ=255)
ctgactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCgg                                                                   >  1:1543188/1‑68 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGtt                                                           >  1:1646095/1‑62 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:1009885/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:240395/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:2703122/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:2828705/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:2045251/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:2859858/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:3149059/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:3215336/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:3426542/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:1389332/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:100727/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:800013/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                    >  1:620525/1‑69 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacc                                                     >  1:1888943/1‑68 (MQ=255)
              cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCAGTAAAGGTTTGcacca                                                    >  1:1203496/1‑69 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:2316930/68‑1 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:1390115/68‑1 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:2215136/68‑1 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:2158267/68‑1 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:2514451/68‑1 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:2859934/68‑1 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:1830649/68‑1 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:2961698/68‑1 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:1735220/68‑1 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:430612/68‑1 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:3224116/68‑1 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:1450599/68‑1 (MQ=255)
                       tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                            <  1:3326833/68‑1 (MQ=255)
                        aaCCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGt                                          <  1:1655416/69‑1 (MQ=255)
                         aCCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGt                                          <  1:2835279/68‑1 (MQ=255)
                            aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                      <  1:588680/69‑1 (MQ=255)
                            aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                      <  1:2879520/69‑1 (MQ=255)
                            aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                      <  1:3088430/69‑1 (MQ=255)
                            aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                      <  1:579678/69‑1 (MQ=255)
                            aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                      <  1:836797/69‑1 (MQ=255)
                            aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                      <  1:966078/69‑1 (MQ=255)
                            aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                      <  1:3435788/69‑1 (MQ=255)
                            aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                      <  1:1014006/69‑1 (MQ=255)
                            aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                      <  1:2081059/69‑1 (MQ=255)
                                                        aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTg            >  1:1706504/1‑67 (MQ=255)
                                                        aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa           >  1:1175672/1‑68 (MQ=255)
                                                        aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa           >  1:2556923/1‑68 (MQ=255)
                                                        aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa           >  1:186823/1‑68 (MQ=255)
                                                        aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa           >  1:1421387/1‑68 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCCCCCCATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:3421277/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:387686/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:3018334/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:2238231/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:1733424/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:1447223/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:1253434/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:1040639/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:304361/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:1247967/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:486170/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:524486/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:1108312/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:80027/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:2382716/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:2974792/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:2139275/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:2956224/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:1957361/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:1972345/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:3015788/1‑69 (MQ=255)
                                                             ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg     >  1:2517992/1‑69 (MQ=255)
                                                               ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt    >  1:2002995/1‑68 (MQ=255)
                                                               ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt    >  1:2361298/1‑68 (MQ=255)
                                                               ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt    >  1:2484543/1‑68 (MQ=255)
                                                               ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt    >  1:2087777/1‑68 (MQ=255)
                                                               ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt    >  1:1219695/1‑68 (MQ=255)
                                                               ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt    >  1:29209/1‑68 (MQ=255)
                                                               ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt    >  1:1585709/1‑68 (MQ=255)
                                                               ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt    >  1:1736566/1‑68 (MQ=255)
                                                               ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt    >  1:3043666/1‑68 (MQ=255)
                                                               ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt    >  1:2998014/1‑68 (MQ=255)
                                                                 cGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAATAGgtgt  <  1:5116/68‑1 (MQ=255)
                                                                 cGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtgt  <  1:2085108/68‑1 (MQ=255)
                                                                 cGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtgt  <  1:2732316/68‑1 (MQ=255)
                                                                 cGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtgt  <  1:1557037/68‑1 (MQ=255)
                                                                 cGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtgt  <  1:1059610/68‑1 (MQ=255)
                                                                  |                                                                  
CTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGT  >  minE/1353336‑1353468

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: