Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 735,840 C→T 6.3% A275V (GCC→GTC)  holB → DNA polymerase III, delta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE735,8400CT 6.3% 246.3 / 4.6 79A275V (GCC→GTC) holBDNA polymerase III, delta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (41/33);  new base T (2/3);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

ATGGATGCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGCCACATTCGTG  >  minE/735773‑735908
                                                                   |                                                                    
atggatgCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCg                                                                    >  1:2983743/1‑70 (MQ=255)
atggatgCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCg                                                                    <  1:1582842/70‑1 (MQ=255)
atggatgCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGa                                                                   <  1:4656001/71‑1 (MQ=255)
atggatgCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGa                                                                   <  1:1736570/71‑1 (MQ=255)
   gatgCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAAc                                                                 <  1:1268392/70‑1 (MQ=255)
   gatgCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACt                                                                >  1:3619217/1‑71 (MQ=255)
        gcTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCaa                                                           <  1:2152927/71‑1 (MQ=255)
          tAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAAc                                                         >  1:4656103/1‑71 (MQ=255)
           aaaaCGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGa                                                                   >  1:159144/1‑60 (MQ=255)
           aaaaCGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAAc                                                         >  1:1052261/1‑70 (MQ=255)
           aaaaCGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAAc                                                         >  1:4457117/1‑70 (MQ=255)
            aaaCGCCATTATGGCGCTGCGCACGTGACCAATGTTGATGTGCCGGGGCTGGTCGTCGAACTGGCAAAcc                                                        <  1:3175774/70‑1 (MQ=255)
            aaaCGCCATTATGGCGCTGCGCACGTGACCAATGTTGATGTGCCGGGGCTGGTCGTCGAACTGGCAAAcc                                                        <  1:606893/70‑1 (MQ=255)
            aaaCGCCATTATGGCGCTGCGCACGTGACCAATGTTGACGTGCCGGGGCTGGTCGTCGAACTGGCAAAcc                                                        <  1:633901/70‑1 (MQ=255)
            aaaCGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAAcc                                                        <  1:2744894/70‑1 (MQ=255)
             aaCGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCAt                                                      <  1:1939982/71‑1 (MQ=255)
               cGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCt                                                    >  1:4343890/1‑71 (MQ=255)
               cGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCt                                                    >  1:3151713/1‑71 (MQ=255)
                gCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAAc                                                                 >  1:2899927/1‑57 (MQ=255)
                  catcatGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTtc                                                 <  1:296016/71‑1 (MQ=255)
                  catcatGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTtc                                                 <  1:3542979/71‑1 (MQ=255)
                  catcatGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTtc                                                 <  1:2809933/71‑1 (MQ=255)
                      atGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTccc                                             >  1:346623/1‑71 (MQ=255)
                      atGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTccc                                             >  1:1975284/1‑71 (MQ=255)
                      atGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTccc                                             >  1:4293849/1‑71 (MQ=255)
                      atGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTccc                                             >  1:3039797/1‑71 (MQ=255)
                         gTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGTACTGGCAAACCATCTTTCTCCCTcg                                          >  1:497023/1‑71 (MQ=255)
                          tgctgcGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTcgc                                         <  1:3449387/71‑1 (MQ=255)
                                    gTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAAc                                                                 >  1:3017146/1‑37 (MQ=255)
                                    gTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAAc                                                                 >  1:2143419/1‑37 (MQ=255)
                                    gTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGc                               <  1:3903185/71‑1 (MQ=255)
                                    gTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGc                               <  1:304103/71‑1 (MQ=255)
                                     tGACCAATGTTGATGTGCCGGGGCTGGTCGTCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAgg                                >  1:2084325/1‑69 (MQ=255)
                                     tGACCAATGTTGATGTGCCGGGGCTGGTCGTCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAgg                                >  1:647897/1‑69 (MQ=255)
                                       aCCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCaaa                                                          >  1:3495557/1‑41 (MQ=255)
                                       aCCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCta                             >  1:2158308/1‑70 (MQ=255)
                                        ccAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCtata                           <  1:4153765/71‑1 (MQ=255)
                                          aaTGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACt                         <  1:334595/71‑1 (MQ=255)
                                          aaTGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACt                         <  1:970830/71‑1 (MQ=255)
                                           aTGTTGATGTGCCGTGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTcgcg                                        <  1:1024932/55‑1 (MQ=255)
                                           aTGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCAt                                                      >  1:1015217/1‑41 (MQ=255)
                                           aTGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCAt                                                      >  1:2028352/1‑41 (MQ=255)
                                           aTGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTg                        <  1:2389435/71‑1 (MQ=255)
                                           aTGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTg                        >  1:594930/1‑71 (MQ=255)
                                            tGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCaaa                                                          <  1:263656/36‑1 (MQ=255)
                                            tGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCtat                            <  1:2016382/66‑1 (MQ=255)
                                             gTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGcc                                      >  1:641820/1‑55 (MQ=255)
                                              ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCttt                                                  >  1:1464179/1‑42 (MQ=255)
                                              ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCta                                           >  1:2203385/1‑48 (MQ=255)
                                              ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTc                                           >  1:4223320/1‑49 (MQ=255)
                                              ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCgcgc                                       >  1:4030272/1‑53 (MQ=255)
                                              ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTgggg                     >  1:1675258/1‑71 (MQ=255)
                                              ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTgggg                     <  1:1889239/71‑1 (MQ=255)
                                              ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTgggg                     >  1:4198250/1‑71 (MQ=255)
                                              ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTgggg                     >  1:2036966/1‑71 (MQ=255)
                                              ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTgggg                     >  1:36322/1‑71 (MQ=255)
                                              ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTgggg                     >  1:1069479/1‑71 (MQ=255)
                                                gATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATAc                          <  1:826098/64‑1 (MQ=255)
                                                gATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGa                   >  1:2578303/1‑71 (MQ=255)
                                                 aTGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGAt                  >  1:3043978/1‑71 (MQ=255)
                                                  tgtgCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATg                 <  1:3382121/71‑1 (MQ=255)
                                                  tgtgCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATg                 <  1:743936/71‑1 (MQ=255)
                                                   gtgCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCtat                            <  1:4332164/59‑1 (MQ=255)
                                                   gtgCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTAGCGCCTGCAGGCTATACTGGGGGATGt                >  1:1850522/1‑71 (MQ=255)
                                                    tgCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGtt               <  1:1802101/71‑1 (MQ=255)
                                                    tgCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGtt               <  1:779975/71‑1 (MQ=255)
                                                      ccGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTg             <  1:1269066/71‑1 (MQ=255)
                                                       cGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGc            <  1:303606/71‑1 (MQ=255)
                                                         ggCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGt                >  1:549724/1‑65 (MQ=255)
                                                         ggCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGt                >  1:555092/1‑65 (MQ=255)
                                                         ggCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGCca          <  1:489606/71‑1 (MQ=255)
                                                               gTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGCCACATTc     <  1:3194654/70‑1 (MQ=255)
                                                               gTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGCCACATTc     <  1:3020322/70‑1 (MQ=255)
                                                               gTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGCCACATTCg    >  1:2843136/1‑71 (MQ=255)
                                                               gTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGCCACATTCg    >  1:4633219/1‑71 (MQ=255)
                                                                tCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCt                                     >  1:4054341/1‑37 (MQ=255)
                                                                tCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGCCACATTCGt   >  1:1867224/1‑71 (MQ=255)
                                                                  gCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGCCACATTCGTg  <  1:3707524/70‑1 (MQ=255)
                                                                  gCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGCCACATTCGTg  >  1:611599/1‑70 (MQ=255)
                                                                   ccGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATAc                          <  1:2781748/45‑1 (MQ=255)
                                                                   |                                                                    
ATGGATGCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGCCACATTCGTG  >  minE/735773‑735908

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: