Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,895 T→A 10.8% P326P (CCA→CCT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8950TA10.8% 172.0 / 12.4 65P326P (CCA→CCTacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (24/34);  new base A (2/5);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCGCAG  >  minE/311825‑311965
                                                                      |                                                                      
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)
     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)
            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)
                              cTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACa                                                               >  1:1661374/1‑50 (MQ=255)
                              cTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGgaag                                          <  1:3082551/71‑1 (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)
                                               ttCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACggg                                             >  1:806478/1‑51 (MQ=255)
                                                   cGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGt                     >  1:1624866/1‑71 (MQ=255)
                                                   cGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGt                     >  1:474727/1‑71 (MQ=255)
                                                    gAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAg                                                       <  1:3096623/36‑1 (MQ=255)
                                                    gAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACgg                                              <  1:3057978/45‑1 (MQ=255)
                                                        ggCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGc                <  1:366013/71‑1 (MQ=255)
                                                            tggtgTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCa                  <  1:338232/65‑1 (MQ=255)
                                                            tggtgTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCa                  <  1:2974299/65‑1 (MQ=255)
                                                            tggtgTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCa                  <  1:1771341/65‑1 (MQ=255)
                                                             ggtgTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTcc                                 <  1:1106588/49‑1 (MQ=255)
                                                              gtgtCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGa           <  1:1678893/70‑1 (MQ=255)
                                                               tgtCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCAt                               >  1:951494/1‑49 (MQ=255)
                                                               tgtCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGa           >  1:3237093/1‑69 (MQ=255)
                                                               tgtCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGa           >  1:2083971/1‑69 (MQ=255)
                                                                gtCGTATGGGTAAACAATTTTCAGACCCGACGGgaag                                          <  1:2888826/37‑1 (MQ=255)
                                                                   gTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACgg                                     >  1:681618/1‑39 (MQ=255)
                                                                    tATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACgg            <  1:565006/63‑1 (MQ=255)
                                                                    tATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACgg            <  1:2300513/63‑1 (MQ=255)
                                                                     aTGGGTAAACAATTTTCAGACCCGACGGGAAGAACgg                                     >  1:3137999/1‑37 (MQ=255)
                                                                      tGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTGCATCTTCGCCAGTTCAGCACGGATTGCCGCAg  >  1:2817060/1‑71 (MQ=255)
                                                                      tGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCa               <  1:3785352/58‑1 (MQ=255)
                                                                      tGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCg     >  1:3380831/1‑68 (MQ=255)
                                                                      tGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCGCa   <  1:2292815/70‑1 (MQ=255)
                                                                      |                                                                      
GGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCGCAG  >  minE/311825‑311965

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: