Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,494,028 A→G 55.2% R710R (CGT→CGC nuoG ← NADH:ubiquinone oxidoreductase, chain G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,494,0280AG55.2% 8.4 / 36.5 49R710R (CGT→CGCnuoGNADH:ubiquinone oxidoreductase, chain G
Reads supporting (aligned to +/- strand):  ref base A (14/8);  new base G (14/13);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

AATATCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATCTTTGAT  >  minE/1493962‑1494096
                                                                  |                                                                    
aaTATCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCGCTGTAACg                                                                    >  1:1496060/1‑69 (MQ=255)
  tatCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATTGCGGTACGACCACTATAGCgat                                                                  <  1:1893086/69‑3 (MQ=255)
  tatCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCgat                                                                  <  1:1069036/69‑3 (MQ=255)
  tatCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCgat                                                                  <  1:927101/69‑3 (MQ=255)
  tatCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCgat                                                                  <  1:3431536/69‑3 (MQ=255)
  tatCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCgat                                                                  <  1:3301429/69‑3 (MQ=255)
  tatCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCgat                                                                  <  1:321256/69‑3 (MQ=255)
  tatCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCgat                                                                  <  1:2912230/69‑3 (MQ=255)
  tatCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCgat                                                                  <  1:290290/69‑3 (MQ=255)
  tatCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCgat                                                                  <  1:2587253/69‑3 (MQ=255)
  tatCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCgat                                                                  <  1:2411675/69‑3 (MQ=255)
  tatCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCgat                                                                  <  1:1766172/69‑3 (MQ=255)
   atCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCgat                                                                  <  1:2860703/68‑3 (MQ=255)
   atCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCgat                                                                  <  1:2569037/68‑3 (MQ=255)
           gCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCa                                                          <  1:310179/68‑1 (MQ=255)
           gCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCa                                                          <  1:753795/68‑1 (MQ=255)
           gCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCa                                                          <  1:592301/68‑1 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGGGGTTCACg                                                        >  1:102948/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCAc                                                         >  1:1290475/1‑68 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                        >  1:1321083/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                        >  1:2905904/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                        >  1:263277/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                        >  1:430518/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                        >  1:1695563/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                        >  1:2398458/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                        >  1:2301270/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                        >  1:700729/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                        >  1:2257405/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                        >  1:2157727/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                        >  1:2101837/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                        >  1:1996465/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                             >  1:596323/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                             >  1:379328/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                             >  1:2219973/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                             >  1:2018859/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                             >  1:340990/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                             >  1:2939931/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                             >  1:3353395/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                             >  1:1360882/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGa                              >  1:1516807/1‑68 (MQ=255)
                                                 gCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATcc                   >  1:3158047/1‑69 (MQ=255)
                                                    gTACGACCGCTGTAACGGTGCGGTTCACGGGCCTGTTTCTGCCCACGAATACGGAAt                            <  1:1393539/57‑1 (MQ=255)
                                                      aCGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCa                      >  1:2266224/1‑61 (MQ=255)
                                                      aCGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGgcgc              >  1:2306726/1‑69 (MQ=255)
                                                      aCGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGgcgc              >  1:74285/1‑69 (MQ=255)
                                                            gCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATc        <  1:2312072/69‑1 (MQ=255)
                                                            gCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATc        <  1:1996866/69‑1 (MQ=255)
                                                            gCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATc        <  1:783004/69‑1 (MQ=255)
                                                             cTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATc        <  1:724653/68‑1 (MQ=255)
                                                                  aCGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATCTTTGAt  <  1:2688970/69‑1 (MQ=255)
                                                                  aCGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATCTTTGAt  <  1:158854/69‑1 (MQ=255)
                                                                  aCGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATCTTTGAt  <  1:2428298/69‑1 (MQ=255)
                                                                  aCGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATCTTTGAt  <  1:171654/69‑1 (MQ=255)
                                                                  |                                                                    
AATATCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATCTTTGAT  >  minE/1493962‑1494096

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