Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 931,721 A→G 31.5% N165N (AAT→AAC yddG ← predicted methyl viologen efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE931,7210AG31.5% 63.9 / 37.2 76N165N (AAT→AACyddGpredicted methyl viologen efflux pump
Reads supporting (aligned to +/- strand):  ref base A (24/28);  new base G (15/9);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAC  >  minE/931662‑931787
                                                           |                                                                  
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:1895249/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:2845365/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                 >  1:1647497/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                 >  1:2199774/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                 >  1:1821838/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                 >  1:2560653/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                 >  1:1779443/1‑67 (MQ=255)
                                            aaaaaCGGTAATTCCGTTAAATCCACGTGCGTATTTATTCGTTACTGTGCAGTAGGCTGCCCAGATa                 >  1:2409733/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:2689989/1‑69 (MQ=255)
                                                        tCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGaa   >  1:1421991/1‑69 (MQ=255)
                                                        tCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGa    >  1:524349/1‑68 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc  >  1:1070237/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc  >  1:1120555/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc  >  1:661163/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc  >  1:2102905/1‑69 (MQ=255)
                                                         ccATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAc  >  1:10550/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:3120090/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)
                                                           |                                                                  
GTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGACAAAAACGGTAATTCCATTAAATCCGCGTGCGTATTTATTCGTTACTGTGCAATAGGCTGCCCAGATAAACGCACCAATGAAC  >  minE/931662‑931787

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: