Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,474,835 C→T 14.3% intergenic (+133/‑101) 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,8350CT14.3% 80.5 / 4.0 28intergenic (+133/‑101)nrdA/nrdBribonucleoside diphosphate reductase 1, alpha subunit/ribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like
Reads supporting (aligned to +/- strand):  ref base C (12/12);  new base T (2/2);  total (14/14)
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

TAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCT  >  minE/1474766‑1474902
                                                                     |                                                                    
tAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCt                                                                     <  1:3644908/71‑1 (MQ=255)
  aGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCg                                                                   <  1:3323169/71‑1 (MQ=255)
     cgcgCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGAt                                                               >  1:2683499/1‑71 (MQ=255)
       cgcCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGAttt                                                             <  1:1382229/71‑1 (MQ=255)
        gcCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTg                                                            <  1:2149011/71‑1 (MQ=255)
        gcCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTg                                                            <  1:2094629/71‑1 (MQ=255)
                      aGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAg                                               <  1:778086/70‑1 (MQ=255)
                                  ccGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTc                                  >  1:2642924/1‑71 (MQ=255)
                                     gATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                               >  1:3807445/1‑71 (MQ=255)
                                     gATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                               <  1:1285296/71‑1 (MQ=255)
                                     gATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                               <  1:1280417/71‑1 (MQ=255)
                                       cgtaGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCAt                              >  1:1274066/4‑71 (MQ=2)
                                             gtaaaCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGca                        <  1:3461498/68‑1 (MQ=2)
                                              taaaCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGAcg                                            <  1:4043089/48‑1 (MQ=2)
                                               gAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCAt                              <  1:1152260/62‑1 (MQ=255)
                                                aaCGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAg                     >  1:3381782/1‑70 (MQ=255)
                                                aaCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGcatct                     >  1:1589145/1‑68 (MQ=2)
                                                    ccTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgc                >  1:3008940/1‑71 (MQ=255)
                                                     cTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcg               >  1:930117/1‑71 (MQ=255)
                                                       tATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcggc             >  1:968122/1‑71 (MQ=255)
                                                       tATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcggc             >  1:2947588/1‑71 (MQ=255)
                                                       tATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcggc             >  1:1572154/1‑71 (MQ=255)
                                                           cGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTaa         <  1:3979128/71‑1 (MQ=255)
                                                           cGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTa          >  1:1836969/1‑70 (MQ=255)
                                                             gCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTaaaa       >  1:1262225/1‑71 (MQ=255)
                                                                  cGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCt  >  1:1734719/1‑71 (MQ=255)
                                                                    gCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTc                                  <  1:4120725/37‑1 (MQ=255)
                                                                    gCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTc                                  <  1:2788727/37‑1 (MQ=255)
                                                                     |                                                                    
TAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCT  >  minE/1474766‑1474902

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