Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,841,416 A→G 40.3% H10H (CAT→CAC ispF ← 2C‑methyl‑D‑erythritol 2,4‑cyclodiphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,841,4160AG40.3% 46.8 / 80.8 114H10H (CAT→CACispF2C‑methyl‑D‑erythritol 2,4‑cyclodiphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (51/17);  new base G (36/10);  total (87/27)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

AATCCTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAA  >  minE/1841354‑1841481
                                                              |                                                                 
aaTCCTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTAc                                                              >  1:1998791/1‑68 (MQ=255)
aaTCCTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTAc                                                              >  1:2962663/1‑68 (MQ=255)
aaTCCTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTAc                                                              >  1:2069536/1‑68 (MQ=255)
   ccTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGt                                                            >  1:1625241/1‑67 (MQ=255)
   ccTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTc                                                           >  1:513510/1‑68 (MQ=255)
   ccTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTc                                                           >  1:2469650/1‑68 (MQ=255)
   ccTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTc                                                           >  1:1927777/1‑68 (MQ=255)
   ccTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTc                                                           >  1:2259945/1‑68 (MQ=255)
   ccTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTc                                                           >  1:2579784/1‑68 (MQ=255)
   ccTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTc                                                           >  1:3212695/1‑68 (MQ=255)
   ccTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTc                                                           >  1:1068656/1‑68 (MQ=255)
        ttCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAAcc                                                     >  1:1236328/1‑69 (MQ=255)
        ttCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAAcc                                                     >  1:709233/1‑69 (MQ=255)
        ttCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAAcc                                                     >  1:2315977/1‑69 (MQ=255)
        ttCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAAcc                                                     >  1:1481950/1‑69 (MQ=255)
        ttCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAAcc                                                     >  1:1843217/1‑69 (MQ=255)
        ttCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAAcc                                                     >  1:2824566/1‑69 (MQ=255)
        ttCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAAc                                                      >  1:2931120/1‑68 (MQ=255)
                   tGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAAtt                                           >  1:1826023/1‑68 (MQ=255)
                   tGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAAtt                                           >  1:400634/1‑68 (MQ=255)
                   tGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAAtt                                           >  1:3155079/1‑68 (MQ=255)
                   tGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAAtt                                           >  1:110879/1‑68 (MQ=255)
                   tGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAAtt                                           >  1:1098323/1‑68 (MQ=255)
                   tGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAAtt                                           >  1:1066940/1‑68 (MQ=255)
                   agCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAAtt                                           >  1:656031/2‑68 (MQ=255)
                    gCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCg                                         >  1:1848008/1‑69 (MQ=255)
                    gCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCg                                         >  1:175550/1‑69 (MQ=255)
                    gCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCg                                         >  1:763268/1‑69 (MQ=255)
                    gCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCg                                         >  1:3101190/1‑69 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:424418/69‑1 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:420436/69‑1 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:410906/69‑1 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:2373069/69‑1 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:2940818/69‑1 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:342167/69‑1 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:3247672/69‑1 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:2186285/69‑1 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:685286/69‑1 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:1374248/69‑1 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:943144/69‑1 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:1130097/69‑1 (MQ=255)
                       tACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCAt                                      <  1:2026098/69‑1 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAAcc                                                     >  1:2017572/1‑51 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCaa                                             >  1:3184086/1‑59 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCa                                              >  1:942512/1‑58 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCa                                       >  1:1701847/1‑65 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:333817/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:3216929/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:3203449/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:1237564/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:3018998/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:378544/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:1240149/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:3336571/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:2937575/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:878895/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:83877/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:2407419/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:2424792/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:2441422/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:1975654/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCAt                                      >  1:1917304/1‑66 (MQ=255)
                          gCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTAGCAt                                      >  1:2516168/1‑66 (MQ=255)
                                       aTTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTga                       >  1:1783897/1‑68 (MQ=255)
                                       aTTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTga                       >  1:864127/1‑68 (MQ=255)
                                       aTTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTg                        >  1:2833578/1‑67 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATg                                  >  1:832747/1‑51 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTgatggatg                 >  1:1843739/1‑68 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTgatggatg                 >  1:1238512/1‑68 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTgatggatg                 >  1:473463/1‑68 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTgatggatg                 >  1:3376052/1‑68 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGatggatgg                >  1:1060907/1‑69 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGatggatgg                >  1:612779/1‑69 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGatggatgg                >  1:966233/1‑69 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGatggatgg                >  1:2339973/1‑69 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGatggatgg                >  1:1103151/1‑69 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGatggatgg                >  1:1700727/1‑69 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGatggatgg                >  1:2010209/1‑69 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGatggatgg                >  1:2608249/1‑69 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGatggatgg                >  1:2644691/1‑69 (MQ=255)
                                             ccTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGatggatgg                >  1:295480/1‑69 (MQ=255)
                                                    cGCCAAATGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCggg          <  1:646117/68‑1 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCggg          <  1:405430/68‑1 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCggg          <  1:2921337/68‑1 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCggg          <  1:2906338/68‑1 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCggg          <  1:2631083/68‑1 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCggg          <  1:369058/68‑1 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCggg          <  1:1142450/68‑1 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCggg          <  1:2104164/68‑1 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCggg          <  1:2036288/68‑1 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCgg           >  1:3349100/1‑67 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCgg           >  1:138177/1‑67 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCgg           >  1:156380/1‑67 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCgg           >  1:3262624/1‑67 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCgg           >  1:3057240/1‑67 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCgg           >  1:2170253/1‑67 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCgg           >  1:2193460/1‑67 (MQ=255)
                                                    cGCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTAGCATTATGTATTCTCCTGATGGATGGTTCgg           >  1:1340719/1‑67 (MQ=255)
                                                     gCCAAAGGCGTGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCggg          <  1:3318343/67‑1 (MQ=255)
                                                           ggCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTaa  <  1:3007635/69‑1 (MQ=255)
                                                              |                                                                 
AATCCTTTTTCGTAAGGAATGCGTACGCCACCAATGATAATTGGGCCTTCACCGCCAAAGGCATGTACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAA  >  minE/1841354‑1841481

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: