Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 394,053 T→C 7.0% V71A (GTA→GCA)  ahpC → alkyl hydroperoxide reductase, C22 subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE394,0530TC 7.0% 215.7 / 6.1 71V71A (GTA→GCA) ahpCalkyl hydroperoxide reductase, C22 subunit
Reads supporting (aligned to +/- strand):  ref base T (23/43);  new base C (2/3);  total (25/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

CCGACCGAACTGGGTGACGTTGCTGACCACTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTAAAATCAAA  >  minE/393983‑394120
                                                                      |                                                                   
ccgaccgaACTGGGTGACGTTGCTGACCACTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGt                                                                     <  1:2693110/71‑1 (MQ=255)
ccgaccgaACTGGGTGACGTTGCTGACCACTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGt                                                                     <  1:503194/71‑1 (MQ=255)
       aaCTGGGTGACGTTGCTGACCACTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTAcc                                                              <  1:2077611/71‑1 (MQ=255)
        aCTGGGTGACGTTGCTGACCACTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTAcc                                                              >  1:2485556/1‑70 (MQ=255)
        aCTGGGTGACGTTGCTGACCACTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTAcc                                                              >  1:1620073/1‑70 (MQ=255)
             gTGACGTTGCTGACCACTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACt                                                        <  1:2701036/71‑1 (MQ=255)
                   ttGCTGACCACTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTc                                                  <  1:993217/71‑1 (MQ=255)
                         aCCACTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCca                                             >  1:3105867/1‑70 (MQ=255)
                         aCCACTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCca                                             >  1:50078/1‑70 (MQ=255)
                             cTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAg                                        >  1:4428/1‑71 (MQ=255)
                             cTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAg                                        <  1:2638933/71‑1 (MQ=255)
                             cTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAg                                        <  1:173518/71‑1 (MQ=255)
                               aCGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTa                                                                    <  1:178006/41‑1 (MQ=255)
                                cGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCAt                                     <  1:561702/71‑1 (MQ=255)
                                       cTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCAcag                              <  1:1784506/71‑1 (MQ=255)
                                         gCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGa                                                            <  1:2925446/39‑1 (MQ=255)
                                         gCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGa                                                            <  1:2773659/39‑1 (MQ=255)
                                         gCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGa                                                            <  1:1436102/39‑1 (MQ=255)
                                         gCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGa                                                            <  1:2051809/39‑1 (MQ=255)
                                         gCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGa                                                            <  1:1621299/39‑1 (MQ=255)
                                          cAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCAcagca                            <  1:665313/70‑1 (MQ=255)
                                          cAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCAcagc                             <  1:1983193/69‑1 (MQ=255)
                                            gAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGct                         <  1:2406335/71‑1 (MQ=255)
                                             aaaCTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGct                         <  1:3135303/70‑1 (MQ=255)
                                              aaCTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGctct                       >  1:2182082/1‑71 (MQ=255)
                                                   ggCGTAGACGTATACGCAGTATCTACCGATACTCTCTTCACCCACAAAGCATGGCACAGCAGCTCTGAAAc                  <  1:1951645/71‑1 (MQ=255)
                                                   ggCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAAc                  <  1:445694/71‑1 (MQ=255)
                                                      gTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAg                                        <  1:344635/46‑1 (MQ=255)
                                                      gTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCAt               <  1:20965/71‑1 (MQ=255)
                                                       tAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACagcagc                          <  1:58211/59‑1 (MQ=255)
                                                       tAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACagcagc                          <  1:3043000/59‑1 (MQ=255)
                                                       tAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACagcagc                          <  1:1305028/59‑1 (MQ=255)
                                                       tAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACagcagc                          <  1:363100/59‑1 (MQ=255)
                                                         gACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGctc                        >  1:2820148/1‑59 (MQ=255)
                                                         gACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGc            >  1:1528994/1‑71 (MQ=255)
                                                          aCGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGc            >  1:2597824/1‑70 (MQ=255)
                                                          aCGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGc            >  1:662188/1‑70 (MQ=255)
                                                          aCGTATACGCAGTATCTACAGATACTCACTTCACCCACAAAGCATg                                    >  1:1434447/1‑46 (MQ=255)
                                                           cGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTa          <  1:2594907/71‑1 (MQ=255)
                                                           cGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTa          <  1:41012/71‑1 (MQ=255)
                                                            gTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAAc                  >  1:1398291/1‑62 (MQ=255)
                                                               taCGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGcaca                               <  1:681711/46‑1 (MQ=255)
                                                               taCGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGcaca                               <  1:409598/46‑1 (MQ=255)
                                                               taCGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGcaca                               <  1:2645020/46‑1 (MQ=255)
                                                               taCGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGcaca                               <  1:1484420/46‑1 (MQ=255)
                                                               taCGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCAcag                              <  1:1670999/47‑1 (MQ=255)
                                                               taCGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCAcag                              <  1:2104181/47‑1 (MQ=255)
                                                               taCGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCAcag                              <  1:3304841/47‑1 (MQ=255)
                                                               taCGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGa                     <  1:2865879/56‑1 (MQ=255)
                                                               taCGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGa                     <  1:1546206/56‑1 (MQ=255)
                                                               taCGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGa                     <  1:2503667/56‑1 (MQ=255)
                                                               taCGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTaaaa       <  1:3064258/70‑1 (MQ=255)
                                                               taCGCAGCATCTACCGATACTCACTTCACCCACAAAGCATgg                                   <  1:2051293/42‑1 (MQ=255)
                                                               taCGCAGCATCTACCGATACTCACTTCACCCACAAAGCATgg                                   <  1:728417/42‑1 (MQ=255)
                                                               taCGCAGCATCTACCGATACTCACTTCACCCACAAAGCATgg                                   <  1:2107622/42‑1 (MQ=255)
                                                                 cGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCAcagc                             >  1:1490578/1‑46 (MQ=255)
                                                                 cGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCAcagc                             >  1:2941860/1‑46 (MQ=255)
                                                                 cGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCAcagc                             >  1:3334328/1‑46 (MQ=255)
                                                                 cGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGaaa                   >  1:1986312/1‑56 (MQ=255)
                                                                 cGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGaaa                   >  1:278409/1‑56 (MQ=255)
                                                                 cGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGaaa                   >  1:684079/1‑56 (MQ=255)
                                                                 cGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTAAAATc     >  1:245566/1‑70 (MQ=255)
                                                                 cGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTAAAATCa    <  1:2875346/71‑1 (MQ=255)
                                                                 cGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTAAAATCa    <  1:6554/71‑1 (MQ=255)
                                                                 cGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTAAAATCa    <  1:2584100/71‑1 (MQ=255)
                                                                 cGCAGTATCTACCGATACTCACTTCACACACAAAGCAt                                     >  1:349283/1‑38 (MQ=255)
                                                                 cGCAGCATCTACCGATACTCACTTCACCCACAAAGCATGGca                                 >  1:3423910/1‑42 (MQ=255)
                                                                 cGCAGCATCTACCGATACTCACTTCACCCACAAAGCATGGca                                 >  1:1936089/1‑42 (MQ=255)
                                                                  gCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTAAAATCaa   <  1:1848852/71‑1 (MQ=255)
                                                                   cAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTAAAATCaaa  <  1:3313717/71‑1 (MQ=255)
                                                                    aGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTaa         >  1:591640/1‑63 (MQ=255)
                                                                    aGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTAAAATCaaa  >  1:1661141/1‑70 (MQ=255)
                                                                    aGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTAAAATCaaa  >  1:1027277/1‑70 (MQ=255)
                                                                      tATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGctct                       <  1:1208900/47‑1 (MQ=255)
                                                                      |                                                                   
CCGACCGAACTGGGTGACGTTGCTGACCACTACGAAGAACTGCAGAAACTGGGCGTAGACGTATACGCAGTATCTACCGATACTCACTTCACCCACAAAGCATGGCACAGCAGCTCTGAAACCATCGCTAAAATCAAA  >  minE/393983‑394120

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: