Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 313,239 C→G 10.9% D284H (GAT→CAT) ‡ acrA ← multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE313,2390CG10.9% 166.2 / 10.6 64D284H (GAT→CAT) ‡acrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base G (4/3);  total (33/31)
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.29e-01

AAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTC  >  minE/313170‑313308
                                                                     |                                                                     
aaaTAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATcc                                                                      >  1:2839355/1‑71 (MQ=255)
aaaTAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATcc                                                                      >  1:1643598/1‑71 (MQ=255)
 aaTAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATcc                                                                      <  1:3018348/70‑1 (MQ=255)
      cGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTc                                                                <  1:1249372/71‑1 (MQ=255)
        tttGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCgg                                                              >  1:2870226/1‑71 (MQ=255)
          tGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCggg                                                             <  1:3586547/70‑1 (MQ=255)
           ggATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGa                                                            <  1:1256628/70‑1 (MQ=255)
            gATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGaa                                                           <  1:1484255/70‑1 (MQ=255)
             aTTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCgg                                                              >  1:782222/1‑66 (MQ=255)
             aTTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGa                                                         <  1:2252637/71‑1 (MQ=255)
              ttAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGa                                                         <  1:3160175/70‑1 (MQ=255)
                  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)
                        cttcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATcc                                                                      >  1:2993819/1‑47 (MQ=255)
                        cttcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAg                                               <  1:3262776/70‑1 (MQ=255)
                         ttcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAggg                                             >  1:38203/1‑71 (MQ=255)
                          tcCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGt                                            >  1:1407321/1‑71 (MQ=255)
                           ccAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGaa                                                           <  1:1627026/55‑1 (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:2972675/71‑1 (MQ=255)
                                 gTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                     >  1:963958/1‑71 (MQ=255)
                                 gTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                     >  1:2161671/1‑71 (MQ=255)
                                    cgcgcACGAACATACCCGGCAGCAGAGTGTGATcc                                                                      <  1:3899538/35‑1 (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:2520459/1‑71 (MQ=255)
                                       gcACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGAGCCGGTg                               >  1:2524937/1‑71 (MQ=255)
                                       gcACAAACATACCCGGCAAAAGCGTATGTTGCGGGTTAGGGAAGACAGCACGTAGGGTTATGGAGCCGGTg                               >  1:2547432/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)
                                             aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGa                                                         <  1:1570178/39‑1 (MQ=255)
                                             aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGAc                                     >  1:3023332/1‑59 (MQ=255)
                                             aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGa                         >  1:602593/1‑71 (MQ=255)
                                             aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGa                         >  1:1242164/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)
                                                                     |                                                                     
AAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTC  >  minE/313170‑313308

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