Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,613,749 G→T 73.1% H158N (CAC→AAC)  acpT ← holo‑(acyl carrier protein) synthase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,613,7490GT73.1% 107.4 / 50.7 89H158N (CAC→AAC) acpTholo‑(acyl carrier protein) synthase 2
Reads supporting (aligned to +/- strand):  ref base G (11/13);  new base T (17/48);  total (28/61)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCAGATCGCAAGGCTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGCTTC  >  minE/2613691‑2613812
                                                          |                                                               
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGATGAGTTAt                                                               <  1:873234/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:1878000/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:243486/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:2942131/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:3221201/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:767860/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:737426/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:2514335/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:554677/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:543717/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:1686133/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:165484/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:411761/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:3329023/61‑1 (MQ=255)
gTGCAGATCGCCAGGCTTAAATTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:2784553/61‑1 (MQ=255)
       tCGCAAGGCTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGa                                                                <  1:1393184/53‑1 (MQ=255)
       tCGCAAGGCTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGa                                                                <  1:2690050/53‑1 (MQ=255)
       tCGCAAGGCTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGa                                                                <  1:2081500/53‑1 (MQ=255)
        cGCAAGGCTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCt                                               >  1:3259192/1‑69 (MQ=255)
        cGCAAGGCTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCt                                               >  1:1170432/1‑69 (MQ=255)
        cGCAAGGCTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCt                                               >  1:818802/1‑69 (MQ=255)
        cGCAAGGCTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCt                                               >  1:2288252/1‑69 (MQ=255)
              gCTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGAt                                         >  1:1007216/1‑69 (MQ=255)
              gCTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGAt                                         >  1:431886/1‑69 (MQ=255)
               cTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGAtt                                        >  1:1963073/1‑69 (MQ=255)
               cTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGAtt                                        >  1:1506480/1‑69 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:3005810/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:3100770/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:2854121/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:2839694/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:3344430/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:275996/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:634831/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:831429/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:975973/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:2274556/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:2168807/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:1079986/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:1143681/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:1430063/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:1769591/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:13001/1‑67 (MQ=255)
                  aaaTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGAttt                                       >  1:1337807/1‑67 (MQ=255)
                    aTTTTCTAGCTGACAATGGCTGACTGACAGCGAGGAGTTAt                                                               <  1:2164811/41‑1 (MQ=255)
                                    tGGCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACg                    >  1:1369513/1‑68 (MQ=255)
                                    tGGCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGc                   >  1:909418/1‑69 (MQ=255)
                                     ggCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCt                  <  1:1830901/69‑1 (MQ=255)
                                     ggCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCt                  <  1:376185/69‑1 (MQ=255)
                                     ggCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCt                  <  1:3066082/69‑1 (MQ=255)
                                     ggCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCt                  <  1:3065751/69‑1 (MQ=255)
                                     ggCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCt                  <  1:1832917/69‑1 (MQ=255)
                                     ggCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCt                  <  1:1912295/69‑1 (MQ=255)
                                     ggCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCt                  <  1:1952823/69‑1 (MQ=255)
                                     ggCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCt                  <  1:279002/69‑1 (MQ=255)
                                     ggCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCt                  <  1:1286203/69‑1 (MQ=255)
                                     ggCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCt                  <  1:815039/69‑1 (MQ=255)
                                      gCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTg                 >  1:879636/1‑69 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:3370202/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:1627164/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:1077460/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:331157/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:539314/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:719111/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:1111738/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:2207075/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:244357/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:2094484/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:2611444/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:1770549/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:274432/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:2021395/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:2948230/69‑1 (MQ=255)
                                             gacAGCGAGGAGTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACg          <  1:1881450/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCTCCACGCTGTTTAACGATGGcttc  <  1:2262159/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCATGCTGTTTAACGATGGcttc  <  1:3275711/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:482464/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:3081745/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:692466/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:2035828/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:2622257/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:1188928/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:3268271/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:1645777/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:2383412/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:2172589/69‑1 (MQ=255)
                                                     ggagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:3306709/69‑1 (MQ=255)
                                                      gagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:1450942/68‑1 (MQ=255)
                                                      gagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:2026728/68‑1 (MQ=255)
                                                      gagTTATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGcttc  <  1:2339899/68‑1 (MQ=255)
                                                          |                                                               
GTGCAGATCGCAAGGCTTAAATTTTCAAGCTGACAATGGCTGACTGACAGCGAGGAGTGATAGGTGCTGTCTACGCTGACGATTTGCCAGGCGCTGCCGCCACGCTGTTTAACGATGGCTTC  >  minE/2613691‑2613812

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