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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5610CT16.4% 288.9 / 40.5 116Q104Q (CAG→CAAacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (58/39);  new base T (9/10);  total (67/49)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 118 total aligned reads displayed.

ACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGT  >  minE/312493‑312631
                                                                    |                                                                      
atccccTGCTGCTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTg                                                                      >  1:3585923/3‑71 (MQ=255)
aCCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                                                                      >  1:240322/1‑71 (MQ=255)
  cccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                                                     <  1:9975/70‑1 (MQ=255)
   ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                                                                   <  1:3724069/71‑1 (MQ=255)
      tgctgcTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGa                                                                 >  1:2814008/1‑70 (MQ=255)
      tgctgcTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGa                                                                 >  1:1081387/1‑70 (MQ=255)
      tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGa                                                                 >  1:114107/1‑70 (MQ=255)
        ctgctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGatat                                                              <  1:1885283/71‑1 (MQ=255)
        ctgctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGata                                                               <  1:3637246/70‑1 (MQ=255)
        ctgctgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCAGCACCTGCACTTGCGCGatat                                                              <  1:2365335/71‑1 (MQ=255)
          gctgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTgcgc                                                                   <  1:2523939/64‑1 (MQ=255)
           ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                                                                      >  1:3640926/1‑60 (MQ=255)
           ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg                                                           >  1:1844797/1‑71 (MQ=255)
           ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg                                                           >  1:734831/1‑71 (MQ=255)
            tgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                                                >  1:1130171/1‑65 (MQ=255)
            tgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATcag                                                           >  1:3561527/1‑68 (MQ=255)
             gAACTTCTTGCGGCAGCAACTGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                          >  1:846843/1‑70 (MQ=255)
             gAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCa                                                         >  1:265643/1‑71 (MQ=255)
              aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                          <  1:849499/69‑1 (MQ=255)
              aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                          <  1:854367/69‑1 (MQ=255)
              aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                          <  1:859700/69‑1 (MQ=255)
              aaCCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATcagg                                                          <  1:598713/69‑4 (MQ=19)
               aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                        >  1:3351000/1‑70 (MQ=255)
               aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                        >  1:2828559/1‑70 (MQ=255)
               aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATc                                                       <  1:1409466/71‑1 (MQ=255)
                cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCt                                                                       >  1:2098160/1‑54 (MQ=255)
                  tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                     >  1:1851380/1‑70 (MQ=255)
                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                     <  1:1706671/69‑1 (MQ=255)
                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                    >  1:2067582/1‑70 (MQ=255)
                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                    <  1:1148549/70‑1 (MQ=255)
                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                    <  1:2299812/70‑1 (MQ=255)
                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                    <  1:1343322/70‑1 (MQ=255)
                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                                   >  1:10732/1‑71 (MQ=255)
                    ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTTCACTTGCGCGATATCAGGATCGGt                                                    <  1:94435/68‑1 (MQ=255)
                    ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtcc                                                  >  1:2054891/2‑69 (MQ=255)
                    ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtcc                                                  >  1:1607135/2‑69 (MQ=255)
                    ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtc                                                   <  1:1321715/69‑2 (MQ=16)
                      gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                  >  1:1500614/1‑69 (MQ=255)
                      gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                                                  >  1:2540908/1‑69 (MQ=255)
                       cGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccgg                                               <  1:2303622/71‑3 (MQ=255)
                       cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                                                <  1:1935227/70‑1 (MQ=255)
                       cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAg                                               <  1:3529200/71‑1 (MQ=255)
                        ggCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                                              >  1:484397/1‑71 (MQ=16)
                         gcagcaTCGGCATCGCCAGCTGCAGTTTGTTGTGCACCTGAACCTGCGCGATATCCGCa                                                         <  1:1993315/59‑1 (MQ=255)
                         gcagcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGAc                                             <  1:352053/71‑1 (MQ=16)
                         gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                             <  1:2401338/71‑1 (MQ=255)
                         gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                             <  1:503449/71‑1 (MQ=255)
                          cagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                             >  1:315622/1‑70 (MQ=255)
                          cagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACt                                            <  1:2520448/71‑1 (MQ=255)
                           agcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATc                                                       >  1:3313306/1‑59 (MQ=255)
                           agcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACt                                            >  1:3771861/1‑70 (MQ=255)
                           agcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACt                                            >  1:2187233/1‑70 (MQ=255)
                            gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTc                                           >  1:1205334/1‑70 (MQ=255)
                            gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                          <  1:1549774/71‑1 (MQ=255)
                            gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                          <  1:3253030/71‑1 (MQ=255)
                               aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                       >  1:962620/1‑71 (MQ=255)
                               aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                       >  1:689156/1‑71 (MQ=255)
                               aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                       >  1:684315/1‑71 (MQ=255)
                                cGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAgg                                      >  1:605973/1‑71 (MQ=255)
                                cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAgg                                      >  1:1005238/1‑71 (MQ=255)
                                cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCACGATATCCGCATCAGTACCAGACTCAAAgg                                      >  1:460897/1‑71 (MQ=255)
                                cGGCATCGCCAGCTGCAGTTTGTTATGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAgg                                      >  1:3402512/1‑71 (MQ=255)
                                 ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                     >  1:651667/1‑71 (MQ=255)
                                  gCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTc                                    <  1:811430/71‑1 (MQ=255)
                                      cGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACt                                            >  1:598168/1‑59 (MQ=255)
                                       gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaaa                                        >  1:3284913/1‑62 (MQ=255)
                                        cgaGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAGGTAAGGGta                              <  1:2853999/69‑2 (MQ=12)
                                        ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg                              <  1:1102180/71‑1 (MQ=255)
                                         gaGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAGGTAAGGGtaa                             <  1:639487/70‑3 (MQ=255)
                                         gaGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAGGTAAGGGtaa                             <  1:2916913/70‑3 (MQ=255)
                                              gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                       >  1:1511077/1‑56 (MQ=255)
                                              gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAAACTCAAAGGTCAGGGTGATCTg                         >  1:2604417/1‑70 (MQ=255)
                                               gaGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAgg                                      >  1:201974/2‑56 (MQ=12)
                                                aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                                   >  1:2509207/1‑42 (MQ=255)
                                                aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCAc                      >  1:1255720/1‑71 (MQ=255)
                                                  tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                     <  1:2747186/54‑1 (MQ=255)
                                                  tttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                     <  1:949551/54‑1 (MQ=255)
                                                   ttGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                     >  1:2051342/1‑37 (MQ=255)
                                                    tGTTCTGCACCTGAACCTGCGCGATATCCGCATCa                                                      >  1:3421359/1‑35 (MQ=255)
                                                    tGTTCTGCACCTGAACCTGCGCGATATCCGCATCa                                                      >  1:1010098/1‑35 (MQ=255)
                                                    tGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                  >  1:1017383/1‑71 (MQ=255)
                                                     gTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                  <  1:188869/70‑1 (MQ=255)
                                                       tCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAGTCAAAGGTCAggg                                <  1:3450658/54‑1 (MQ=255)
                                                        cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc               >  1:927970/1‑70 (MQ=255)
                                                        cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc               >  1:3303864/1‑70 (MQ=255)
                                                        cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg              >  1:541457/1‑71 (MQ=255)
                                                         tGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGt             >  1:54975/1‑71 (MQ=255)
                                                          gCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg              <  1:2500555/69‑1 (MQ=255)
                                                          gCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg              <  1:2697790/69‑1 (MQ=255)
                                                            aCCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                             <  1:3729467/36‑1 (MQ=255)
                                                            aCCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg            >  1:63853/1‑69 (MQ=255)
                                                            aCCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg            >  1:1255845/1‑69 (MQ=255)
                                                              cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                          <  1:1536096/37‑1 (MQ=255)
                                                              cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                  >  1:1609936/1‑61 (MQ=255)
                                                              cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt        >  1:3256021/1‑71 (MQ=255)
                                                                gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa      >  1:1004204/1‑71 (MQ=255)
                                                                 aaCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTc       <  1:1145281/69‑1 (MQ=255)
                                                                 aaCCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCACGGTGATCTGCACGGTACCCGTGGAGTCa      <  1:2806951/70‑1 (MQ=255)
                                                                    cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt  <  1:2635633/71‑1 (MQ=255)
                                                                    cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt  <  1:3169355/71‑1 (MQ=255)
                                                                    |                                                                      
ACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGT  >  minE/312493‑312631

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: