Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,494,025 G→A 55.1% Y711Y (TAC→TAT nuoG ← NADH:ubiquinone oxidoreductase, chain G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,494,0250GA55.1% 4.7 / 37.9 49Y711Y (TAC→TATnuoGNADH:ubiquinone oxidoreductase, chain G
Reads supporting (aligned to +/- strand):  ref base G (14/8);  new base A (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 = 1.00e+00

AATATCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATC  >  minE/1493962‑1494090
                                                               |                                                                 
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:592301/68‑1 (MQ=255)
           gCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCa                                                    <  1:310179/68‑1 (MQ=255)
           gCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCa                                                    <  1:753795/68‑1 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGGGGTTCACg                                                  >  1:102948/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCAc                                                   >  1:1290475/1‑68 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                  >  1:2301270/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                  >  1:430518/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                  >  1:1321083/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                  >  1:1996465/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                  >  1:2905904/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                  >  1:2101837/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                  >  1:263277/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                  >  1:1695563/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                  >  1:2157727/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                  >  1:2398458/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                  >  1:700729/1‑69 (MQ=255)
            cTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCACTATAGCGATGCGGTTCACg                                                  >  1:2257405/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                       >  1:2219973/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                       >  1:2939931/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                       >  1:2018859/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                       >  1:1360882/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                       >  1:596323/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                       >  1:3353395/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                       >  1:340990/1‑69 (MQ=255)
                                       ggCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGaa                       >  1:379328/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:783004/69‑1 (MQ=255)
                                                            gCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATc  <  1:1996866/69‑1 (MQ=255)
                                                             cTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATc  <  1:724653/68‑1 (MQ=255)
                                                               |                                                                 
AATATCCTGCGGCTGACGCGGCTCATGAACGCTGATATTGGCGCGCATGGCGGTACGACCGCTGTAACGGTGCGGTTCACGGGCCAGTTTCTGCCCACGAATACGGAATGTCGCATCCGGCGCAGCATC  >  minE/1493962‑1494090

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: