Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,250,566 C→T 54.5% I16I (ATC→ATT thiH → thiamin biosynthesis ThiGH complex subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,250,5660CT54.5% 3.8 / 54.7 66I16I (ATC→ATTthiHthiamin biosynthesis ThiGH complex subunit
Reads supporting (aligned to +/- strand):  ref base C (23/7);  new base T (24/12);  total (47/19)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

GGCATCGGCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACT  >  minE/2250510‑2250628
                                                        |                                                              
ggCATCGGCATGAAAACCTTCGGCGATCGCTGGCGACAACTGGACTGGGACGACATTCGCCTGCGTATc                                                    <  1:232378/69‑1 (MQ=255)
ggCATCGGCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATTCGCCTGCGTATc                                                    <  1:991789/69‑1 (MQ=255)
ggCATCGGCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATTCGCCTGCGTATc                                                    <  1:2460125/69‑1 (MQ=255)
ggCATCGGCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATTCGCCTGCGTATc                                                    <  1:1749148/69‑1 (MQ=255)
ggCATCGGCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATTCGCCTGCGTATc                                                    <  1:1700952/69‑1 (MQ=255)
ggCATCGGCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATTCGCCTGCGTATc                                                    <  1:3388591/69‑1 (MQ=255)
ggCATCGGCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATTCGCCTGCGTATc                                                    <  1:47745/69‑1 (MQ=255)
ggCATCGGCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATTCGCCTGCGTATc                                                    <  1:560113/69‑1 (MQ=255)
ggCATCGGCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATTCGCCTGCGTATc                                                    <  1:139545/69‑1 (MQ=255)
ggCATCGGCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATTCGCCTGCGTATc                                                    <  1:946597/69‑1 (MQ=255)
  cATCGGCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATTCGCCTGCGTATc                                                    <  1:2919074/67‑1 (MQ=255)
      ggCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACgg                                               >  1:1437042/1‑68 (MQ=255)
      ggCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACgg                                               >  1:2793335/1‑68 (MQ=255)
      ggCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACgg                                               >  1:2714004/1‑68 (MQ=255)
      ggCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACg                                                >  1:3339256/1‑67 (MQ=255)
      ggCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACg                                                >  1:3335022/1‑67 (MQ=255)
            aaaaCCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATcc                                                               <  1:3404456/46‑1 (MQ=255)
            aaaaCCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATcc                                                               <  1:2980032/46‑1 (MQ=255)
            aaaaCCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATc                                                    <  1:2733630/57‑1 (MQ=255)
            aaaaCCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATc                                                    <  1:2675802/57‑1 (MQ=255)
            aaaaCCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATc                                                    <  1:2826217/57‑1 (MQ=255)
                           cgatGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGAGGTagag                         >  1:1683500/4‑69 (MQ=255)
                           cGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGCCGTagc                          >  1:369539/1‑67 (MQ=255)
                           cGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTagag                         >  1:1591575/1‑69 (MQ=255)
                           cGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTagag                         >  1:3309131/1‑69 (MQ=255)
                           cGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTagag                         >  1:1942964/1‑69 (MQ=255)
                           cGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTagag                         >  1:724567/1‑69 (MQ=255)
                           cGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTagag                         >  1:1687910/1‑69 (MQ=255)
                           cGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTagag                         >  1:927874/1‑69 (MQ=255)
                               ggCTACAACTGGACTGGGACGACATCCGCCTGCGTATc                                                    <  1:3300125/38‑1 (MQ=255)
                               ggCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCggg                     <  1:3104559/69‑1 (MQ=255)
                                  gACAACTGGACTGGCACGACATTCGCCTGCGTATc                                                    <  1:2183571/35‑1 (MQ=38)
                                    cAACTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTa                <  1:1745478/69‑1 (MQ=255)
                                     aaCTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTa                <  1:1511196/68‑1 (MQ=255)
                                        tggactggGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATg            >  1:149274/1‑69 (MQ=255)
                                               ggACGACATCCGCCTGCGTTTCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCa     >  1:3208241/1‑69 (MQ=255)
                                               ggACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCa     >  1:1722271/1‑69 (MQ=255)
                                               ggACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCa     >  1:922033/1‑69 (MQ=255)
                                               ggACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCa     >  1:1937857/1‑69 (MQ=255)
                                               ggACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCa     >  1:2706997/1‑69 (MQ=255)
                                               ggACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCa     >  1:38993/1‑69 (MQ=255)
                                               ggACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCa     >  1:2449932/1‑69 (MQ=255)
                                               ggACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCa     >  1:3328447/1‑69 (MQ=255)
                                               ggACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCa     >  1:2385118/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGgcg                   >  1:2282388/1‑52 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAAc   >  1:3037849/1‑68 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:1096384/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:553437/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:976883/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:1172123/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:577495/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:642352/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:648636/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:93379/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:910262/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:1544848/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:165974/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:3282443/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:3242929/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:3202005/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:2943952/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:1004737/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:1827721/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:2135487/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:2418421/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:2152931/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACt  >  1:2322237/1‑69 (MQ=255)
                                                  cgacATTCGCCTGCGTATCAACAGCAAAACGGCTGCTGACGTAGAGCCGGCGCTAAATGCCTCGCAACt  >  1:3032552/1‑69 (MQ=255)
                                                        |                                                              
GGCATCGGCATGAAAACCTTCAGCGATCGCTGGCGACAACTGGACTGGGACGACATCCGCCTGCGTATCAACGGCAAAACGGCTGCTGACGTAGAGCGGGCGCTAAATGCCTCGCAACT  >  minE/2250510‑2250628

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: