Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 745,050 C→T 13.2% intergenic (+110/‑290) 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,0500CT13.2% 93.8 / 3.7 52intergenic (+110/‑290)ycfP/ndhconserved hypothetical protein/respiratory NADH dehydrogenase 2/cupric reductase
Reads supporting (aligned to +/- strand):  ref base C (18/27);  new base T (3/4);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

CCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAA  >  minE/744985‑745111
                                                                 |                                                             
ccGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGTGct                                                            >  1:878421/1‑69 (MQ=255)
ccGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGTGct                                                            >  1:1855860/1‑69 (MQ=255)
ccGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGTGct                                                            >  1:672090/1‑69 (MQ=255)
                     aGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACg                                       <  1:2384987/69‑1 (MQ=255)
                     aGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACg                                       <  1:1460074/69‑1 (MQ=255)
                       gCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGcc                                     <  1:1152572/69‑1 (MQ=255)
                       gCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGcc                                     <  1:1313805/69‑1 (MQ=255)
                       gCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGcc                                     <  1:1881942/69‑1 (MQ=255)
                       gCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGcc                                     <  1:1878002/69‑1 (MQ=255)
                       gCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGcc                                     <  1:2908361/69‑1 (MQ=255)
                            cACGCCGCATCCGACAGTCGAGCATGAATGCCTGATGTGCTTCTTATCAGGCCTACCGAACGCCCTGCa                                <  1:3170313/69‑1 (MQ=255)
                            cACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGTGCTTCTTATCAGGCCTACCGAACGCCCTGCa                                <  1:122114/69‑1 (MQ=255)
                            cACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGTGCTTCTTATCAGGCCTACCGAACGCCCTGCa                                <  1:1267154/69‑1 (MQ=255)
                            cACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGTGCTTCTTATCAGGCCTACCGAACGCCCTGCa                                <  1:3010651/69‑1 (MQ=255)
                              cgccgcATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATa                              <  1:1205216/69‑1 (MQ=255)
                              cgccgcATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATa                              <  1:2171373/69‑1 (MQ=255)
                              cgccgcATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATa                              <  1:767147/69‑1 (MQ=255)
                              cgccgcATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATa                              <  1:1458066/69‑1 (MQ=255)
                              cgccgcATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATa                              <  1:440980/69‑1 (MQ=255)
                              cgccgcATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATa                              <  1:632879/69‑1 (MQ=255)
                                  gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                          >  1:725454/1‑68 (MQ=38)
                                  gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                          >  1:2830414/1‑68 (MQ=38)
                                  gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                          >  1:2698100/1‑68 (MQ=38)
                                  gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                          >  1:3181768/1‑68 (MQ=38)
                                  gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                          >  1:37517/1‑68 (MQ=38)
                                  gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                          >  1:2400509/1‑68 (MQ=38)
                                  gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                          >  1:2307458/1‑68 (MQ=38)
                                  gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                          >  1:1132169/1‑68 (MQ=38)
                                  gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                          >  1:1251182/1‑68 (MQ=38)
                                  gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                          >  1:1441345/1‑68 (MQ=38)
                                  gcATCCGACAGTCGAGCATCAATGCCTGATGCGCCTCTTATCAGGCCTACCGAATGCCCTCCATacaca                          >  1:2384500/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:2906240/69‑1 (MQ=255)
                                          cAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCActc                  <  1:549887/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:14381/59‑1 (MQ=255)
                                                               tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct       <  1:1767813/59‑1 (MQ=255)
                                                               tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct       <  1:18435/59‑1 (MQ=255)
                                                               tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct       <  1:2468715/59‑1 (MQ=255)
                                                               tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct       <  1:226736/59‑1 (MQ=255)
                                                               tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct       <  1:795516/59‑1 (MQ=255)
                                                               tGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCActct       <  1:1004564/59‑1 (MQ=255)
                                                                 |                                                             
CCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAA  >  minE/744985‑745111

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: