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‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACT  >  minE/1474797‑1474934
                                                                           |                                                                      
gTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGc                                                                     >  1:139728/1‑71 (MQ=255)
gTGCCGGATGCAGCGTGAACGCATTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGc                                                                     >  1:1755111/1‑71 (MQ=255)
     ggATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCg                                                                 >  1:2861943/1‑70 (MQ=255)
      gATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                               >  1:3642708/1‑71 (MQ=255)
      gATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                               <  1:2643893/71‑1 (MQ=255)
       aTGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCAt                                                              <  1:1398706/71‑1 (MQ=255)
           aGCGTAAACGCCTGATCCGTCCAAC‑‑‑‑GGTTCGG‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAg                                                           >  1:3457726/1‑71 (MQ=2)
            gCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCa                                                               <  1:2546502/71‑1 (MQ=0)
             cGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCAt                                                              >  1:3206446/1‑71 (MQ=0)
               taaaCGCCTGATCCGTCCAAC‑‑‑‑GGTTCGG‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGcat                                                       <  1:3519389/69‑2 (MQ=2)
                 aaCGCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAgg                                                    <  1:1702689/71‑1 (MQ=255)
                 aaCGCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAg                                                     <  1:2833867/70‑1 (MQ=255)
                 aaCGCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAg                                                     <  1:3225702/70‑1 (MQ=255)
                 aaCGCCTGATCCGTCCAAC‑‑‑‑GGTTCGG‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGcatct                                                     <  1:1790819/71‑4 (MQ=2)
                 aaCGCCTGATCCGTCCAAC‑‑‑‑GGTTCGG‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGcatc                                                      <  1:3415668/70‑3 (MQ=2)
                  aCGCCTGATCCGTCCAAC‑‑‑‑GGTTCGG‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCAt                                                              >  1:2284845/1‑61 (MQ=2)
                    gCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATg                                                 <  1:3217889/71‑1 (MQ=255)
                    gCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATg                                                 <  1:2401244/71‑1 (MQ=255)
                        tATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcggc                                             >  1:3370307/1‑71 (MQ=255)
                        tATCCGGCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcggc                                             >  1:1032332/1‑71 (MQ=255)
                              gCCTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATg                                                 <  1:2240725/61‑1 (MQ=255)
                               ccTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAAt                                      >  1:237384/1‑71 (MQ=255)
                                cTAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATg                                     <  1:2909558/71‑1 (MQ=255)
                                 tAC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTaa                                         <  1:1571474/66‑1 (MQ=255)
                                  aC‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAg                                                     <  1:3692041/53‑1 (MQ=255)
                                          cTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTAt                               <  1:2901458/71‑1 (MQ=255)
                                           tCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATc                              >  1:788138/1‑71 (MQ=255)
                                             gg‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATcc                             <  1:669431/70‑1 (MQ=255)
                                             gg‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCg                            <  1:3873256/71‑1 (MQ=255)
                                               g‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCgg                           <  1:294594/71‑1 (MQ=255)
                                               g‑‑ATTTGTAGGCCTGATAAGACGAGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCgg                           <  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‑‑‑‑GGCTCGG‑‑ATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACT  >  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: