Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,492,939 A→G 50.8% L164L (TTG→CTG)  nuoH ← NADH:ubiquinone oxidoreductase, membrane subunit H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,492,9390AG50.8% ‑2.7 / 53.7 61L164L (TTG→CTG) nuoHNADH:ubiquinone oxidoreductase, membrane subunit H
Reads supporting (aligned to +/- strand):  ref base A (26/4);  new base G (13/18);  total (39/22)
Fisher's exact test for biased strand distribution p-value = 4.18e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

ATGCGCCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAG  >  minE/1492874‑1493006
                                                                 |                                                                   
aTGCGCCTGGCTGTTGACGATGTCGGTCATGTTAATGAACCGGCCTCGCCACCACGCCCATCAGGGaaa                                                                <  1:127607/69‑1 (MQ=255)
  gcgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:1104069/69‑1 (MQ=255)
  gcgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:756753/69‑1 (MQ=255)
  gcgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:459127/69‑1 (MQ=255)
  gcgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:447202/69‑1 (MQ=255)
  gcgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:2760770/69‑1 (MQ=255)
  gcgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:2645774/69‑1 (MQ=255)
  gcgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:2217639/69‑1 (MQ=255)
  gcgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:1899355/69‑1 (MQ=255)
  gcgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:1768387/69‑1 (MQ=255)
  gcgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:1477840/69‑1 (MQ=255)
  gcgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:1358760/69‑1 (MQ=255)
  gcgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:1204414/69‑1 (MQ=255)
   cgcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:1546384/68‑1 (MQ=255)
   cgcCTGGCTGTTGACGATGTCAGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAg                                                               <  1:1364951/69‑1 (MQ=255)
    gcCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:1198836/67‑1 (MQ=255)
     ccTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:233800/66‑1 (MQ=255)
            gTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:2751701/59‑1 (MQ=255)
               gACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAacac                                                   >  1:3382885/1‑69 (MQ=255)
               gACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAacac                                                   >  1:2779067/1‑69 (MQ=255)
               gACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAacac                                                   >  1:3173269/1‑69 (MQ=255)
               gACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAacac                                                   >  1:2398742/1‑69 (MQ=255)
               gACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAacac                                                   >  1:2239756/1‑69 (MQ=255)
               gACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAacac                                                   >  1:1141370/1‑69 (MQ=255)
               gACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAacac                                                   >  1:559265/1‑69 (MQ=255)
               gACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAacac                                                   >  1:1068401/1‑69 (MQ=255)
               gACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAacac                                                   >  1:909854/1‑69 (MQ=255)
               gACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAacac                                                   >  1:1435499/1‑69 (MQ=255)
                 cGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGaaa                                                                <  1:1642008/54‑1 (MQ=255)
                      tCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAACACTTCGTAg                                            >  1:2275527/1‑69 (MQ=255)
                      tCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAACACTTCGTAg                                            >  1:1285207/1‑69 (MQ=255)
                      tCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAACACTTCGTAg                                            >  1:651434/1‑69 (MQ=255)
                           cATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAg                                       >  1:2363741/1‑69 (MQ=255)
                           cATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAg                                       >  1:3075943/1‑69 (MQ=255)
                            aTGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:96464/1‑65 (MQ=255)
                            aTGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:3318524/1‑65 (MQ=255)
                            aTGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:3310308/1‑65 (MQ=255)
                            aTGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:3277299/1‑65 (MQ=255)
                            aTGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:2572541/1‑65 (MQ=255)
                            aTGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:1959475/1‑65 (MQ=255)
                            aTGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:2561877/1‑65 (MQ=255)
                            aTGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:2284996/1‑65 (MQ=255)
                            aTGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:2533127/1‑65 (MQ=255)
                            aTGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:2416452/1‑65 (MQ=255)
                            aTGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:2503369/1‑65 (MQ=255)
                            aTGTTGAATGAACCGGCCTGCGACACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:2758540/1‑65 (MQ=255)
                            aTGTTGAATGAACCGGACTGCGCCACCACGCCCATCAGGGAAAGCCCGAGGAACACTTCGTAGCt                                          >  1:789388/1‑65 (MQ=255)
                                                       acGCCCATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCAccc           >  1:2054496/1‑69 (MQ=255)
                                                       acGCCCATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCAccc           >  1:546617/1‑69 (MQ=255)
                                                       acGCCCATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCAccc           >  1:698851/1‑69 (MQ=255)
                                                       acGCCCATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCAccc           >  1:3112439/1‑69 (MQ=255)
                                                       acGCCCATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCAccc           >  1:2975107/1‑69 (MQ=255)
                                                          cccATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTgcgc                              >  1:871093/1‑47 (MQ=255)
                                                          cccATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTgcgc                              >  1:466042/1‑47 (MQ=255)
                                                          cccATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTgcgc                              >  1:2231905/1‑47 (MQ=255)
                                                          cccATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCAca                     >  1:622895/1‑55 (MQ=255)
                                                          cccATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACg                     >  1:768869/1‑56 (MQ=255)
                                                          cccATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACg                     >  1:1444612/1‑56 (MQ=255)
                                                          cccATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGc        <  1:1204671/69‑1 (MQ=255)
                                                           ccATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGc        <  1:2582290/68‑1 (MQ=255)
                                                           ccATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGc        <  1:2432977/68‑1 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCAt                  >  1:931459/1‑53 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:701919/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  <  1:2537578/69‑1 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  <  1:149751/69‑1 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:2863698/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:820330/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:2867432/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:3096963/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  <  1:1618380/69‑1 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:755209/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:1715236/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:626097/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  <  1:3221334/69‑1 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:1244126/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:3243204/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:3315868/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:335548/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  >  1:451995/1‑69 (MQ=255)
                                                                aaGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACCCATCGCACCCAGCAACGAg  <  1:784936/69‑1 (MQ=255)
                                                                 aGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAg  <  1:1234071/68‑1 (MQ=255)
                                                                 |                                                                   
ATGCGCCTGGCTGTTGACGATGTCGGTCATGTTGAATGAACCGGCCTGCGCCACCACGCCCATCAAGGAAAGCCCGAGGAACACTTCGTAGCTCAGGGTCTGCGCAGAAGCACGCATCGCACCCAGCAACGAG  >  minE/1492874‑1493006

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: