Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 931,730 G→A 31.9% R162R (CGC→CGT yddG ← predicted methyl viologen efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE931,7300GA31.9% 65.7 / 35.5 75R162R (CGC→CGTyddGpredicted methyl viologen efflux pump
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base A (15/9);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAACGCCAGGAAAT  >  minE/931662‑931797
                                                                    |                                                                   
gTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCcg                                                                     <  1:2711859/69‑1 (MQ=255)
          aGTAAACCCACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGtatt                                                           <  1:2853104/69‑1 (MQ=255)
             aaaCCCACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGtattta                                                         <  1:1834078/68‑1 (MQ=255)
                 cctcaGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGt                                                    <  1:2100134/66‑1 (MQ=255)
                 ccACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGt                                                    <  1:397379/69‑1 (MQ=255)
                 ccACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGt                                                    <  1:2205873/69‑1 (MQ=255)
                 ccACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGt                                                    <  1:1815197/69‑1 (MQ=255)
                 ccACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGt                                                    <  1:2810786/69‑1 (MQ=255)
                 ccACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGt                                                    <  1:2806860/69‑1 (MQ=255)
                 ccACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGt                                                    <  1:2619338/69‑1 (MQ=255)
                  cacaGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGt                                                    <  1:534170/68‑1 (MQ=255)
                   acaGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTa                                                  <  1:845488/69‑1 (MQ=255)
                   acaGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTa                                                  <  1:2091070/69‑1 (MQ=255)
                   acaGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTa                                                  <  1:2032990/69‑1 (MQ=255)
                   acaGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTa                                                  <  1:198512/69‑1 (MQ=255)
                   acaGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTa                                                  <  1:2559399/69‑1 (MQ=255)
                   acaGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTa                                                  <  1:2636689/69‑1 (MQ=255)
                   acaGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTa                                                  <  1:3076420/69‑1 (MQ=255)
                   acaGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTa                                                  <  1:1326098/69‑1 (MQ=255)
                   acaGACTAGCTCCCGTTAGCAGGAGAAAAACGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTa                                                  <  1:2158599/69‑1 (MQ=255)
                         ttGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGc                                            <  1:376657/69‑1 (MQ=255)
                         ttGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGc                                            <  1:2117419/69‑1 (MQ=255)
                                 gTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCt                                    <  1:2401520/69‑1 (MQ=255)
                                 gTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCt                                    <  1:833516/69‑1 (MQ=255)
                                           cAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAg                              >  1:3422534/1‑65 (MQ=255)
                                           cAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATaa                          >  1:2458729/1‑69 (MQ=255)
                                           cAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATaa                          >  1:3136359/1‑69 (MQ=255)
                                           cAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATaa                          >  1:3287022/1‑69 (MQ=255)
                                           cAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATaa                          >  1:106817/1‑69 (MQ=255)
                                           cAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATaa                          >  1:1644469/1‑69 (MQ=255)
                                           cAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATaa                          >  1:1762503/1‑69 (MQ=255)
                                           cAAAAACGGTAATTCCATTAAATCAGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATaa                          >  1:1667716/1‑69 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTATGCTGCCCAGATa                           >  1:2793897/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:2845365/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:470702/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:3424525/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:614558/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:3146589/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:8306/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:1480235/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:1647497/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:1821838/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:1895249/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:2199774/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:1779443/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:2409733/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                           >  1:2560653/1‑67 (MQ=255)
                                              aaaCGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACg                       <  1:2654154/69‑1 (MQ=255)
                                              aaaCGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACg                       <  1:2502361/69‑1 (MQ=255)
                                               aaCTGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACg                       <  1:1778174/68‑1 (MQ=255)
                                                  ggTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCAcc                   <  1:325339/69‑1 (MQ=255)
                                                  ggTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCAcc                   <  1:458434/69‑1 (MQ=255)
                                                  ggTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCAcc                   <  1:1207131/69‑1 (MQ=255)
                                                  ggTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCAcc                   <  1:1209938/69‑1 (MQ=255)
                                                  ggTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCAcc                   <  1:2171030/69‑1 (MQ=255)
                                                  ggTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCAcc                   <  1:1293095/69‑1 (MQ=255)
                                                  ggTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCAcc                   <  1:33881/69‑1 (MQ=255)
                                                  ggTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCAcc                   <  1:1580137/69‑1 (MQ=255)
                                                  ggTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCAcc                   <  1:687422/69‑1 (MQ=255)
                                                  ggTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCAcc                   <  1:1485809/69‑1 (MQ=255)
                                                        tCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGaa             >  1:3377981/1‑69 (MQ=255)
                                                        tCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGaa             >  1:1421991/1‑69 (MQ=255)
                                                        tCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGaa             >  1:2689989/1‑69 (MQ=255)
                                                        tCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGa              >  1:524349/1‑68 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc            >  1:661163/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc            >  1:1120555/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc            >  1:1070237/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc            >  1:10550/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc            >  1:3120090/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc            >  1:1392092/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc            >  1:3235810/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc            >  1:3193263/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc            >  1:2102905/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc            >  1:3109886/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc            >  1:1709393/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc            >  1:227830/1‑69 (MQ=255)
                                                                   cgcgTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAACGCCAGGAAAt  <  1:3411708/69‑1 (MQ=255)
                                                                   cgcgTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAACGCCAGGAAAt  <  1:982956/69‑1 (MQ=255)
                                                                    |                                                                   
GTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAACGCCAGGAAAT  >  minE/931662‑931797

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: