Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,474,864 C→T 15.5% intergenic (+162/‑72) 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,8640CT15.5% 128.4 / 10.2 45intergenic (+162/‑72)nrdA/nrdBribonucleoside diphosphate reductase 1, alpha subunit/ribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like
Reads supporting (aligned to +/- strand):  ref base C (16/22);  new base T (3/4);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACT  >  minE/1474797‑1474934
                                                                           |                                                                      
gTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGc                                                                     >  1:139728/1‑71 (MQ=255)
gTGCCGGATGCAGCGTGAACGCATTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGc                                                                     >  1:1755111/1‑71 (MQ=255)
     ggATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCg                                                                 >  1:2861943/1‑70 (MQ=255)
      gATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                               >  1:3642708/1‑71 (MQ=255)
      gATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                               <  1:2643893/71‑1 (MQ=255)
       aTGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCAt                                                              <  1:1398706/71‑1 (MQ=255)
           aGCGTAAACGCCTGATCCGTCCAAC‑‑‑‑GGTTCG‑G‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAg                                                           >  1:3457726/1‑71 (MQ=2)
            gCGTGAACGCCTTATCCGGCCTACGGATGGCGCG‑GGAA‑TTTGTAGGCCTGATAAGACGCGTTAGCGTCGCa                                                               <  1:2546502/71‑1 (MQ=0)
             cGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGA‑A‑TTTGTAGGCCTGATAAGACGCGTTAGCGTCGCAt                                                              >  1:3206446/1‑71 (MQ=0)
               taaaCGCCTGATCCGTCCAAC‑‑‑‑GGTTCG‑G‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGcat                                                       <  1:3519389/69‑2 (MQ=2)
                 aaCGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAgg                                                    <  1:1702689/71‑1 (MQ=255)
                 aaCGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAg                                                     <  1:2833867/70‑1 (MQ=255)
                 aaCGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAg                                                     <  1:3225702/70‑1 (MQ=255)
                 aaCGCCTGATCCGTCCAAC‑‑‑‑GGTTCG‑G‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGcatct                                                     <  1:1790819/71‑4 (MQ=2)
                 aaCGCCTGATCCGTCCAAC‑‑‑‑GGTTCG‑G‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGcatc                                                      <  1:3415668/70‑3 (MQ=2)
                  aCGCCTGATCCGTCCAAC‑‑‑‑GGTTCG‑G‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCAt                                                              >  1:2284845/1‑61 (MQ=2)
                    gCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATg                                                 <  1:3217889/71‑1 (MQ=255)
                    gCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATg                                                 <  1:2401244/71‑1 (MQ=255)
                        tATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcggc                                             >  1:3370307/1‑71 (MQ=255)
                        tATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcggc                                             >  1:1032332/1‑71 (MQ=255)
                              gCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATg                                                 <  1:2240725/61‑1 (MQ=255)
                               ccTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAAt                                      >  1:237384/1‑71 (MQ=255)
                                cTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATg                                     <  1:2909558/71‑1 (MQ=255)
                                 tAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTaa                                         <  1:1571474/66‑1 (MQ=255)
                                  aC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAg                                                     <  1:3692041/53‑1 (MQ=255)
                                          cTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTAt                               <  1:2901458/71‑1 (MQ=255)
                                           tCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATc                              >  1:788138/1‑71 (MQ=255)
                                             g‑g‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATcc                             <  1:669431/70‑1 (MQ=255)
                                             g‑g‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCg                            <  1:3873256/71‑1 (MQ=255)
                                               g‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCgg                           <  1:294594/71‑1 (MQ=255)
                                               g‑‑A‑TTTGTAGGCCTGATAAGACGAGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCgg                           <  1:604547/71‑1 (MQ=255)
                                                     ttGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCAt                        >  1:475180/1‑71 (MQ=255)
                                                        tAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTa                      <  1:901173/70‑1 (MQ=255)
                                                          ggCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAAc                   >  1:816354/1‑71 (MQ=255)
                                                           gCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACt                  <  1:910614/71‑1 (MQ=255)
                                                            ccTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTc                 >  1:1699658/1‑71 (MQ=255)
                                                                 tAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCaa             >  1:1346502/1‑70 (MQ=255)
                                                                   aGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAg          <  1:1082570/71‑1 (MQ=255)
                                                                     aCGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCa                         <  1:545070/54‑1 (MQ=255)
                                                                      cgcgCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGc                                    >  1:2492379/1‑42 (MQ=255)
                                                                      cgcgCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTa                                >  1:409811/1‑46 (MQ=255)
                                                                      cgcgCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTa                                >  1:2663659/1‑46 (MQ=255)
                                                                      cgcgCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACt                  <  1:2138498/60‑1 (MQ=255)
                                                                        cgcCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCaa             <  1:3782427/63‑1 (MQ=255)
                                                                         gcCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAAt                                      >  1:1750236/1‑37 (MQ=255)
                                                                           cAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACt  >  1:3345321/1‑71 (MQ=255)
                                                                           |                                                                      
GTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACT  >  minE/1474797‑1474934

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: ‑