Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,529 A→G 11.7% M115T (ATG→ACG)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5290AG11.7% 207.2 / 17.4 77M115T (ATG→ACG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (37/31);  new base G (6/3);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAG  >  minE/312459‑312599
                                                                      |                                                                      
cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa                                                                        >  1:3013521/1‑71 (MQ=255)
cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa                                                                        <  1:3217671/71‑1 (MQ=255)
cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa                                                                        <  1:129090/71‑1 (MQ=255)
  tCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATc                                                                      <  1:4006811/71‑1 (MQ=255)
  tCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATc                                                                      <  1:339736/71‑1 (MQ=255)
      gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                                    >  1:1051993/1‑69 (MQ=255)
      gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                                    >  1:2661605/1‑69 (MQ=255)
      gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                                    >  1:2656583/1‑69 (MQ=255)
      gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                                    >  1:2771831/1‑69 (MQ=255)
      gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                                    >  1:3842665/1‑69 (MQ=255)
      gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                                    >  1:3920429/1‑69 (MQ=255)
        aGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGccct                                                                 >  1:2937306/1‑68 (MQ=255)
          ctgctgTATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg                                                              <  1:2824893/71‑1 (MQ=255)
          ctgctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg                                                              <  1:1574427/71‑1 (MQ=255)
                  tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                                       <  1:2469988/70‑1 (MQ=255)
                  tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                                       <  1:46484/70‑1 (MQ=255)
                  tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt                                                      <  1:263968/71‑1 (MQ=255)
                  tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt                                                      <  1:2443771/71‑1 (MQ=255)
                    atTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGtt                                                     <  1:3458726/70‑1 (MQ=255)
                       tctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCAt                                                                       >  1:1570324/1‑49 (MQ=255)
                       tctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTg                                                  <  1:3633194/70‑1 (MQ=255)
                        ctcAACGCTCACCCCTTGCTGCTGGACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCa                                                >  1:3768052/1‑71 (MQ=255)
                        ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc                                                 <  1:1285477/70‑1 (MQ=255)
                        ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc                                                 <  1:2309462/70‑1 (MQ=255)
                        ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCa                                                >  1:3244722/1‑71 (MQ=255)
                           aaCGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCt                                             >  1:737035/1‑71 (MQ=255)
                            aCGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTg                                            <  1:2484024/71‑1 (MQ=255)
                             cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg                                                              >  1:1001038/1‑52 (MQ=255)
                             cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                                           <  1:3679459/71‑1 (MQ=255)
                                 gatccccTGCTGCTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCAcc                                              >  1:2370368/4‑64 (MQ=255)
                                 cACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGtt                                                     >  1:2061855/1‑57 (MQ=255)
                                   ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                     >  1:2161056/1‑71 (MQ=255)
                                   ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                     >  1:3584377/1‑71 (MQ=255)
                                   ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                     <  1:3909276/71‑1 (MQ=255)
                                      cTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                                   >  1:144054/1‑70 (MQ=255)
                                       ttGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGttt                                                        <  1:2385447/48‑1 (MQ=255)
                                        tgctgcTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                                       <  1:983296/48‑1 (MQ=255)
                                        tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTg                                                       <  1:2691196/48‑1 (MQ=19)
                                         gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCt                                                   >  1:2949376/1‑51 (MQ=255)
                                         gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                <  1:1439362/70‑1 (MQ=255)
                                         gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACATGCGCGata                               >  1:4121072/1‑71 (MQ=255)
                                         gctgctGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGt                                                          >  1:4055741/1‑44 (MQ=16)
                                         gctgctGAACATCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGt                                                          >  1:1610271/1‑44 (MQ=16)
                                          ctgctgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGat                                <  1:3652408/69‑1 (MQ=255)
                                            gctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATcc                            >  1:3016428/1‑71 (MQ=255)
                                             ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                         >  1:3629191/1‑57 (MQ=255)
                                             ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                <  1:882196/66‑1 (MQ=255)
                                             ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg                           >  1:2171676/1‑71 (MQ=255)
                                               gAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                                      <  1:510482/58‑1 (MQ=255)
                                                aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                          <  1:1552607/69‑1 (MQ=255)
                                                aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                        <  1:3180559/71‑1 (MQ=255)
                                                 aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCAc                                               >  1:3997405/1‑47 (MQ=255)
                                                 aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                        >  1:4026884/1‑70 (MQ=255)
                                                  cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                >  1:1272280/1‑61 (MQ=255)
                                                  cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                >  1:1338591/1‑61 (MQ=255)
                                                  cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                >  1:2182012/1‑61 (MQ=255)
                                                   ttcttGCGGCAGCAACGGCATCGCCAGCTGCAGtt                                                         <  1:3660103/35‑1 (MQ=255)
                                                    tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                     >  1:35390/1‑70 (MQ=255)
                                                     cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                         >  1:1240221/1‑49 (MQ=255)
                                                     cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                         >  1:1235457/1‑49 (MQ=255)
                                                     cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                   >  1:92574/1‑71 (MQ=255)
                                                        gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                  >  1:1210264/1‑69 (MQ=255)
                                                        gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                  >  1:2360078/1‑69 (MQ=255)
                                                        gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                  >  1:3028572/1‑69 (MQ=255)
                                                         cGGCATCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                <  1:737645/70‑2 (MQ=12)
                                                         cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                <  1:1260947/70‑1 (MQ=255)
                                                           gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc             <  1:304415/71‑1 (MQ=255)
                                                              gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg           >  1:3535914/1‑69 (MQ=12)
                                                              gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg           >  1:3531102/1‑69 (MQ=12)
                                                              gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg           >  1:2110872/1‑69 (MQ=12)
                                                              gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGata                               <  1:1416572/50‑1 (MQ=255)
                                                              gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa          <  1:3570354/71‑1 (MQ=255)
                                                              gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa          <  1:2788957/71‑1 (MQ=255)
                                                               caACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa         <  1:3336376/71‑1 (MQ=255)
                                                                 aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg       >  1:2376829/1‑71 (MQ=255)
                                                                 aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg       >  1:2351372/1‑71 (MQ=255)
                                                                  cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTTAACCTGCGCGATATCCGCAt                        <  1:791388/53‑1 (MQ=255)
                                                                   gggaTCGACAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt     <  1:3347239/68‑1 (MQ=255)
                                                                   ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGa              <  1:2990753/62‑1 (MQ=255)
                                                                     cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt     >  1:2012322/1‑69 (MQ=255)
                                                                     cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt     >  1:3728487/1‑69 (MQ=255)
                                                                      aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg  <  1:2178774/71‑1 (MQ=255)
                                                                      aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg  <  1:841145/71‑1 (MQ=255)
                                                                      aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg  <  1:1578154/71‑1 (MQ=255)
                                                                      |                                                                      
CATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAG  >  minE/312459‑312599

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: