Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 931,645 T→C 31.1% I191V (ATA→GTA)  yddG ← predicted methyl viologen efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE931,6450TC31.1% 65.2 / 37.4 77I191V (ATA→GTA) yddGpredicted methyl viologen efflux pump
Reads supporting (aligned to +/- strand):  ref base T (38/15);  new base C (6/18);  total (44/33)
Fisher's exact test for biased strand distribution p-value = 1.73e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCATAAGCAAATCCTAAGGTAAATGCCGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGAC  >  minE/931579‑931705
                                                                  |                                                            
cagcaTAAGCAAATCCTAAGGTAAATGCCGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCa                                                            >  1:2035728/1‑69 (MQ=255)
cagcaTAAGCAAATCCTAAGGTAAATGCCGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCa                                                            >  1:3026087/1‑69 (MQ=255)
cagaaTAAGCAAATCCTAAGGTAAATGCCGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCa                                                            >  1:960868/1‑69 (MQ=255)
                     tAAATGCCGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGaa                                       <  1:1951561/69‑1 (MQ=255)
                     tAAATGCCGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGaa                                       <  1:3419719/69‑1 (MQ=255)
                     tAAATGCCGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGaa                                       <  1:770760/69‑1 (MQ=255)
                     tAAATGCCGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGaa                                       <  1:706429/69‑1 (MQ=255)
                           ccGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAAt                                    >  1:48857/1‑66 (MQ=255)
                           ccGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATa                                   >  1:337048/1‑67 (MQ=255)
                           ccGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAg                                  >  1:1425446/1‑68 (MQ=255)
                           ccGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAg                                  >  1:692968/1‑68 (MQ=255)
                           ccGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAg                                  >  1:1005065/1‑68 (MQ=255)
                           ccGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAg                                  >  1:814129/1‑68 (MQ=255)
                           ccGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAg                                  >  1:3348907/1‑68 (MQ=255)
                           acGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAg                                  >  1:1517952/2‑68 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:2233779/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:2537602/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:2205859/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:3087883/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:3108513/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:3292892/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:2074635/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:2728369/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:3328804/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:188738/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:1868001/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:377343/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:1806323/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:1635869/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:645106/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:1214176/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATTTCTGGt                                                    <  1:1190779/41‑1 (MQ=255)
                                    tGAGTTTAATCATGACGGGCGTGCTAAATACCATATCTGGt                                                    <  1:436358/41‑1 (MQ=255)
                                     gAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAg                       >  1:3310420/1‑69 (MQ=255)
                                     gAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAg                       >  1:386506/1‑69 (MQ=255)
                                     gAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAg                       >  1:452916/1‑69 (MQ=255)
                                     gAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAg                       >  1:2723888/1‑69 (MQ=255)
                                     gAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAg                       >  1:1073570/1‑69 (MQ=255)
                                     gAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAg                       >  1:1384425/1‑69 (MQ=255)
                                     gAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAg                       >  1:1500221/1‑69 (MQ=255)
                                     gAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAg                       >  1:2183203/1‑69 (MQ=255)
                                      aGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGa                      <  1:1859440/69‑1 (MQ=255)
                                      aGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGa                      >  1:642609/1‑69 (MQ=255)
                                      aGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGa                      <  1:536331/69‑1 (MQ=255)
                                      aGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGa                      >  1:531451/1‑69 (MQ=255)
                                      aGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGa                      <  1:1666495/69‑1 (MQ=255)
                                      aGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGa                      >  1:2492265/1‑69 (MQ=255)
                                      aGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGa                      >  1:1876001/1‑69 (MQ=255)
                                      aGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGa                      >  1:201680/1‑69 (MQ=255)
                                      aGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGa                      >  1:2036988/1‑69 (MQ=255)
                                      aGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGa                      >  1:2150697/1‑69 (MQ=255)
                                      aGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGa                      <  1:27612/69‑1 (MQ=255)
                                        tttAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGaa                                       >  1:242466/1‑50 (MQ=255)
                                         ttAATCATGACGGGCGTGCTAAATACCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACtt                   >  1:903691/1‑69 (MQ=255)
                                         ttAATCATGACGGGCGTGCTAAATACCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACtt                   >  1:1240082/1‑69 (MQ=255)
                                         ttAATCATGACGGGCGTGCTAAATACCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACtt                   >  1:193424/1‑69 (MQ=255)
                                         ttAATCATGACGGGCGTGCTAAATACCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACtt                   >  1:393691/1‑69 (MQ=255)
                                         ttAATCATGACGGGCGTGCTAAATACCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACtt                   >  1:397672/1‑69 (MQ=255)
                                         ttAATCATGACGGGCGTGCTAAATACCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACtt                   >  1:409392/1‑69 (MQ=255)
                                                 gACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGt           <  1:1394433/69‑1 (MQ=255)
                                                 gACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGt           <  1:2246635/69‑1 (MQ=255)
                                                 gACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGt           <  1:3041162/69‑1 (MQ=255)
                                                    gggCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAg        <  1:2483470/69‑1 (MQ=255)
                                                    gggCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAg        <  1:852930/69‑1 (MQ=255)
                                                        gTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAgg    <  1:731814/69‑1 (MQ=255)
                                                        gTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACACACAGACTTGCTCCCGTTAGCAgg    <  1:3187167/69‑1 (MQ=255)
                                                           cTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTa                                >  1:2026973/1‑38 (MQ=255)
                                                           cTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGAc  >  1:1719340/1‑68 (MQ=255)
                                                           cTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGAc  >  1:1723745/1‑68 (MQ=255)
                                                           cTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGAc  >  1:2961209/1‑68 (MQ=255)
                                                           cTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGAc  >  1:1831256/1‑68 (MQ=255)
                                                           cTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGAc  >  1:3045380/1‑68 (MQ=255)
                                                           cTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGAc  >  1:3373461/1‑68 (MQ=255)
                                                           cTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGAc  >  1:2426667/1‑68 (MQ=255)
                                                           cTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGAc  >  1:773600/1‑68 (MQ=255)
                                                           cTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGAc  >  1:1119961/1‑68 (MQ=255)
                                                           cTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGAc  >  1:2090210/1‑68 (MQ=255)
                                                                  tCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAg                       <  1:3293418/40‑1 (MQ=255)
                                                                  |                                                            
CAGCATAAGCAAATCCTAAGGTAAATGCCGCAGAGATGAGTTTAATCATGACGGGCGTGCTAAATATCATTTCTGGTTGTGGCGTAAGAAAATAGTAAACCCACAGACTTGCTCCCGTTAGCAGGAC  >  minE/931579‑931705

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: