Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 310,794 T→A 8.5% E693D (GAA→GAT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE310,7940TA 8.5% 158.6 / 6.4 59E693D (GAA→GATacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (23/31);  new base A (3/2);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

GTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTG  >  minE/310724‑310864
                                                                      |                                                                      
gTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                                                        <  1:2909431/71‑1 (MQ=255)
gTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                                                        <  1:5607459/71‑1 (MQ=255)
   cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                                                        <  1:4161674/68‑1 (MQ=255)
   cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                                     >  1:3586773/1‑71 (MQ=255)
   cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                                     >  1:3162654/1‑71 (MQ=255)
    gCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                                    <  1:1106219/71‑1 (MQ=255)
    gCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                                    <  1:4482940/71‑1 (MQ=255)
    gCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                                    <  1:5023481/71‑1 (MQ=255)
        gTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                                <  1:295943/71‑1 (MQ=255)
          cAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCGTGAGTCAGTTTTTCGTGAcc                                                                >  1:280878/1‑69 (MQ=255)
          cAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                                <  1:3969444/69‑1 (MQ=255)
          cAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                                >  1:4977237/1‑69 (MQ=255)
             cATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGc                                                           <  1:83119/71‑1 (MQ=255)
             cATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGc                                                           <  1:5482814/71‑1 (MQ=255)
               tatCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                                         <  1:194227/71‑1 (MQ=255)
                      gTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                                <  1:5345261/57‑1 (MQ=255)
                         cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                                               >  1:5678351/1‑71 (MQ=255)
                         cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                                               <  1:4248127/71‑1 (MQ=255)
                         cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCa                                                <  1:2166150/70‑1 (MQ=255)
                         cTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCa                                                <  1:4569373/70‑1 (MQ=255)
                            cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                              >  1:1054075/1‑69 (MQ=255)
                            cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                              >  1:1785531/1‑69 (MQ=255)
                            cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                              >  1:2705732/1‑69 (MQ=255)
                              ctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                                                      <  1:4555066/43‑1 (MQ=255)
                               tgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                                         >  1:3442614/1‑71 (MQ=255)
                                gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGt                                                                    >  1:1755266/1‑43 (MQ=255)
                                gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                                         >  1:1410835/1‑70 (MQ=255)
                                  ttCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                                       >  1:4169904/1‑70 (MQ=255)
                                  ttCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                                       >  1:2339927/1‑70 (MQ=255)
                                   tCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGcc                                                      >  1:1967700/1‑54 (MQ=255)
                                    cTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAc                                                                 <  1:1587239/42‑1 (MQ=255)
                                    cTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa                                     <  1:1220595/70‑1 (MQ=255)
                                    cTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                                      <  1:817166/69‑1 (MQ=255)
                                     tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                                     >  1:3649713/1‑37 (MQ=255)
                                     tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                                     >  1:1733112/1‑37 (MQ=255)
                                     tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                                <  1:1351816/42‑1 (MQ=255)
                                     tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                                <  1:70925/42‑1 (MQ=255)
                                     tGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAcc                                                                <  1:5243633/42‑1 (MQ=255)
                                      gcaagcaaCTTGTTACGCGCCTGAGTCAGTTTTTCGTGACCa                                                               >  1:1807289/1‑42 (MQ=255)
                                      gcaagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCa                                                               >  1:2489477/1‑42 (MQ=255)
                                      ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtc                                  <  1:1626494/70‑1 (MQ=255)
                                        aagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcga                                <  1:999838/71‑3 (MQ=255)
                                        aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaa                                <  1:2247475/71‑1 (MQ=255)
                                        aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaa                                <  1:1283897/71‑1 (MQ=255)
                                         agcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcga                                >  1:1523076/1‑68 (MQ=255)
                                          gcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcgaaa                              >  1:1711755/1‑70 (MQ=255)
                                          gcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcgaaa                              >  1:4088067/1‑70 (MQ=255)
                                          gcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag                              >  1:425769/1‑71 (MQ=255)
                                           caaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc                             >  1:4575506/1‑71 (MQ=255)
                                           caaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGc                             >  1:2708401/1‑71 (MQ=255)
                                            aaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                            <  1:5341281/71‑1 (MQ=255)
                                              cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                                                    <  1:5583754/45‑1 (MQ=255)
                                              cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                              >  1:436058/1‑51 (MQ=255)
                                                 gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                                      <  1:835662/56‑1 (MQ=255)
                                                 gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa                                      <  1:2018286/56‑1 (MQ=255)
                                                 gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                       <  1:4700272/71‑1 (MQ=255)
                                                 gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                       <  1:4797019/71‑1 (MQ=255)
                                                   tACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                       >  1:4778008/1‑69 (MQ=255)
                                                         ccTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTAccc               <  1:577185/71‑1 (MQ=255)
                                                         ccTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTAccc               <  1:3974914/71‑1 (MQ=255)
                                                             aGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGtt           <  1:2198278/71‑1 (MQ=255)
                                                                  gTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtc                                        <  1:3457405/37‑1 (MQ=255)
                                                                   tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                            <  1:2744478/48‑1 (MQ=255)
                                                                   tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                            <  1:3020273/48‑1 (MQ=255)
                                                                   tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                            <  1:4875957/48‑1 (MQ=255)
                                                                   tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                            <  1:383746/48‑1 (MQ=255)
                                                                   tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                       <  1:1198118/53‑1 (MQ=255)
                                                                   tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGc                      <  1:1731220/54‑1 (MQ=255)
                                                                    ttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGa     >  1:2639938/1‑70 (MQ=255)
                                                                     tttCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca                                 >  1:2049266/1‑41 (MQ=255)
                                                                     tttCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca                                 >  1:1446223/1‑41 (MQ=255)
                                                                     tttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCa                     >  1:4645013/1‑53 (MQ=255)
                                                                     tttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGt            >  1:4281216/1‑62 (MQ=255)
                                                                      ttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTTCCCAGTTCCACGATTg  >  1:3998606/1‑71 (MQ=255)
                                                                      ttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGAtt   <  1:2025623/70‑1 (MQ=255)
                                                                      |                                                                      
GTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTG  >  minE/310724‑310864

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: