Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,573 C→A 11.1% A100A (GCG→GCT) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5730CA11.1% 337.4 / 16.6 108A100A (GCG→GCT) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (57/39);  new base A (5/7);  total (62/46)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

Only 100 of 110 total aligned reads displayed.

CTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATG  >  minE/312504‑312642
                                                                     |                                                                     
ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg                                                                      >  1:1844797/1‑71 (MQ=255)
ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg                                                                      >  1:734831/1‑71 (MQ=255)
  gAACTTCTTGCGGCAGCAACTGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                                     >  1:846843/1‑70 (MQ=255)
  gAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCa                                                                    >  1:265643/1‑71 (MQ=255)
   aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                                     <  1:859700/69‑1 (MQ=255)
   aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                                     <  1:854367/69‑1 (MQ=255)
   aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                                                     <  1:849499/69‑1 (MQ=255)
    aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                                   >  1:3482184/1‑70 (MQ=255)
    aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                                   >  1:2828559/1‑70 (MQ=255)
    aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                                   >  1:3351000/1‑70 (MQ=255)
    aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATc                                                                  <  1:1409466/71‑1 (MQ=255)
     cttcttGCGGCAGCAACGGTATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                                                   >  1:2097907/1‑69 (MQ=255)
       tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                                >  1:1851380/1‑70 (MQ=255)
        cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                                                                <  1:1706671/69‑1 (MQ=255)
        cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                               <  1:1343322/70‑1 (MQ=255)
        cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                               <  1:2299812/70‑1 (MQ=255)
        cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGt                                                               >  1:2067582/1‑70 (MQ=255)
        cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                                              >  1:10732/1‑71 (MQ=255)
         ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTTCACTTGCGCGATATCAGGATCGGt                                                               <  1:94435/68‑1 (MQ=255)
         ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtcc                                                             >  1:1607135/2‑69 (MQ=255)
         ctGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGtcc                                                             >  1:2054891/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)
              gcagcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGAc                                                        <  1:352053/71‑1 (MQ=16)
              gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                                        <  1:503449/71‑1 (MQ=255)
              gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                                        <  1:2401338/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:2187233/1‑70 (MQ=255)
                agcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACt                                                       >  1:3771861/1‑70 (MQ=255)
                 gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTc                                                      >  1:1205334/1‑70 (MQ=255)
                 gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                     <  1:3253030/71‑1 (MQ=255)
                 gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                                                     <  1:1549774/71‑1 (MQ=255)
                    aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                  >  1:684315/1‑71 (MQ=255)
                    aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                  >  1:962620/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)
                      ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt                                                >  1:651667/1‑71 (MQ=255)
                       gCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTc                                               <  1:811430/71‑1 (MQ=255)
                        cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCa                                              <  1:3712403/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)
                                 cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                                                    <  1:3310785/56‑1 (MQ=255)
                                   gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                                                  >  1:1511077/1‑56 (MQ=255)
                                   gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg                                    >  1:633201/1‑70 (MQ=255)
                                   gCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAAACTCAAAGGTCAGGGTGATCTg                                    >  1:2604417/1‑70 (MQ=255)
                                    gaGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACTGGAAgg                                                 >  1:201974/2‑56 (MQ=12)
                                    cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCa                                  >  1:1096764/1‑71 (MQ=255)
                                     aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                                                              >  1:2509207/1‑42 (MQ=255)
                                     aGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCAc                                 >  1:1640108/1‑71 (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)
                                             cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGc                                   >  1:1241334/1‑61 (MQ=255)
                                             cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc                          >  1:3303864/1‑70 (MQ=255)
                                             cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTAccc                          >  1:927970/1‑70 (MQ=255)
                                             cTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg                         >  1:541457/1‑71 (MQ=255)
                                              tGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGt                        >  1:54975/1‑71 (MQ=255)
                                              tGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGt                        >  1:1424638/1‑71 (MQ=255)
                                               gCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCg                         <  1:2500555/69‑1 (MQ=255)
                                                 aCCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                                                        <  1:3729467/36‑1 (MQ=255)
                                                 aCCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg                       >  1:1255845/1‑69 (MQ=255)
                                                 aCCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTg                       >  1:63853/1‑69 (MQ=255)
                                                   cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTa                             >  1:1609936/1‑61 (MQ=255)
                                                   cTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                   >  1:3256021/1‑71 (MQ=255)
                                                    tGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGt                   >  1:1926290/1‑70 (MQ=255)
                                                     gAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTc                  >  1:2730804/1‑70 (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:3169355/71‑1 (MQ=255)
                                                         cTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt             <  1:2635633/71‑1 (MQ=255)
                                                           gcgcGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa                 >  1:2965798/1‑65 (MQ=255)
                                                            cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGa                     <  1:2305820/60‑1 (MQ=255)
                                                            cgcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGt             >  1:3093525/1‑68 (MQ=255)
                                                             gcgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCa                 <  1:1074987/63‑1 (MQ=255)
                                                              cgATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTagag        <  1:2077407/71‑1 (MQ=255)
                                                                atatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa      >  1:2873614/1‑71 (MQ=255)
                                                                atatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGa      >  1:10652/1‑71 (MQ=255)
                                                                 tatCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt                                       <  1:2708622/37‑1 (MQ=255)
                                                                  atCCGCATCAGTACGAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACa    >  1:33773/1‑71 (MQ=255)
                                                                  atCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCAc                                 >  1:2586712/1‑42 (MQ=255)
                                                                     cGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATg  >  1:2575464/1‑70 (MQ=255)
                                                                     |                                                                     
CTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTGCACGGTACCCGTGGAGTCACTGTTAGAGGACATG  >  minE/312504‑312642

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: