Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,487,774 T→G 53.1% V418V (GTA→GTC nuoM ← NADH:ubiquinone oxidoreductase, membrane subunit M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,487,7740TG53.1% 9.2 / 49.3 49V418V (GTA→GTCnuoMNADH:ubiquinone oxidoreductase, membrane subunit M
Reads supporting (aligned to +/- strand):  ref base T (13/10);  new base G (20/6);  total (33/16)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

TAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGA  >  minE/1487708‑1487817
                                                                  |                                           
tAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTAc                                           <  1:2369600/69‑1 (MQ=255)
tAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTAc                                           <  1:2238524/69‑1 (MQ=255)
 aaCATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTAc                                           <  1:2607628/68‑1 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACTAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCt                                      >  1:1759892/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAAcc                                       >  1:652305/1‑66 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCt                                      >  1:459179/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCt                                      >  1:1642332/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCt                                      >  1:3320848/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCt                                      >  1:3020440/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCt                                      >  1:3010183/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCt                                      >  1:1939102/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCt                                      >  1:12824/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCt                                      >  1:1191441/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCt                                      >  1:33550/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCt                                      >  1:936146/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCt                                      >  1:3027587/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCt                                      >  1:2650528/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCt                                      >  1:899007/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCt                                      >  1:2426305/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCt                                      >  1:224742/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCt                                      >  1:173391/1‑67 (MQ=255)
       gCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACAGGGACAACCt                                      >  1:2698951/1‑67 (MQ=255)
                       aTGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGa                          >  1:3079824/1‑63 (MQ=255)
                       aTGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGa                          >  1:121941/1‑63 (MQ=255)
                       aTGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGa                          >  1:818340/1‑63 (MQ=255)
                       aTGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGa                          >  1:616769/1‑63 (MQ=255)
                       aTGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGa                          >  1:513879/1‑63 (MQ=255)
                       aTGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGa                          >  1:3360629/1‑63 (MQ=255)
                       aTGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGa                          >  1:1737997/1‑63 (MQ=255)
                       aTGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGa                          >  1:1881203/1‑63 (MQ=255)
                       aTGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGa                          >  1:2957852/1‑63 (MQ=255)
                       aTGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGa                          >  1:2969332/1‑63 (MQ=255)
                       aTGCAAAGACAAGCCCAAAGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGa                          >  1:2577091/1‑63 (MQ=255)
                                ccAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAAtt           >  1:2313499/1‑69 (MQ=255)
                                ccAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAAtt           >  1:1844567/1‑69 (MQ=255)
                                        aaGGTTGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCcgacg   <  1:838333/69‑1 (MQ=255)
                                        aaGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCcgacg   <  1:3307732/69‑1 (MQ=255)
                                        aaGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCcgacg   <  1:2767688/69‑1 (MQ=255)
                                        aaGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCcgacg   <  1:1786582/69‑1 (MQ=255)
                                        aaGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCcgacg   <  1:2427368/69‑1 (MQ=255)
                                        aaGGTAGAGATAACGGTAATCACCGGGACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCcgacg   <  1:1358274/69‑1 (MQ=255)
                                         aGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCcgacga  <  1:1382102/69‑1 (MQ=255)
                                         aGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCcgacga  <  1:711648/69‑1 (MQ=255)
                                          ggTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCcgacga  <  1:2467996/68‑1 (MQ=255)
                                                             aCCGGTACAACCTGGAAGCTGCCGAACAGAATCATaa              >  1:2327103/1‑37 (MQ=255)
                                                             aCCGGTACAACCTGGAAGCTGCCGAACAGAATCATaa              <  1:3155886/37‑1 (MQ=255)
                                                             aCCGGTACAACCTGGAAGCTGCCGAACAGAATCATaa              <  1:2652323/37‑1 (MQ=255)
                                                             aCCGGTACAACCTGGAAGCTGCCGAACAGAATCATaa              <  1:1887029/37‑1 (MQ=255)
                                                             aCCGGTACAACCTGGAAGCTGCCGAACAGAATCATaa              <  1:1565880/37‑1 (MQ=255)
                                                                  |                                           
TAACATCGCCAGCGAATAAACAGATGCAAAGACCAGCCCAAAGGTAGAGATAACGGTAATCACCGGTACAACCTGGAAGCTGCCGAACAGAATCATAAATTCGCCGACGA  >  minE/1487708‑1487817

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: