Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,841,484 C→T 40.4% E224E (GAG→GAA ispD ← 4‑diphosphocytidyl‑2C‑methyl‑D‑erythritol synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,841,4840CT40.4% 34.1 / 54.3 79E224E (GAG→GAAispD4‑diphosphocytidyl‑2C‑methyl‑D‑erythritol synthase
Reads supporting (aligned to +/- strand):  ref base C (32/15);  new base T (23/9);  total (55/24)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

TACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCACGGCCTTCGACCAA  >  minE/1841419‑1841546
                                                                 |                                                              
tACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTTATGGATGGTTCGGGTGAGGTAAAACTCg                                                             >  1:213045/1‑69 (MQ=255)
     cAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAg                                                        <  1:3231323/69‑1 (MQ=255)
     cAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAg                                                        <  1:1642985/69‑1 (MQ=255)
      aaaaCCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAg                                                        <  1:2356042/68‑1 (MQ=255)
                 cAATTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCt                                              >  1:878949/1‑67 (MQ=255)
                 cAATTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCt                                              >  1:840747/1‑67 (MQ=255)
                 cAATTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCt                                              >  1:2474971/1‑67 (MQ=255)
                 cAATTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCt                                              >  1:2777025/1‑67 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:2964404/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:3031172/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:2823016/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:1837430/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:267424/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:1991012/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:202976/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:25846/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:2459188/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:3318173/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:379844/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:1560011/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:649195/1‑64 (MQ=255)
                   aTTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATc                                               >  1:1289065/1‑64 (MQ=255)
                            aTGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGa                                 >  1:2805568/1‑69 (MQ=255)
                            aTGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGa                                 >  1:145882/1‑69 (MQ=255)
                            aTGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGa                                 >  1:758886/1‑69 (MQ=255)
                            aTGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGAGTGa                                 >  1:630264/1‑69 (MQ=255)
                                 ttCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACttt                             <  1:3030428/68‑1 (MQ=255)
                                 ttCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACttt                             <  1:2512018/68‑1 (MQ=255)
                                 ttCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACttt                             <  1:572033/68‑1 (MQ=255)
                                 ttCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACttt                             <  1:2762430/68‑1 (MQ=255)
                                 ttCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACttt                             <  1:656536/68‑1 (MQ=255)
                                 ttCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACttt                             <  1:861917/68‑1 (MQ=255)
                                 ttCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACttt                             <  1:2313441/68‑1 (MQ=255)
                                 ttCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACttt                             <  1:1834461/68‑1 (MQ=255)
                                 ttCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACttt                             <  1:133869/68‑1 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACcttga                           >  1:1222211/1‑66 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGAt                          >  1:520910/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGAt                          >  1:328644/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGAt                          >  1:2838479/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGAt                          >  1:211163/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGAt                          >  1:1903170/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGAt                          >  1:2697066/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGAt                          >  1:2775300/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGAt                          >  1:1928250/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTatt                          >  1:1911508/1‑67 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAAt                          >  1:2754830/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAAt                          >  1:1569500/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAAt                          >  1:1032340/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAAt                          >  1:1199881/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAAt                          >  1:715024/1‑69 (MQ=255)
                                   ctcCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAAt                          >  1:2127385/1‑69 (MQ=255)
                                        gatggatgGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTTATGt                        <  1:1934872/66‑1 (MQ=255)
                                        gatggatgGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGATGt                        <  1:1118960/66‑1 (MQ=255)
                                        gatggatgGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGATGt                        <  1:3420620/66‑1 (MQ=255)
                                        gatggatgGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGATGt                        <  1:1629730/66‑1 (MQ=255)
                                        gatggatgGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGATGt                        <  1:2220846/66‑1 (MQ=255)
                                        gatggatgGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGATGt                        <  1:173657/66‑1 (MQ=255)
                                        gatggatgGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGATGt                        <  1:1785901/66‑1 (MQ=255)
                                        gatggatgGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGATGt                        <  1:652165/66‑1 (MQ=255)
                                        gatggatgGTTCGGGTGAGGTAAAATTCGGCCAGCGCCAAATCCTCCGGGCGCGTTACCTTGATGt                        <  1:695253/66‑1 (MQ=255)
                                        gatggatgGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGt                        <  1:3056518/66‑1 (MQ=255)
                                        gatggatgGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGt                        <  1:1696548/66‑1 (MQ=255)
                                        gatggatgGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGt                        <  1:1416965/66‑1 (MQ=255)
                                              tggTTCGGGTGAGGTAAACCTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCa                >  1:2091343/1‑68 (MQ=255)
                                              tggTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCa                >  1:570844/1‑68 (MQ=255)
                                              tggTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCAc               >  1:1759091/1‑69 (MQ=255)
                                              tggTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCAc               >  1:2534456/1‑69 (MQ=255)
                                              tggTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCAc               >  1:868597/1‑69 (MQ=255)
                                                  tCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATg                         >  1:323935/1‑55 (MQ=255)
                                                  tCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGtt                       >  1:2826722/1‑57 (MQ=255)
                                                  tCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGtt                       >  1:2394605/1‑57 (MQ=255)
                                                  tCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGtt                       >  1:1269606/1‑57 (MQ=255)
                                                             aaaaCTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCACg              >  1:1609878/1‑55 (MQ=255)
                                                             aaaaCTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCACGGCCTTCGACCaa  >  1:3092610/1‑67 (MQ=255)
                                                             aaaaCTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCACGGCCTTCGACCaa  >  1:824110/1‑67 (MQ=255)
                                                             aaaaCTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCACGGCCTTCGACCaa  >  1:1534892/1‑67 (MQ=255)
                                                             aaaaCTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCACGGCCTTCGACCaa  >  1:398428/1‑67 (MQ=255)
                                                             aaaaCTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCACGGCCTTCGACCaa  >  1:2153879/1‑67 (MQ=255)
                                                             aaaaCTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCACGGCCGTCGACCaa  >  1:2112882/1‑67 (MQ=255)
                                                                 |                                                              
TACGTCAAAACCGTGTCCAATTCGCATTATGTATTCTCCTGATGGATGGTTCGGGTGAGGTAAAACTCGGCCAGTGCCAAATCTTCCGGGCGCGTGACTTTAATGTTATCCGCACGGCCTTCGACCAA  >  minE/1841419‑1841546

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: