Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,513 T→C 9.8% Q120Q (CAA→CAG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5130TC 9.8% 252.7 / 13.9 92Q120Q (CAA→CAGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (50/33);  new base C (6/3);  total (56/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 106 total aligned reads displayed.

GATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAC  >  minE/312444‑312583
                                                                     |                                                                      
gATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtctt                                                                       >  1:5520150/1‑71 (MQ=255)
gATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtctt                                                                       >  1:2008611/1‑71 (MQ=255)
gATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtct                                                                        >  1:3463587/1‑70 (MQ=255)
 aTAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtctt                                                                       <  1:3206255/70‑1 (MQ=255)
  tAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGc                                                                     <  1:724190/71‑1 (MQ=255)
    aCGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCgg                                                                   <  1:3627125/71‑1 (MQ=255)
     cGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCgg                                                                   >  1:816479/1‑70 (MQ=255)
      gCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCgg                                                                   >  1:83276/1‑69 (MQ=255)
      gCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAAAGCTCACCCCTTGCTGCTGAACTTCTTGCgg                                                                   >  1:465464/1‑69 (MQ=255)
            aaCCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGc                                                                     <  1:4396384/61‑1 (MQ=255)
              ccATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCgg                                                                   >  1:3797976/1‑61 (MQ=255)
               cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa                                                        >  1:4944481/1‑71 (MQ=255)
               cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa                                                        >  1:4614853/1‑71 (MQ=255)
                aTCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCAt                                                       <  1:4390367/71‑1 (MQ=255)
                aTCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCAt                                                       >  1:1760925/1‑71 (MQ=255)
                   aGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtctt                                                                       >  1:4292374/1‑52 (MQ=255)
                   aGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCg                                                     <  1:1136844/70‑1 (MQ=255)
                    ggAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGcc                                                   <  1:1255017/71‑1 (MQ=255)
                     gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCa                                                  <  1:4812912/71‑1 (MQ=255)
                        gctgctGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCaa                                                             <  1:395066/57‑1 (MQ=255)
                           gctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACg                                                           >  1:1248655/1‑56 (MQ=255)
                            ctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATc                                                      >  1:899414/1‑60 (MQ=255)
                            ctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATc                                                      >  1:4102938/1‑60 (MQ=255)
                             tgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGc                                                                     >  1:2005831/1‑44 (MQ=255)
                             tgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGc                                                                     >  1:1138464/1‑44 (MQ=255)
                             tgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGc                                                                     >  1:4150828/1‑44 (MQ=255)
                             tgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGt                                          >  1:3063409/1‑71 (MQ=255)
                             tgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGt                                          >  1:563107/1‑71 (MQ=255)
                              ggATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAg                                           <  1:1328137/69‑1 (MQ=255)
                              ggATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGtt                                         >  1:2160592/1‑71 (MQ=255)
                               gatgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTg                                                                      <  1:2282503/41‑1 (MQ=255)
                               gatgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTg                                                                      <  1:4504649/41‑1 (MQ=255)
                               gatgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTg                                                                      <  1:2077533/41‑1 (MQ=255)
                               gatgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTg                                                                      <  1:2900383/41‑1 (MQ=255)
                               gatgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGttt                                        >  1:1003350/1‑71 (MQ=255)
                                atgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                       <  1:168728/71‑1 (MQ=255)
                                 tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt                                      <  1:2535880/71‑1 (MQ=255)
                                  gatTTCTCAACGCTCACCCCTTGCTGCTGAActtct                                                                        >  1:155326/1‑36 (MQ=255)
                                  gatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGtt                                     >  1:2136240/1‑71 (MQ=255)
                                   atTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGgcagc                                                               <  1:949662/44‑1 (MQ=255)
                                   atTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCt                                               >  1:4356838/1‑60 (MQ=255)
                                   atTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCt                                               >  1:4721322/1‑60 (MQ=255)
                                   atTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTc                                    >  1:3648995/1‑71 (MQ=255)
                                    tttCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCt                                   >  1:2938909/1‑71 (MQ=255)
                                      tctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTg                                  >  1:1653928/1‑70 (MQ=255)
                                      tctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTg                                  >  1:4675381/1‑70 (MQ=255)
                                       ttcAACACTGATCCCCTGCTGCTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCa                                >  1:2213787/2‑71 (MQ=255)
                                       ttcAACACTGATCCCCTGCTGCTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCa                                >  1:3473372/2‑71 (MQ=255)
                                       ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc                                 <  1:2942190/70‑1 (MQ=255)
                                          aaCGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCt                             >  1:5522174/1‑71 (MQ=255)
                                           aCGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCAt                                                       >  1:3613632/1‑44 (MQ=255)
                                            cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTg                            <  1:381890/70‑1 (MQ=255)
                                            cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                           <  1:5650579/71‑1 (MQ=255)
                                            cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                           <  1:4243671/71‑1 (MQ=255)
                                            cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                           <  1:1460206/71‑1 (MQ=255)
                                             gCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGaa                          >  1:2302376/1‑71 (MQ=255)
                                              cTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGt                                          <  1:5464836/54‑1 (MQ=255)
                                              cTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                         >  1:5181787/1‑71 (MQ=255)
                                              cTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                         >  1:5188030/1‑71 (MQ=255)
                                              cTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                         >  1:4415100/1‑71 (MQ=255)
                                               tCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAcc                        >  1:2304822/1‑71 (MQ=255)
                                                 aCCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCt                       <  1:1199429/70‑1 (MQ=255)
                                                 aCCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                      >  1:124613/1‑71 (MQ=255)
                                                 aCCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                      >  1:444578/1‑71 (MQ=255)
                                                   cccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                       <  1:195482/52‑1 (MQ=255)
                                                   cccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                       <  1:1997438/52‑1 (MQ=255)
                                                   cccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcg                    >  1:5373523/1‑71 (MQ=255)
                                                   cccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcg                    >  1:3268437/1‑71 (MQ=255)
                                                    ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                   >  1:1201616/1‑71 (MQ=255)
                                                    ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                   <  1:1253011/71‑1 (MQ=255)
                                                    ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                   >  1:3131442/1‑71 (MQ=255)
                                                    ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcg                    >  1:630402/1‑70 (MQ=255)
                                                     cTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCAcc                              >  1:2674422/1‑59 (MQ=255)
                                                      ttGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGtt                                         <  1:3360678/47‑1 (MQ=255)
                                                      ttGCTGCTCAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                           <  1:2550526/61‑1 (MQ=255)
                                                       tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTg                                              <  1:238437/41‑1 (MQ=19)
                                                       tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTg                                              <  1:914114/41‑1 (MQ=19)
                                                       tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGa                 >  1:2052532/1‑70 (MQ=255)
                                                        gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGt                                          >  1:3854876/1‑44 (MQ=255)
                                                        gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCt                             >  1:442104/1‑57 (MQ=255)
                                                        gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGata               >  1:3822746/1‑71 (MQ=255)
                                                        gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGata               >  1:5065724/1‑71 (MQ=255)
                                                        gctgctGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGtt                                         <  1:2431612/45‑1 (MQ=16)
                                                        gctgctGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTc                                    >  1:3609236/1‑50 (MQ=19)
                                                        gctgctGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTc                                    >  1:3360919/1‑50 (MQ=19)
                                                         ctgctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGatat              >  1:3947693/1‑71 (MQ=255)
                                                           gctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGt                                          <  1:3670099/41‑1 (MQ=255)
                                                           gctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGt                                          <  1:4092937/41‑1 (MQ=255)
                                                           gctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGt                                          <  1:4430184/41‑1 (MQ=255)
                                                           gctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                         <  1:4098675/58‑1 (MQ=255)
                                                            ctgAACTTCTTGCGGCAGCACCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGCACCTGCGCGATATCCg           >  1:2945349/1‑71 (MQ=255)
                                                                aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt        >  1:448469/1‑70 (MQ=255)
                                                                aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt        >  1:4642006/1‑70 (MQ=255)
                                                                 cttcttGCGGCTGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc                                 <  1:277137/44‑1 (MQ=255)
                                                                 cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc                                 <  1:3342115/44‑1 (MQ=255)
                                                                  ttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                           <  1:2553563/49‑1 (MQ=255)
                                                                   tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                         >  1:1774050/1‑50 (MQ=255)
                                                                   tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                         >  1:1260472/1‑50 (MQ=255)
                                                                    cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg     <  1:4296053/69‑1 (MQ=255)
                                                                     ttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc  >  1:5200020/1‑71 (MQ=255)
                                                                     |                                                                      
GATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAC  >  minE/312444‑312583

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