Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 384,355 A→G 40.4% L53L (CTA→CTG ybdB → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE384,3550AG40.4% 15.4 / 31.0 47L53L (CTA→CTGybdBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/12);  new base G (13/6);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGAT  >  minE/384287‑384414
                                                                    |                                                           
cttggCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtacta                                                             <  1:1060457/66‑1 (MQ=255)
cTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtacta                                                             <  1:110161/69‑1 (MQ=255)
cTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtacta                                                             <  1:1305483/69‑1 (MQ=255)
cTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtacta                                                             <  1:2161729/69‑1 (MQ=255)
cTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtacta                                                             <  1:1357877/69‑1 (MQ=255)
cTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtacta                                                             <  1:1410400/69‑1 (MQ=255)
       atgatgTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATggcg                                                      >  1:621782/1‑69 (MQ=255)
       atgatgTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATggcg                                                      >  1:3185648/1‑69 (MQ=255)
                             gCCGGTTGATACCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                <  1:25466/69‑1 (MQ=255)
                             gCCGGTTGATACCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                <  1:2352750/69‑1 (MQ=255)
                             gCCGGTTGATACCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                <  1:2796177/69‑1 (MQ=255)
                             gCCGGTTGATACCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                <  1:2895583/69‑1 (MQ=255)
                             gCCGGTTGATACCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                <  1:1624058/69‑1 (MQ=255)
                             gCCGGTTGATACCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                <  1:1294107/69‑1 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGTCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:1311923/1‑69 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                >  1:2864896/1‑59 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:164285/1‑69 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:1770295/1‑69 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:3020643/1‑69 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:2048183/1‑69 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:2932245/1‑69 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:2108540/1‑69 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:698418/1‑69 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:2179374/1‑69 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:2566593/1‑69 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:820897/1‑69 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTAGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:2052158/1‑69 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTggg                       <  1:996928/68‑1 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTggg                       <  1:3398627/68‑1 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTggg                       <  1:1144426/68‑1 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTggg                       <  1:2465067/68‑1 (MQ=255)
                                       aCCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTggg                       <  1:1911974/68‑1 (MQ=255)
                                              cTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATgg               >  1:1710759/1‑69 (MQ=255)
                                              cTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATgg               >  1:57925/1‑69 (MQ=255)
                                              cTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATgg               >  1:2070544/1‑69 (MQ=255)
                                                catcaGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGc              >  1:2233600/1‑68 (MQ=255)
                                                catcaGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGc              >  1:1276006/1‑68 (MQ=255)
                                                catcaGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGc              >  1:902590/1‑68 (MQ=255)
                                                catcaGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGc              >  1:914172/1‑68 (MQ=255)
                                                   caGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:2125408/38‑1 (MQ=255)
                                                   caGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:2799516/38‑1 (MQ=255)
                                                   caGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:1905624/38‑1 (MQ=255)
                                                   caGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:65806/38‑1 (MQ=255)
                                                   caGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:1364017/38‑1 (MQ=255)
                                                   caGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:898145/38‑1 (MQ=255)
                                                    aGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:2978114/37‑1 (MQ=255)
                                                       cGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGAttt       >  1:3247668/1‑68 (MQ=255)
                                                       cGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGAttt       >  1:918300/1‑68 (MQ=255)
                                                       cGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGAttt       >  1:610729/1‑68 (MQ=255)
                                                       cGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGAttt       >  1:1425223/1‑68 (MQ=255)
                                                           cGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTatgat  >  1:2420593/1‑69 (MQ=255)
                                                           cGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTatgat  >  1:2284065/1‑69 (MQ=255)
                                                           cGGTTTACTACATGGCGGCGCGTAGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTatgat  >  1:2196675/1‑69 (MQ=255)
                                                                    |                                                           
CTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGAT  >  minE/384287‑384414

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: