Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 889,796 T→G 47.9% A66A (GCA→GCC fabI ← enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE889,7960TG47.9% ‑3.3 / 85.3 69A66A (GCA→GCCfabIenoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base T (19/17);  new base G (25/8);  total (44/25)
Fisher's exact test for biased strand distribution p-value = 7.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCGTCAAATTTCGGCCAAACTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAATTCTTCTACGCGGCCT  >  minE/889732‑889859
                                                                |                                                               
aaCCGTCAAATTTCGGCCAAACTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCaa                                                             <  1:869366/69‑1 (MQ=255)
 aCCGTCAAATTTCGGCCAAACTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAAc                                                            <  1:2917772/69‑1 (MQ=255)
 aCCGTCAAATTTCGGCCAAACTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAAc                                                            <  1:464317/69‑1 (MQ=255)
       aaaTTTCGGCCAAACTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCa                                                      <  1:1006793/69‑1 (MQ=255)
       aaaTTTCGGCCAAACTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCa                                                      <  1:2110659/69‑1 (MQ=255)
       aaaTTTCGGCCAAACTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCa                                                      <  1:352618/69‑1 (MQ=255)
       aaaTTTCGGCCAAACTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCa                                                      <  1:2016248/69‑1 (MQ=255)
       aaaTTTCGGCCAAACTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCa                                                      <  1:1896667/69‑1 (MQ=255)
       aaaTTTCGGCCAAACTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCa                                                      <  1:116362/69‑1 (MQ=255)
                  aaaCTTTCCCCAGTTCATCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:3130281/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGATCATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:1522266/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACAt                                                          >  1:1242833/1‑54 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACt                                                    >  1:276150/1‑60 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTTCAGAACGa                                           >  1:2563975/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:561661/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:3050890/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:2850506/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:2953023/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           <  1:2363410/69‑1 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:2410700/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:2831509/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           <  1:242215/69‑1 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:2807076/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:3256744/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:2636983/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           <  1:2671601/69‑1 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:595883/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           <  1:1236764/69‑1 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           <  1:3330384/69‑1 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           <  1:3353945/69‑1 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:164207/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:867180/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           <  1:1687085/69‑1 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:1708393/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           <  1:1785693/69‑1 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:1871937/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:831865/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:1995353/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:2110739/1‑69 (MQ=255)
                  aaaCTTTCCCCAGTTCAGCGAACATGGTGTCGATGCAGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:1614268/1‑69 (MQ=255)
                  aaaCTTTCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGAt                                          >  1:1995042/1‑69 (MQ=255)
                  aaaCTTACCCCAGTTCAGCGAACATGGTGTCGATGCTGGTATCTTCGGCAACATCGCACTGCAGAACGa                                           >  1:4953/1‑69 (MQ=255)
                     cTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGAt                                          >  1:3108377/1‑67 (MQ=255)
                     cTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGAt                                          >  1:2936172/1‑67 (MQ=255)
                     cTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGAt                                          >  1:2325537/1‑67 (MQ=255)
                     cTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGAt                                          >  1:2874749/1‑67 (MQ=255)
                     cTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGAt                                          >  1:2501789/1‑67 (MQ=255)
                     cTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGAt                                          >  1:2261933/1‑67 (MQ=255)
                     cTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGAt                                          >  1:2251660/1‑67 (MQ=255)
                     cTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGAt                                          >  1:2073323/1‑67 (MQ=255)
                     cTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGAt                                          >  1:1597131/1‑67 (MQ=255)
                     cTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGAt                                          >  1:12453/1‑67 (MQ=255)
                                tCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACAt                                                          >  1:3023272/1‑40 (MQ=255)
                                tCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCATAACGATGTCAGAACCCaa                              >  1:997852/1‑68 (MQ=255)
                                tCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCaa                              >  1:1466271/1‑68 (MQ=255)
                                 cAGCGAACATGGTTTCGATGCTGTGATCTTCTGCAACATGGCa                                                      <  1:2899841/43‑1 (MQ=25)
                                                 gATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAAttcttc            >  1:2601049/1‑69 (MQ=255)
                                                 gATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAAttcttc            >  1:3435165/1‑69 (MQ=255)
                                                 gATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAAttcttc            >  1:3226571/1‑69 (MQ=255)
                                                 gATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAAttcttc            >  1:296156/1‑69 (MQ=255)
                                                 gATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAAttcttc            >  1:2474899/1‑69 (MQ=255)
                                                 gATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAAttcttc            >  1:2078304/1‑69 (MQ=255)
                                                          aTCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAATTCTTCTACGCGGc    <  1:3170132/68‑1 (MQ=255)
                                                          aTCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAATTCTTCTACGCGGc    <  1:3154693/68‑1 (MQ=255)
                                                          aTCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAATTCTTCTACGCGGc    <  1:4855/68‑1 (MQ=255)
                                                          aTCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAATTCTTCTACGCGGc    <  1:2124527/68‑1 (MQ=255)
                                                          aTCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAATTCTTCTACGCGGc    <  1:797009/68‑1 (MQ=255)
                                                          aTCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAATTCTTCTACGCGGc    <  1:1936643/68‑1 (MQ=255)
                                                          aTCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAATTCTTCTACGCGGc    <  1:1054898/68‑1 (MQ=255)
                                                             ttctGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCaaa                  <  1:2952464/51‑1 (MQ=255)
                                                             ttctGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAATTCTTCTACGCGGCCt  >  1:1666301/1‑67 (MQ=255)
                                                             ttctGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAATTCTTCTACGCGGCCt  >  1:1083064/1‑67 (MQ=255)
                                                                |                                                               
AACCGTCAAATTTCGGCCAAACTTTCCCCAGTTCAGCGAACATGGTGTCGATGCTGGCATCTTCTGCAACATCGCACTGCAGAACGATGTCAGAACCCAATTGAGCGGCAAATTCTTCTACGCGGCCT  >  minE/889732‑889859

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