Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 745,079 G→C 36.8% intergenic (+139/‑261) ycfP → / → ndh conserved hypothetical protein/respiratory NADH dehydrogenase 2/cupric reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE745,0790GC36.8% 43.3 / 47.6 57intergenic (+139/‑261)ycfP/ndhconserved hypothetical protein/respiratory NADH dehydrogenase 2/cupric reductase
Reads supporting (aligned to +/- strand):  ref base G (4/32);  new base C (14/7);  total (18/39)
Fisher's exact test for biased strand distribution p-value = 2.34e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAAC  >  minE/745013‑745130
                                                                  |                                                   
cACGCCGCATCCGACAGTCGAGCATGAATGCCTGATGTGCTTCTTATCAGGCCTACCGAACGCCCTGCa                                                   <  1:3170313/69‑1 (MQ=255)
cACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGTGCTTCTTATCAGGCCTACCGAACGCCCTGCa                                                   <  1:3010651/69‑1 (MQ=255)
cACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGTGCTTCTTATCAGGCCTACCGAACGCCCTGCa                                                   <  1:122114/69‑1 (MQ=255)
cACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGTGCTTCTTATCAGGCCTACCGAACGCCCTGCa                                                   <  1:1267154/69‑1 (MQ=255)
  cgccgcATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATa                                                 <  1:632879/69‑1 (MQ=255)
  cgccgcATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATa                                                 <  1:440980/69‑1 (MQ=255)
  cgccgcATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATa                                                 <  1:1458066/69‑1 (MQ=255)
  cgccgcATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATa                                                 <  1:2171373/69‑1 (MQ=255)
  cgccgcATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATa                                                 <  1:1205216/69‑1 (MQ=255)
  cgccgcATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATa                                                 <  1:767147/69‑1 (MQ=255)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                                             >  1:2400509/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                                             >  1:2384500/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                                             >  1:2307458/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                                             >  1:2830414/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                                             >  1:3181768/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                                             >  1:725454/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                                             >  1:37517/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                                             >  1:2698100/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                                             >  1:1132169/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                                             >  1:1441345/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                                             >  1:1251182/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                                             >  1:1347711/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacac                                              >  1:1025343/1‑68 (MQ=38)
      gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacac                                              >  1:1327001/1‑68 (MQ=38)
              cAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCActc                                     <  1:549887/69‑1 (MQ=255)
              cAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCActc                                     <  1:2906240/69‑1 (MQ=255)
                 tCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCActc                                     <  1:1204271/66‑1 (MQ=255)
                 tCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCtat                                  <  1:2869593/69‑1 (MQ=255)
                  cGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCtat                                  <  1:1001204/68‑1 (MQ=255)
                     gCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCa                                                   <  1:2971018/48‑1 (MQ=255)
                      cATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCAc                             >  1:535170/1‑69 (MQ=255)
                      cATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCAc                             >  1:3404765/1‑69 (MQ=255)
                              ccTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACaa                     >  1:1375572/1‑69 (MQ=255)
                              ccTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACaa                     >  1:1814275/1‑69 (MQ=255)
                                   tGCTCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct                          <  1:3309866/59‑1 (MQ=255)
                                   tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct                          <  1:1004564/59‑1 (MQ=255)
                                   tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct                          <  1:2468715/59‑1 (MQ=255)
                                   tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct                          <  1:795516/59‑1 (MQ=255)
                                   tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct                          <  1:226736/59‑1 (MQ=255)
                                   tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct                          <  1:18435/59‑1 (MQ=255)
                                   tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct                          <  1:1767813/59‑1 (MQ=255)
                                   tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct                          <  1:14381/59‑1 (MQ=255)
                                           ttATCATGACTACCGAACGCCCTGCATACACCCCTCActc                                     <  1:641559/40‑1 (MQ=255)
                                                aGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATaataaa   <  1:2516337/69‑1 (MQ=255)
                                                aGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATaataaa   <  1:886163/69‑1 (MQ=255)
                                                aGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATaataaa   <  1:2055170/69‑1 (MQ=255)
                                                aGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATaataaa   <  1:2049674/69‑1 (MQ=255)
                                                aGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATaataaa   <  1:1936844/69‑1 (MQ=255)
                                                aGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATaataaa   <  1:140745/69‑1 (MQ=255)
                                                 ggCCTACCGAACGCCGTGCATACACCGCTCTCTCTATATCACTCTCACAAATTCGCTCAAATaataaa   <  1:1861257/68‑1 (MQ=255)
                                                  gCCTACCGAATGCCCTCCATACACACCTCACTCTATATCACTCTCACAAATCAGCTCAAATaataaac  <  1:932385/68‑1 (MQ=255)
                                                  gCCTACCGAATGCCCTCCATACACACCTCACTCTATATCACTCTCACAAATCAGCTCAAATaataaac  <  1:598571/68‑1 (MQ=255)
                                                  gCCTACCGAATGCCCTCCATACACACCTCACTCTATATCACTCTCACAAATCAGCTCAAATaataaac  <  1:2619117/68‑1 (MQ=255)
                                                  gCCTACCGAATGCCCTCCATACACACCTCACTCTATATCACTCTCACAAATCAGCTCAAATaataaac  <  1:2567040/68‑1 (MQ=255)
                                                  gCCTACCGAATGCCCTCCATACACACCTCACTCTATATCACTCTCACAAATCAGCTCAAATaataaac  <  1:2340632/68‑1 (MQ=255)
                                                  gCCTACCGAATGCCCTCCATACACACCTCACTCTATATCACTCTCACAAATCAGCTCAAATaataaac  <  1:2187524/68‑1 (MQ=255)
                                                  gCCTACCGAATGCCCTCCATACACACCTCACTCTATATCACTCTCACAAATCAGCTCAAATaataaac  <  1:1597807/68‑1 (MQ=255)
                                                                  |                                                   
CACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAAC  >  minE/745013‑745130

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: