Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,871 T→C 9.7% K334K (AAA→AAG) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8710TC 9.7% 223.5 / 12.7 82K334K (AAA→AAG) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (33/41);  new base C (3/5);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCG  >  minE/311802‑311941
                                                                     |                                                                      
cAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAttt                                                                       <  1:2916420/71‑1 (MQ=255)
 aGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                                                      <  1:2519929/71‑1 (MQ=255)
 aGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                                                      >  1:2720296/1‑71 (MQ=255)
  gATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                                                      <  1:2376201/70‑1 (MQ=255)
  gATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTc                                                                     <  1:16216/71‑1 (MQ=255)
   aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                                                                    >  1:2419244/1‑71 (MQ=255)
      cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                                                                  >  1:3199439/1‑70 (MQ=255)
      cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                                                                  >  1:3619781/1‑70 (MQ=255)
       aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa                                                                <  1:458249/71‑1 (MQ=255)
        tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                                                                  <  1:1422026/68‑1 (MQ=255)
         aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg                                                              >  1:3135286/1‑71 (MQ=255)
         aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg                                                              >  1:3267093/1‑71 (MQ=255)
           ccAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGc                                                            <  1:1785394/71‑1 (MQ=255)
            cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACggag                                                           <  1:1311492/71‑3 (MQ=255)
            cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCg                                                           >  1:1722786/1‑71 (MQ=255)
             aGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                                          <  1:1272663/71‑1 (MQ=255)
             aGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                                          <  1:553765/71‑1 (MQ=255)
              ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                                         <  1:1604772/71‑1 (MQ=255)
              ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                                         <  1:1442507/71‑1 (MQ=255)
                aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt                                                       <  1:2208850/71‑1 (MQ=255)
                  cacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt                                                       <  1:2805917/69‑1 (MQ=255)
                    cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACgg                                                             >  1:392723/1‑61 (MQ=255)
                    cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                   <  1:1542576/71‑1 (MQ=255)
                    cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                   >  1:1162280/1‑71 (MQ=255)
                    cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                   >  1:1150993/1‑71 (MQ=255)
                     gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGt                                                  <  1:2697793/71‑1 (MQ=255)
                     gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGt                                                  <  1:1905565/71‑1 (MQ=255)
                      aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                 >  1:1894549/1‑71 (MQ=255)
                      aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                 <  1:2057586/71‑1 (MQ=255)
                      aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                 >  1:2392995/1‑71 (MQ=255)
                      aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                 <  1:2483885/71‑1 (MQ=255)
                       ggATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                                                                    >  1:2842723/1‑51 (MQ=255)
                       ggATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt                                                <  1:3598567/71‑1 (MQ=255)
                       ggATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt                                                >  1:1509933/1‑71 (MQ=255)
                       ggATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt                                                >  1:609365/1‑71 (MQ=255)
                       ggATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAAGGCGTGGTGTCGTAt                                                >  1:3645292/1‑71 (MQ=255)
                        gatgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATg                                               <  1:3214708/71‑1 (MQ=255)
                            atggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATa                                                 >  1:3762878/4‑65 (MQ=255)
                            atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGt                                            <  1:1875796/70‑1 (MQ=255)
                            atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTa                                           <  1:2937811/71‑1 (MQ=255)
                             tCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTa                                           <  1:85513/70‑1 (MQ=255)
                               gCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAAc                                        <  1:2808217/71‑1 (MQ=255)
                               gCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTagag                                        <  1:1033738/71‑4 (MQ=16)
                               gCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTaga                                         <  1:2746405/70‑3 (MQ=255)
                                cTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGt                                            <  1:2505299/66‑1 (MQ=255)
                                cTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTa                                           <  1:3251511/67‑1 (MQ=16)
                                 ttCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                                         <  1:2039181/52‑1 (MQ=255)
                                  tCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtt                                                                        >  1:586921/1‑36 (MQ=255)
                                  tCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt                                                       >  1:1157676/1‑53 (MQ=255)
                                  tCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGt                                                  >  1:2126767/1‑58 (MQ=255)
                                   cGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACaa                                      <  1:1306529/69‑1 (MQ=255)
                                   cGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACaa                                      <  1:1509709/69‑1 (MQ=255)
                                    gACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAttt                                   <  1:1820918/71‑1 (MQ=255)
                                      ccAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTc                                 <  1:806699/71‑1 (MQ=255)
                                      acAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAgg                                              <  1:2884987/57‑1 (MQ=12)
                                        aGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                  >  1:1361903/1‑68 (MQ=255)
                                         gCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                              >  1:3029372/1‑71 (MQ=255)
                                          cGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAc                             >  1:2309036/1‑71 (MQ=255)
                                           gTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCa                                >  1:1919865/1‑67 (MQ=255)
                                            ttttTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAgg                                              >  1:2817605/1‑52 (MQ=255)
                                              ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATggg                                             >  1:2337074/1‑51 (MQ=255)
                                              ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                              >  1:3565858/1‑66 (MQ=255)
                                              ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGa                         >  1:3821041/1‑71 (MQ=255)
                                                taCCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTa                                           >  1:1690817/2‑51 (MQ=255)
                                                aaCCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGAc                        <  1:266319/70‑1 (MQ=255)
                                                  ccACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACggg                     <  1:2267577/71‑1 (MQ=255)
                                                  ccACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACgg                      >  1:1022388/1‑70 (MQ=255)
                                                    aCTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                 >  1:288281/1‑41 (MQ=255)
                                                     cTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACa                                       >  1:1661374/1‑50 (MQ=255)
                                                     cTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGgaag                  <  1:3082551/71‑1 (MQ=255)
                                                      ttCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                   <  1:3633181/37‑1 (MQ=255)
                                                      ttCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                   <  1:3714020/37‑1 (MQ=255)
                                                       tCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                 >  1:2727521/1‑38 (MQ=255)
                                                        cGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                              >  1:3899604/1‑56 (MQ=255)
                                                        cGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGgaagaa                <  1:3499436/70‑1 (MQ=255)
                                                         gTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAttt                                   <  1:402317/50‑1 (MQ=255)
                                                         gTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAttt                                   <  1:3302175/50‑1 (MQ=255)
                                                          tGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAAc               >  1:1198088/1‑69 (MQ=255)
                                                          tGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACg              <  1:2551473/70‑1 (MQ=255)
                                                           gAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAATCAATTTTCAGa                              >  1:2399734/1‑53 (MQ=255)
                                                           gAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                              >  1:2828827/1‑53 (MQ=255)
                                                            aaTAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTa                                           >  1:2677196/1‑39 (MQ=255)
                                                                gagaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                  <  1:2773972/44‑1 (MQ=255)
                                                                 agaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGgaag                  <  1:979477/59‑1 (MQ=255)
                                                                  gaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCt     <  1:205173/71‑1 (MQ=255)
                                                                   aTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAttt                                   <  1:2456976/40‑1 (MQ=255)
                                                                   aTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCtt    <  1:846757/71‑1 (MQ=255)
                                                                    ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACggg                     <  1:440361/53‑1 (MQ=255)
                                                                    ttttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCtt    <  1:3107271/70‑1 (MQ=255)
                                                                     tttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAg                               >  1:3763265/1‑42 (MQ=255)
                                                                     tttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCg  >  1:644488/1‑71 (MQ=255)
                                                                     |                                                                      
CAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCG  >  minE/311802‑311941

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: