Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,474,863 C→G 8.4% intergenic (+161/‑73) nrdA → / → nrdB ribonucleoside diphosphate reductase 1, alpha subunit/ribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,474,8630CG 8.4% 144.0 / 2.5 47intergenic (+161/‑73)nrdA/nrdBribonucleoside diphosphate reductase 1, alpha subunit/ribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like
Reads supporting (aligned to +/- strand):  ref base C (20/23);  new base G (2/2);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACA  >  minE/1474794‑1474930
                                                                      |                                                                   
cGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACgcgc                                                                     >  1:2051645/1‑70 (MQ=255)
         gATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                           >  1:77001/1‑71 (MQ=255)
         gATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                           >  1:5542865/1‑71 (MQ=255)
          aTGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCAt                                                          >  1:661098/1‑71 (MQ=255)
             cAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAg                                                       >  1:2206749/1‑71 (MQ=255)
              aGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCAt                                                          >  1:5051680/1‑68 (MQ=2)
                  taaaCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTc                                                              >  1:3925360/3‑60 (MQ=2)
                  taaaCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGcat                                                   <  1:5142320/69‑2 (MQ=2)
                   gAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAg                                                 >  1:5173600/1‑71 (MQ=255)
                       gCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATg                                             >  1:1205345/1‑71 (MQ=255)
                        ccTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgc                                            >  1:5169314/1‑71 (MQ=255)
                         cTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcg                                           <  1:2114059/71‑1 (MQ=255)
                           tATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCg                                                                <  1:2806915/48‑1 (MQ=255)
                            aTCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACgcgc                                                                     >  1:5502626/1‑42 (MQ=255)
                            aTCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCAt                                                          >  1:862763/1‑53 (MQ=255)
                            aTCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcggcg                                        >  1:302316/1‑71 (MQ=255)
                            aTCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcggcg                                        >  1:4687063/1‑71 (MQ=255)
                               cGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTaa                                     <  1:445863/71‑1 (MQ=255)
                                  ccTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAAt                                  <  1:2538107/71‑1 (MQ=255)
                                   cTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                           <  1:3402006/45‑1 (MQ=255)
                                   cTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATg                                 >  1:4818100/1‑71 (MQ=255)
                                    tACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                           <  1:3046209/44‑1 (MQ=255)
                                    tACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAgg                                                <  1:3971003/55‑1 (MQ=255)
                                    aaCGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGc                                                     <  1:544935/50‑1 (MQ=11)
                                          tCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATc                          >  1:3569158/1‑71 (MQ=255)
                                            ggATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATcc                         <  1:3769936/70‑1 (MQ=255)
                                                   gTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCAtt                   >  1:3595493/1‑70 (MQ=255)
                                                      ggCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTaaa                <  1:3489876/70‑1 (MQ=255)
                                                      ggCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTaaa                <  1:2608719/70‑1 (MQ=255)
                                                       gCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACt              >  1:5246982/1‑71 (MQ=255)
                                                          tGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTaa                                     <  1:4666042/45‑1 (MQ=255)
                                                          tGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTaa                                     <  1:1805369/45‑1 (MQ=255)
                                                          tGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTaa                                     <  1:1866718/45‑1 (MQ=255)
                                                          tGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGcc                               <  1:3551191/51‑1 (MQ=255)
                                                           gATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAATGCCTTATCCGGCATTAAACTCCCaa         >  1:5427689/1‑71 (MQ=255)
                                                            aTAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGt                                       >  1:5213091/1‑41 (MQ=255)
                                                            aTAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGcc                               >  1:4217496/1‑49 (MQ=255)
                                                            aTAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGcc                               >  1:2805848/1‑49 (MQ=255)
                                                            aTAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTc             <  1:1711169/67‑1 (MQ=255)
                                                            aTAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCaa         <  1:9259/71‑1 (MQ=255)
                                                            aTAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCaa         >  1:1071461/1‑71 (MQ=255)
                                                             tAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCaa         <  1:1281463/70‑1 (MQ=255)
                                                               aGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAAc        <  1:2226338/69‑1 (MQ=255)
                                                               aGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAAc        <  1:4181682/69‑1 (MQ=255)
                                                                gACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTaa                                     <  1:3835883/39‑1 (MQ=255)
                                                                gACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAgg     <  1:5164120/71‑1 (MQ=255)
                                                                 aCGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTa                            >  1:224610/1‑47 (MQ=255)
                                                                  cgcgCCAGCGTCGCATCAGGCACAGGATGCGGCGTaa                                     <  1:5470319/37‑1 (MQ=255)
                                                                  cgcgCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCaa         <  1:1491326/65‑1 (MQ=255)
                                                                   gcgcCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGaca  <  1:3456300/71‑1 (MQ=255)
                                                                     gcCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAg      >  1:3188449/1‑65 (MQ=255)
                                                                      ccAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTa                            <  1:345226/42‑1 (MQ=255)
                                                                      |                                                                   
CGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACA  >  minE/1474794‑1474930

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