Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,353,423 A→G 43.2% N30N (AAT→AAC thiM ← hydoxyethylthiazole kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,353,4230AG43.2% 19.2 / 62.7 81N30N (AAT→AACthiMhydoxyethylthiazole kinase
Reads supporting (aligned to +/- strand):  ref base A (29/17);  new base G (22/13);  total (51/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCG  >  minE/1353359‑1353474
                                                                |                                                   
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:1390115/68‑1 (MQ=255)
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:3224116/68‑1 (MQ=255)
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:3326833/68‑1 (MQ=255)
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:1735220/68‑1 (MQ=255)
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:1450599/68‑1 (MQ=255)
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:1830649/68‑1 (MQ=255)
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:2961698/68‑1 (MQ=255)
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:430612/68‑1 (MQ=255)
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:2859934/68‑1 (MQ=255)
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:2158267/68‑1 (MQ=255)
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:2215136/68‑1 (MQ=255)
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:2316930/68‑1 (MQ=255)
tAACCATCGCTGGCGATGCGCCGAGCGCCAGCAAGGTATTGGCAGTAAAAGTTTGCACCACATCGTTg                                                  <  1:2514451/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:1014006/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:966078/69‑1 (MQ=255)
     aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                            <  1:836797/69‑1 (MQ=255)
     aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                            <  1:588680/69‑1 (MQ=255)
     aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATg                                            <  1:579678/69‑1 (MQ=255)
                                 aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTg                  >  1:1706504/1‑67 (MQ=255)
                                 aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa                 >  1:186823/1‑68 (MQ=255)
                                 aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa                 >  1:2556923/1‑68 (MQ=255)
                                 aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa                 >  1:1421387/1‑68 (MQ=255)
                                 aGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa                 >  1:1175672/1‑68 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCCCCCCATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:3421277/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:2382716/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:2517992/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:80027/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:524486/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:2956224/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:486170/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:2974792/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:3015788/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:3018334/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:304361/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:387686/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1040639/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1108312/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:2238231/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1247967/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1253434/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1447223/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:2139275/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1733424/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1957361/1‑69 (MQ=255)
                                      ttGGCAGTAAAAGTTTGCACCACATCGTTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAgg           >  1:1972345/1‑69 (MQ=255)
                                        ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:1736566/1‑68 (MQ=255)
                                        ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:1219695/1‑68 (MQ=255)
                                        ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:1585709/1‑68 (MQ=255)
                                        ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:2484543/1‑68 (MQ=255)
                                        ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:3043666/1‑68 (MQ=255)
                                        ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:2002995/1‑68 (MQ=255)
                                        ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:2998014/1‑68 (MQ=255)
                                        ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:2361298/1‑68 (MQ=255)
                                        ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:29209/1‑68 (MQ=255)
                                        ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgt          >  1:2087777/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:207942/1‑69 (MQ=255)
                                              aaaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgc   >  1:2348966/1‑69 (MQ=255)
                                              aaaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgc   >  1:785950/1‑69 (MQ=255)
                                              aaaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgc   >  1:702387/1‑69 (MQ=255)
                                              aaaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgc   >  1:277925/1‑69 (MQ=255)
                                              aaaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcg    >  1:2012010/1‑68 (MQ=255)
                                               aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgcg  >  1:2687504/1‑69 (MQ=255)
                                               aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgcg  >  1:1392179/1‑69 (MQ=255)
                                               aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgcg  >  1:1670932/1‑69 (MQ=255)
                                               aaGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgcg  >  1:3300576/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)
                                                                |                                                   
TAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCG  >  minE/1353359‑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: