Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 735,832 C→G 6.2% G272G (GGC→GGG holB → DNA polymerase III, delta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE735,8320CG 6.2% 235.1 / 4.0 80G272G (GGC→GGGholBDNA polymerase III, delta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (43/32);  new base G (2/3);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

ACGTTGCTGATGGATGCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGCCA  >  minE/735764‑735900
                                                                    |                                                                    
aCGTTGCTGATGGATGCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCt                                                                    >  1:2475812/1‑71 (MQ=255)
aCGTTGCTGATGGATGCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCt                                                                    >  1:4458884/1‑71 (MQ=255)
    tGCTGATGGATGCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTc                                                                >  1:3213893/1‑71 (MQ=255)
     gCTGATGGATGCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCTGGCCTGGTCg                                                               >  1:149844/1‑71 (MQ=255)
         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)
                  cTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTg                                                                   <  1:2869841/54‑1 (MQ=255)
                   tAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAAc                                                 >  1:4656103/1‑71 (MQ=255)
                    aaaaCGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGt                                                                 >  1:3254970/1‑54 (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)
                          ccATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTc                                                                <  1:4026932/49‑1 (MQ=255)
                           catcatGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTtc                                         <  1:3542979/71‑1 (MQ=255)
                           catcatGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTtc                                         <  1:296016/71‑1 (MQ=255)
                           catcatGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTtc                                         <  1:2809933/71‑1 (MQ=255)
                             tcatGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGt                                                                 >  1:3903279/1‑45 (MQ=255)
                               atGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGt                                                                 >  1:1151239/1‑43 (MQ=255)
                               atGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTccc                                     >  1:3039797/1‑71 (MQ=255)
                               atGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTccc                                     >  1:346623/1‑71 (MQ=255)
                               atGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTccc                                     >  1:4293849/1‑71 (MQ=255)
                               atGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTccc                                     >  1:1975284/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:1889239/71‑1 (MQ=255)
                                                       ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTgggg             >  1:1069479/1‑71 (MQ=255)
                                                       ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTgggg             >  1:1675258/1‑71 (MQ=255)
                                                       ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTgggg             >  1:36322/1‑71 (MQ=255)
                                                       ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTgggg             >  1:4198250/1‑71 (MQ=255)
                                                       ttGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTgggg             >  1:2036966/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:743936/71‑1 (MQ=255)
                                                           tgtgCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATg         <  1:3382121/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)
                                                                    |                                                                    
ACGTTGCTGATGGATGCGCTAAAACGCCATCATGGTGCTGCGCAGGTGACCAATGTTGATGTGCCGGGCCTGGTCGCCGAACTGGCAAACCATCTTTCTCCCTCGCGCCTGCAGGCTATACTGGGGGATGTTTGCCA  >  minE/735764‑735900

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: