Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,222,098 A→T 20.3% T300T (ACA→ACT accC → acetyl‑CoA carboxylase, biotin carboxylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,222,0980AT20.3% 105.9 / 21.1 59T300T (ACA→ACTaccCacetyl‑CoA carboxylase, biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base T (7/5);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

GTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAACCGCTGTC  >  minE/2222032‑2222163
                                                                  |                                                                 
gTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGa                                                                 <  1:1126292/69‑1 (MQ=255)
gTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGa                                                                 <  1:495871/69‑1 (MQ=255)
gTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGa                                                                 <  1:3315824/69‑1 (MQ=255)
gTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGa                                                                 <  1:2542474/69‑1 (MQ=255)
gTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGa                                                                 <  1:2499430/69‑1 (MQ=255)
gTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGa                                                                 <  1:1392706/69‑1 (MQ=255)
gTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGa                                                                 <  1:2341364/69‑1 (MQ=255)
gTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGa                                                                 <  1:2220227/69‑1 (MQ=255)
gTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGa                                                                 <  1:2012703/69‑1 (MQ=255)
 ttCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGa                                                                 <  1:98797/68‑1 (MQ=255)
     aaaaCGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             <  1:1630676/68‑1 (MQ=255)
     aaaaCGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             <  1:1676233/68‑1 (MQ=255)
     aaaaCGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             <  1:965375/68‑1 (MQ=255)
     aaaaCGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             <  1:909010/68‑1 (MQ=255)
     aaaaCGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             <  1:2695945/68‑1 (MQ=255)
     aaaaCGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             <  1:1877511/68‑1 (MQ=255)
     aaaaCGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             <  1:2999803/68‑1 (MQ=255)
     aaaaCGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             <  1:3075110/68‑1 (MQ=255)
     aaaaCGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             <  1:397023/68‑1 (MQ=255)
        aCGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             >  1:986480/1‑65 (MQ=255)
        aCGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             >  1:3134974/1‑65 (MQ=255)
            cgtgTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             <  1:3426152/58‑1 (MQ=255)
                ttCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGtt                                                 >  1:462951/1‑69 (MQ=255)
                ttCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGtt                                                 >  1:884456/1‑69 (MQ=255)
                 tCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTg                                                >  1:3410777/1‑69 (MQ=255)
                  cTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             >  1:668650/1‑55 (MQ=255)
                  cTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             >  1:85360/1‑55 (MQ=255)
                  cTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                                                             >  1:875355/1‑55 (MQ=255)
                  cTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGa                                               >  1:1543685/1‑69 (MQ=255)
                  cTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGa                                               >  1:2292664/1‑69 (MQ=255)
                       tCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACTGAAATGATCACCGGCGTTGAc                                              <  1:1593493/65‑1 (MQ=255)
                       tCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACTGAAATGATCACCGGCGTTGAc                                              <  1:2628246/65‑1 (MQ=255)
                       tCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACTGAAATGATCACCGGCGTTGAc                                              <  1:1218052/65‑1 (MQ=255)
                       tCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACTGAAATGATCACCGGCGTTGAc                                              <  1:1544815/65‑1 (MQ=255)
                       tCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACTGAAATGATCACCGGCGTTGAc                                              <  1:1479699/65‑1 (MQ=255)
                             aaaTGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCAcc                                                       <  1:619525/50‑1 (MQ=255)
                                            ttCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCg                     >  1:1806315/1‑69 (MQ=255)
                                            ttCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCg                     >  1:2144010/1‑69 (MQ=255)
                                            ttCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCg                     >  1:215931/1‑69 (MQ=255)
                                                   aGAACACCCGGTTACTGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGc                 >  1:403797/1‑66 (MQ=255)
                                                   aGAACACCCGGTTACTGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGc                 >  1:1580953/1‑66 (MQ=255)
                                                   aGAACACCCGGTTACTGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGc                 >  1:1763119/1‑66 (MQ=255)
                                                   aGAACACCCGGTTACTGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGc                 >  1:624258/1‑66 (MQ=255)
                                                   aGAACACCCGGTTACTGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGc                 >  1:2160413/1‑66 (MQ=255)
                                                   aGAACACCCGGTTACTGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGc                 >  1:2275258/1‑66 (MQ=255)
                                                   aGAACACCCGGTTACTGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGc                 >  1:2283724/1‑66 (MQ=255)
                                                        acCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAAc         <  1:3359791/69‑1 (MQ=255)
                                                        acCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAAc         <  1:113156/69‑1 (MQ=255)
                                                        acCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAAc         <  1:8640/69‑1 (MQ=255)
                                                        acCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAAc         <  1:2761091/69‑1 (MQ=255)
                                                               tACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAACCGCTGTc  >  1:2578053/1‑69 (MQ=255)
                                                               tACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAACCGCTGTc  >  1:1175655/1‑69 (MQ=255)
                                                               tACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAACCGCTGTc  >  1:2081371/1‑69 (MQ=255)
                                                               tACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAACCGCTGTc  >  1:2420004/1‑69 (MQ=255)
                                                               tACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAACCGCTGTc  >  1:1885500/1‑69 (MQ=255)
                                                               tACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAACCGCTGTc  >  1:2730798/1‑69 (MQ=255)
                                                               tACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAACCGCTGTc  >  1:173666/1‑69 (MQ=255)
                                                               tACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAACCGCTGTc  >  1:1296764/1‑69 (MQ=255)
                                                               tACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAACCGCTGTc  >  1:1123217/1‑69 (MQ=255)
                                                                  |                                                                 
GTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGAACAGCTGCGTATCGCTGCCGGTCAACCGCTGTC  >  minE/2222032‑2222163

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: