Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 313,267 G→T 8.4% I274I (ATC→ATA acrA ← multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE313,2670GT 8.4% 229.3 / 5.0 71I274I (ATC→ATAacrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (36/29);  new base T (4/2);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAAT  >  minE/313198‑313335
                                                                     |                                                                    
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:342478/70‑1 (MQ=255)
        cgcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAccc                                                              <  1:347168/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)
                 aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                                                >  1:3023332/1‑59 (MQ=255)
                 aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGa                                                    >  1:1242164/1‑71 (MQ=255)
                 aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGa                                                    >  1:602593/1‑71 (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)
                                                         aGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAgaga                                 >  1:1565714/1‑50 (MQ=255)
                                                         aGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAgaga                                 >  1:1962884/1‑50 (MQ=255)
                                                          gcgcGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCgg           >  1:3242416/1‑71 (MQ=255)
                                                            gcgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCggg          >  1:2847050/1‑70 (MQ=255)
                                                            gcgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCggg          >  1:2252918/1‑70 (MQ=255)
                                                            gcgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGa         <  1:382082/71‑1 (MQ=255)
                                                             cgTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAg                          <  1:2906554/53‑1 (MQ=255)
                                                                aGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCa                           <  1:1542628/49‑1 (MQ=255)
                                                                 gtgTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTa    <  1:2061526/69‑1 (MQ=255)
                                                                 gggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCg            >  1:645354/1‑63 (MQ=255)
                                                                 gggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACtt     <  1:541315/70‑1 (MQ=255)
                                                                 gggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACtt     <  1:2527133/70‑1 (MQ=255)
                                                                 gggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACtt     >  1:2741882/1‑70 (MQ=255)
                                                                  ggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGc                         >  1:228888/1‑49 (MQ=255)
                                                                  ggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTaa   <  1:2448983/71‑1 (MQ=255)
                                                                   gTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAAt  <  1:2996727/71‑1 (MQ=255)
                                                                     |                                                                    
CAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAAT  >  minE/313198‑313335

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: