Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,292,502 G→A 7.3% T219M (ACG→ATG)  yeeO ← predicted multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,292,5020GA 7.3% 216.7 / 4.8 68T219M (ACG→ATG) yeeOpredicted multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (29/34);  new base A (3/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTCCCTGCACCACGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGCGACGAAATCAATAAT  >  minE/1292434‑1292572
                                                                    |                                                                      
tATTCCCTGCACCACGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTc                                                                      <  1:1979358/71‑1 (MQ=255)
     ccTGCACCACGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTc                                                                 >  1:2118557/1‑71 (MQ=255)
      cTGCACCACGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTc                                                                 >  1:891679/1‑70 (MQ=255)
        gCACCACGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGCTAACTGAGTACCGTCAGCTc                                                                 >  1:3705217/1‑68 (MQ=255)
        gCACCACGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTc                                                                 >  1:1091388/1‑68 (MQ=255)
          accacGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGAt                                                            <  1:2114775/71‑1 (MQ=255)
          accacGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGAt                                                            <  1:1258798/71‑1 (MQ=255)
           ccacGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGAt                                                            <  1:3102917/70‑1 (MQ=255)
            cacGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGAt                                                            >  1:1314528/1‑69 (MQ=255)
             acGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAg                                                         >  1:2163651/1‑71 (MQ=255)
              cGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCAAGATAAGt                                                        >  1:479809/1‑71 (MQ=255)
               gAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGt                                                        >  1:1034213/1‑70 (MQ=255)
                aaGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTc                                                       >  1:1084408/1‑70 (MQ=255)
                  gTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAAc                                                    <  1:1470061/71‑1 (MQ=255)
                       ccGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCg                                                                        <  1:309810/46‑1 (MQ=255)
                       ccGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAg                                               <  1:901180/71‑1 (MQ=255)
                        cGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGa                                                             <  1:983536/56‑1 (MQ=255)
                        cGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATa                                                           <  1:2456365/58‑1 (MQ=255)
                        cGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATa                                                           <  1:723331/58‑1 (MQ=255)
                         gCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTg                                             <  1:138232/71‑1 (MQ=255)
                             ccAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAATTCAACGCCAGTGCtt                                          >  1:1454218/1‑70 (MQ=255)
                             ccAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCtt                                          >  1:318667/1‑70 (MQ=255)
                             ccAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCtt                                          >  1:3426678/1‑70 (MQ=255)
                               aaTAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCaa                                                     >  1:3199500/1‑57 (MQ=255)
                                 tAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACt                                     <  1:273682/71‑1 (MQ=255)
                                   agagTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACtt                                    <  1:2057985/70‑1 (MQ=255)
                                   agagTGATGGCAGCTGCTGGATAACTGAGTACCGTAAGCTCCAGATAAGTCAAc                                                    >  1:3348300/1‑54 (MQ=255)
                                    gagTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGAt                                                            <  1:1906704/45‑1 (MQ=255)
                                    gagTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGAt                                                            <  1:2934683/45‑1 (MQ=255)
                                    gagTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCt                                  <  1:10183/71‑1 (MQ=255)
                                       tgttGGCAGCTGCTGGATAACTGAGTACCATCAGCTCCAGATAAGTCAACGCCAGTGCTTTAAc                                      <  1:1655978/61‑1 (MQ=255)
                                       tGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTc                               <  1:1001788/71‑1 (MQ=255)
                                       tGATGGCAGCTGCTGGATAACTGAGTACCATCAGCTCCAGATAAGTCAACGCCAGTGCTTTAAc                                      <  1:1663211/64‑1 (MQ=255)
                                         aTGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTc                                                       <  1:539189/45‑1 (MQ=255)
                                         aTGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAg                                               >  1:1695584/1‑53 (MQ=255)
                                         aTGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAg                                               >  1:1185273/1‑53 (MQ=255)
                                         aTGGCAGCTGCTGGATAACTGAGTACCATCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACtt                                    >  1:3353226/1‑64 (MQ=255)
                                         aTGGCAGCTGCTGGATAACTGAGTACCATCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACtt                                    >  1:3500911/1‑64 (MQ=255)
                                         aTGGCAGCTGCTGGATAACTGAGTACCATCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACtt                                    >  1:2866786/1‑64 (MQ=255)
                                          tGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCttt                                         <  1:144710/58‑1 (MQ=255)
                                          tGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCttt                                         <  1:1871339/58‑1 (MQ=255)
                                          tGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCttt                                         <  1:2198012/58‑1 (MQ=255)
                                          tGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCttt                                         <  1:6293/58‑1 (MQ=255)
                                               gctgctGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCa                         >  1:2552099/1‑69 (MQ=255)
                                               gctgctGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCAt                        >  1:481428/1‑70 (MQ=255)
                                               gctgctGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCAGTCGTGGCAt                        >  1:1551821/1‑70 (MQ=255)
                                                 tgctgGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCAc                     >  1:2230481/1‑71 (MQ=255)
                                                 tgctgGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCAAGTGCTTTAACTTCTGTCGTGGCATCAc                     >  1:3220444/1‑71 (MQ=255)
                                                  gctgGATAACTGAGTACCGTCAGCTCCAGATAAGTCaa                                                     >  1:3612307/1‑38 (MQ=255)
                                                  gctgGATAACTGAGTACCGTCAGCTCCAGATAAGTCaa                                                     >  1:3114295/1‑38 (MQ=255)
                                                   ctgGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCAcc                    >  1:694806/1‑70 (MQ=255)
                                                      gATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCcgc                 <  1:3603417/70‑1 (MQ=255)
                                                       aTAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTg                            >  1:203542/1‑58 (MQ=255)
                                                        tAACTGAGTACCGTCAGCTCAAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGcga               >  1:1630334/1‑70 (MQ=255)
                                                        tAACTGAGTACCGTCAGCACCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCg                              >  1:3689731/1‑55 (MQ=255)
                                                         aaCTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTa                                        <  1:3670803/44‑1 (MQ=255)
                                                         aaCTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTa                                        <  1:3241638/44‑1 (MQ=255)
                                                         aaCTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACtt                                    <  1:1276192/48‑1 (MQ=255)
                                                         aaCTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCcgc                 >  1:3520485/1‑67 (MQ=255)
                                                          aCTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGcgacga            <  1:281971/71‑1 (MQ=255)
                                                           cTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGCGACGaa           <  1:1403084/71‑1 (MQ=255)
                                                           cTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGCGACGaa           <  1:2295069/71‑1 (MQ=255)
                                                           cTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCATTAACt                                     >  1:1287034/1‑45 (MQ=255)
                                                              aGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGCGACGAAAt         <  1:1215818/70‑1 (MQ=255)
                                                              aGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGCGACGAAAt         <  1:2901828/70‑1 (MQ=255)
                                                              aGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGCGACGAAATc        <  1:13548/71‑1 (MQ=255)
                                                                   cGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGCGACGAAATCaataa   <  1:1547008/71‑1 (MQ=255)
                                                                   cGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGCGACGAAATCaataa   <  1:880061/71‑1 (MQ=255)
                                                                   cGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGACGTGGCATCACCCGCGACGAAATCaataa   <  1:662949/71‑1 (MQ=255)
                                                                    gTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGCGACGAAATCaataat  >  1:2803880/1‑71 (MQ=255)
                                                                    gTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGCGACGAAATCaataa   >  1:1655522/1‑70 (MQ=255)
                                                                    |                                                                      
TATTCCCTGCACCACGAAGTGCCCCGCTACCAATAAGAGTGATGGCAGCTGCTGGATAACTGAGTACCGTCAGCTCCAGATAAGTCAACGCCAGTGCTTTAACTTCTGTCGTGGCATCACCCGCGACGAAATCAATAAT  >  minE/1292434‑1292572

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: