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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,353,4080GA36.0% 66.3 / 56.7 97T35T (ACC→ACTthiMhydoxyethylthiazole kinase
Reads supporting (aligned to +/- strand):  ref base G (45/17);  new base A (22/13);  total (67/30)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

CGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCG  >  minE/1353350‑1353474
                                                          |                                                                  
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGtt                                                                 >  1:1646095/1‑62 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:1009885/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:800013/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:620525/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:3426542/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:3215336/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:3149059/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:2859858/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:2828705/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:2703122/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:240395/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:2045251/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:1389332/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                          >  1:100727/1‑69 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacc                                                           >  1:1888943/1‑68 (MQ=255)
cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCAGTAAAGGTTTGcacca                                                          >  1:1203496/1‑69 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:1390115/68‑1 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:430612/68‑1 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:3326833/68‑1 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:3224116/68‑1 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:2961698/68‑1 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:2859934/68‑1 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:2514451/68‑1 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:2316930/68‑1 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:2215136/68‑1 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:2158267/68‑1 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:1830649/68‑1 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:1735220/68‑1 (MQ=255)
         tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:1450599/68‑1 (MQ=255)
          aaCCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGt                                                <  1:1655416/69‑1 (MQ=255)
           aCCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGt                                                <  1:2835279/68‑1 (MQ=255)
              aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                            <  1:2081059/69‑1 (MQ=255)
              aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                            <  1:2879520/69‑1 (MQ=255)
              aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                            <  1:3088430/69‑1 (MQ=255)
              aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                            <  1:3435788/69‑1 (MQ=255)
              aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                            <  1:579678/69‑1 (MQ=255)
              aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                            <  1:588680/69‑1 (MQ=255)
              aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                            <  1:1014006/69‑1 (MQ=255)
              aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                            <  1:836797/69‑1 (MQ=255)
              aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                            <  1:966078/69‑1 (MQ=255)
                                          aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTg                  >  1:1706504/1‑67 (MQ=255)
                                          aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa                 >  1:1421387/1‑68 (MQ=255)
                                          aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa                 >  1:186823/1‑68 (MQ=255)
                                          aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa                 >  1:1175672/1‑68 (MQ=255)
                                          aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa                 >  1:2556923/1‑68 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCCCCCCATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:3421277/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:387686/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:304361/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:3018334/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:3015788/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:2974792/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1108312/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:2956224/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1447223/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:2139275/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1040639/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1733424/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:1253434/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:2238231/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:2517992/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1972345/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:80027/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1957361/1‑69 (MQ=255)
                                               ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:2382716/1‑69 (MQ=255)
                                                 ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:3043666/1‑68 (MQ=255)
                                                 ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:2002995/1‑68 (MQ=255)
                                                 ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:2087777/1‑68 (MQ=255)
                                                 ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:1736566/1‑68 (MQ=255)
                                                 ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:2998014/1‑68 (MQ=255)
                                                 ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:2361298/1‑68 (MQ=255)
                                                 ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:1219695/1‑68 (MQ=255)
                                                 ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:2484543/1‑68 (MQ=255)
                                                 ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:1585709/1‑68 (MQ=255)
                                                 ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:29209/1‑68 (MQ=255)
                                                   cGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAATAGgtgt        <  1:5116/68‑1 (MQ=255)
                                                   cGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtgt        <  1:2085108/68‑1 (MQ=255)
                                                   cGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtgt        <  1:1059610/68‑1 (MQ=255)
                                                   cGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtgt        <  1:1557037/68‑1 (MQ=255)
                                                   cGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtgt        <  1:2732316/68‑1 (MQ=255)
                                                     gTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtgt        <  1:2204905/66‑1 (MQ=255)
                                                       aaaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgc   >  1:2348966/1‑69 (MQ=255)
                                                       aaaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgc   >  1:207942/1‑69 (MQ=255)
                                                       aaaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgc   >  1:277925/1‑69 (MQ=255)
                                                       aaaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgc   >  1:785950/1‑69 (MQ=255)
                                                       aaaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgc   >  1:702387/1‑69 (MQ=255)
                                                       aaaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcg    >  1:2012010/1‑68 (MQ=255)
                                                        aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgcg  >  1:1670932/1‑69 (MQ=255)
                                                        aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgcg  >  1:2687504/1‑69 (MQ=255)
                                                        aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgcg  >  1:3300576/1‑69 (MQ=255)
                                                        aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgcg  >  1:1392179/1‑69 (MQ=255)
                                                        aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgcg  >  1:3111185/1‑69 (MQ=255)
                                                        aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgcg  >  1:965142/1‑69 (MQ=255)
                                                        aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgcg  >  1:2953350/1‑69 (MQ=255)
                                                        aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgc   >  1:1187357/1‑68 (MQ=255)
                                                          |                                                                  
CGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCG  >  minE/1353350‑1353474

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: