Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 313,254 T→C 10.0% I279V (ATC→GTC)  acrA ← multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE313,2540TC10.0% 193.7 / 9.9 70I279V (ATC→GTC) acrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (37/26);  new base C (4/3);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTG  >  minE/313188‑313323
                                                                  |                                                                     
gCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAgcg                                                                   >  1:1054775/1‑71 (MQ=255)
  ccTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAg                                                                     >  1:2731923/1‑67 (MQ=255)
  ccTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGcgcg                                                                 >  1:2621553/1‑71 (MQ=255)
   cttcttCCAGACGTGCGCGCACGAACATACCCGGAAGCAGAGTGTGATCCGGGTTCGGGAAGATAGcgcg                                                                 >  1:2840704/1‑70 (MQ=255)
      cttcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAg                                                              <  1:3262776/70‑1 (MQ=255)
       ttcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAggg                                                            >  1:38203/1‑71 (MQ=255)
        tcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGt                                                           >  1:1407321/1‑71 (MQ=255)
         ccAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAggg                                                            <  1:1014023/69‑1 (MQ=255)
         ccAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGt                                                           <  1:1319217/70‑1 (MQ=255)
          cAGACGTGCGCGCACGAACATACCCGGCAGCAGTGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTg                                                          <  1:177122/70‑1 (MQ=255)
          cAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTg                                                          <  1:1392199/70‑1 (MQ=255)
           atccGTGCACGCACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTAt                                                        >  1:1440171/4‑71 (MQ=255)
              cGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGa                                                     >  1:2506084/1‑71 (MQ=255)
              cGTGCACGCACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGa                                                     <  1:2558874/71‑1 (MQ=255)
               gTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                                    >  1:963958/1‑71 (MQ=255)
               gTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                                    <  1:2972675/71‑1 (MQ=255)
               gTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                                    >  1:2161671/1‑71 (MQ=255)
                  cgcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAccc                                                  <  1:347168/70‑1 (MQ=255)
                  cgcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAccc                                                  <  1:342478/70‑1 (MQ=255)
                  cgcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCa                                                 <  1:968922/71‑1 (MQ=255)
                   gcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAg                                                <  1:2866077/71‑1 (MQ=255)
                     gcACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGAGCCGGTg                                              >  1:2547432/1‑71 (MQ=255)
                     gcACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGAGCCGGTg                                              >  1:2524937/1‑71 (MQ=255)
                     gcACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGAGCCGGTg                                              >  1:2520459/1‑71 (MQ=255)
                      cACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTgg                                             <  1:2958123/71‑1 (MQ=255)
                       aCGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGt                                            <  1:1805484/71‑1 (MQ=255)
                         gAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCt                                          >  1:642806/1‑71 (MQ=255)
                         aaaCATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGag                                                    <  1:1671618/60‑2 (MQ=255)
                          aaCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGa                                                         >  1:378178/1‑55 (MQ=255)
                          aaCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTg                                         <  1:3167200/71‑1 (MQ=255)
                          aaCATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGcacg                                                                 <  1:1086704/47‑1 (MQ=255)
                           aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                                    >  1:3023332/1‑59 (MQ=255)
                           aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGa                                        >  1:1242164/1‑71 (MQ=255)
                           aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGa                                        >  1:602593/1‑71 (MQ=255)
                             aTACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATa                                                                      >  1:2628368/1‑39 (MQ=255)
                              tACCCGTCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGt                                                                >  1:3479220/1‑44 (MQ=255)
                               aCCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATa                                                                      >  1:1083030/1‑37 (MQ=255)
                                   ggCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACgg                                 >  1:1445506/1‑70 (MQ=255)
                                   ggCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACg                                  >  1:3948476/1‑69 (MQ=255)
                                   ggCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGt                                <  1:688501/71‑1 (MQ=255)
                                        cagAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGt                           >  1:3223652/1‑71 (MQ=255)
                                        cagAGTGTGATCCGGGTTCGGGAAGATAGAGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGt                           >  1:12574/1‑71 (MQ=255)
                                          gagTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAAc                             >  1:954253/1‑67 (MQ=255)
                                           agTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGt                           >  1:1175292/1‑68 (MQ=255)
                                           agTGTGATCCAGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCag                        <  1:2877761/71‑1 (MQ=255)
                                            gtgtGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATc                                      <  1:1593406/56‑1 (MQ=255)
                                              gtgATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCaga                       >  1:1108409/1‑69 (MQ=255)
                                               tgATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGa                                                     >  1:1712043/1‑38 (MQ=255)
                                               tgATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAcc                                                   <  1:2816754/40‑1 (MQ=255)
                                                gATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAAt                   >  1:1660834/1‑71 (MQ=255)
                                                 aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCa                                     <  1:3886428/52‑1 (MQ=255)
                                                 aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCa                                     <  1:3697904/52‑1 (MQ=255)
                                                 aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAAtt                  >  1:2336518/1‑71 (MQ=255)
                                                 aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAAt                   <  1:1027793/70‑1 (MQ=255)
                                                   ccGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTc                 >  1:968542/1‑70 (MQ=255)
                                                         tCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGAt                                       >  1:3796927/1‑42 (MQ=255)
                                                         tCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGAt                                       >  1:3545643/1‑42 (MQ=255)
                                                         tCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTa          >  1:3653820/1‑71 (MQ=255)
                                                          cGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTa          <  1:394430/70‑1 (MQ=255)
                                                          cGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTAc         >  1:2428621/1‑71 (MQ=255)
                                                          cGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTAc         >  1:1996961/1‑71 (MQ=255)
                                                            ggAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTg                                         <  1:2761731/37‑1 (MQ=255)
                                                            ggAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCg       <  1:724886/71‑1 (MQ=255)
                                                             gAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGt      >  1:2879162/1‑71 (MQ=255)
                                                               aGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTAc         >  1:2847490/1‑66 (MQ=255)
                                                               aGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTcc    >  1:123678/1‑71 (MQ=255)
                                                               aGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTc     >  1:2122436/1‑70 (MQ=255)
                                                                gATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCt   <  1:1398938/71‑1 (MQ=255)
                                                                 aTAGCGCGTAGGGTGATAGACCCAGTGGTCTGATc                                      <  1:415541/35‑1 (MQ=255)
                                                                 aTAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTg  <  1:661747/71‑1 (MQ=255)
                                                                  |                                                                     
GCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTG  >  minE/313188‑313323

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: