Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 746,152 A→G 59.6% V271V (GTA→GTG ndh → respiratory NADH dehydrogenase 2/cupric reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE746,1520AG59.6% 13.2 / 21.9 37V271V (GTA→GTGndhrespiratory NADH dehydrogenase 2/cupric reductase
Reads supporting (aligned to +/- strand):  ref base A (6/9);  new base G (14/8);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCAGTGCTGATGAAGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGAAAGATATCGGTGGTCTTGAAACTAACCGTATCA  >  minE/746086‑746219
                                                                  |                                                                   
cACCAGTGCTGATGAAGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTATg                                                                   >  1:2329544/1‑69 (MQ=255)
cACCAGTGCTGATGAAGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTATg                                                                   >  1:614081/1‑69 (MQ=255)
              aaGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCa                                                     <  1:987766/69‑1 (MQ=255)
              aaGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCa                                                     <  1:2577765/69‑1 (MQ=255)
              aaGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCa                                                     <  1:2778160/69‑1 (MQ=255)
              aaGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCa                                                     <  1:1973977/69‑1 (MQ=255)
              aaGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCa                                                     <  1:1935824/69‑1 (MQ=255)
              aaGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCa                                                     <  1:3256437/69‑1 (MQ=255)
              aaGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCa                                                     <  1:352719/69‑1 (MQ=255)
              aaGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTATGGGCAGCCGGGAt                                                       >  1:1267218/1‑67 (MQ=255)
               aGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCa                                                     <  1:1949596/68‑1 (MQ=255)
                ggcggcCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTATGGGCAGCCGTGATCaa                                                    <  1:6401/68‑1 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCa                                                               >  1:448720/1‑42 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGCCTTCCt                                     >  1:2795012/1‑68 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGACTTcc                                      >  1:2985171/1‑67 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGACTTCCt                                     >  1:3002597/1‑68 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGACTTCCt                                     >  1:2977802/1‑68 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGACTTCCt                                     >  1:2937578/1‑68 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGACTTCCt                                     >  1:1110469/1‑68 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGACTTCCt                                     >  1:185524/1‑68 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGACTTCCt                                     >  1:1091745/1‑68 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGACTTCCt                                     >  1:1413857/1‑68 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGACTTCCt                                     >  1:1722197/1‑68 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGACTTCCt                                     >  1:254716/1‑68 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGACTTCCt                                     >  1:2089383/1‑68 (MQ=255)
                               aaaGATGGCGAATATATTGAGGCTGATCTGATGGTGTGGGCAGCCGGGATCAAAGCGCCTGACTTCCt                                     >  1:2148212/1‑68 (MQ=255)
                                aaGATGGCGAATATATTGAGGCTGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCt                                     >  1:1857969/1‑67 (MQ=255)
                                aaGATGGCGAATATATTGAGGCTGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCt                                     >  1:2992182/1‑67 (MQ=255)
                                aaGATGGCGAATATATTGAGGCTGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCt                                     >  1:2622822/1‑67 (MQ=255)
                                      gCGAATATATTGAGGCTGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGaaa                                 <  1:1288606/65‑1 (MQ=255)
                                           tatatTGAGGCTGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGAAAGATATCggt                        <  1:3178234/69‑1 (MQ=255)
                                           tatatTGAGGCTGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGAAAGATATCggt                        <  1:1742472/69‑1 (MQ=255)
                                           tatatTGAGGCTGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGAAAGATATCggt                        <  1:3298685/69‑1 (MQ=255)
                                           tatatTGAGGCTGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGAAAGATATCggt                        <  1:2251173/69‑1 (MQ=255)
                                           tatatTGAGGCTGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGAAAGATATCggt                        <  1:781198/69‑1 (MQ=255)
                                           tatatTGAGGATGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGAAAGATATCggt                        <  1:1664316/69‑1 (MQ=255)
                                            atatTGAGGCTGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGAAAGATATCggt                        <  1:674574/68‑1 (MQ=255)
                                                                  aTGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGAAAGATATCGGTGGTCTTGAAACTAACCGTATCa  <  1:1964760/68‑1 (MQ=255)
                                                                  aTGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGAAAGATATCGGTGGTCTTGAAACTAACCGTATCa  <  1:2782783/68‑1 (MQ=255)
                                                                  aTGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGAAAGATATCGGTGGTCTTGAAACTAACCGTATCa  <  1:2759206/68‑1 (MQ=255)
                                                                  |                                                                   
CACCAGTGCTGATGAAGGCGGCCTGCACACTAAAGATGGCGAATATATTGAGGCTGATCTGATGGTATGGGCAGCCGGGATCAAAGCGCCAGACTTCCTGAAAGATATCGGTGGTCTTGAAACTAACCGTATCA  >  minE/746086‑746219

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: