Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,840 C→A 17.8% V345F (GTC→TTC)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8400CA17.8% 367.9 / 40.8 134V345F (GTC→TTC) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (54/56);  new base A (16/8);  total (70/64)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 135 total aligned reads displayed.

AATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTT  >  minE/311772‑311909
                                                                    |                                                                     
aaTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                                                                     <  1:5031656/71‑1 (MQ=255)
 aTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAg                                                                    >  1:1071591/1‑71 (MQ=255)
 aTCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAg                                                                    >  1:2390574/1‑71 (MQ=255)
    aaCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                                                                  <  1:3908644/70‑1 (MQ=255)
      cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                                                               >  1:5371946/1‑71 (MQ=255)
      cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                                                               >  1:2754883/1‑71 (MQ=255)
      cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                                                               >  1:2748303/1‑71 (MQ=255)
      cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                                                               >  1:1417996/1‑71 (MQ=255)
       gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                                              <  1:950230/71‑1 (MQ=255)
       gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                                              <  1:4426564/71‑1 (MQ=255)
       gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                                              <  1:4208913/71‑1 (MQ=255)
       gTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                                              <  1:1484380/71‑1 (MQ=255)
        tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                                                                  <  1:1937636/66‑1 (MQ=255)
        tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                                             <  1:3918965/71‑1 (MQ=255)
        tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                                             <  1:1121907/71‑1 (MQ=255)
        tCGCGCGGAAGTTCTGAAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                                                                  <  1:2248828/66‑1 (MQ=255)
          gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                            >  1:3370189/1‑70 (MQ=255)
          gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                           <  1:2932654/71‑1 (MQ=255)
             cgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCAc                                                        <  1:4399059/71‑1 (MQ=255)
              ggAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACt                                                       >  1:3698759/1‑71 (MQ=255)
                  gTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGt                                                   >  1:3065788/1‑71 (MQ=255)
                   ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                                                                     <  1:4535045/52‑1 (MQ=255)
                   ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                                                                     <  1:5266851/52‑1 (MQ=255)
                   ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCa                                                                     <  1:4895700/52‑1 (MQ=255)
                   ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                                  <  1:1266517/71‑1 (MQ=255)
                   ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                                  <  1:4650009/71‑1 (MQ=255)
                   ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                                  <  1:1448128/71‑1 (MQ=255)
                   ttCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGt                                                   >  1:4502896/1‑70 (MQ=18)
                    tCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                                  >  1:4024720/1‑70 (MQ=255)
                    tCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTg                                                  >  1:3099907/1‑70 (MQ=18)
                     cTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGc                                                                   >  1:2516665/1‑52 (MQ=255)
                     cTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                                                >  1:2686001/1‑71 (MQ=255)
                        cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAAt                                               <  1:3558846/69‑1 (MQ=255)
                        cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAAt                                               <  1:207952/69‑1 (MQ=255)
                         aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                            <  1:836204/71‑1 (MQ=255)
                         aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                            <  1:1660673/71‑1 (MQ=255)
                             aCAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCa                                                         <  1:2979806/54‑1 (MQ=255)
                             aCAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtt                                        >  1:4570445/1‑71 (MQ=255)
                                 aTACATAACGAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                                    >  1:1882947/1‑71 (MQ=255)
                                 aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                                                <  1:1074778/59‑1 (MQ=255)
                                   aCATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGttttt                                                             <  1:5185825/44‑2 (MQ=16)
                                   aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCAc                                   >  1:4986943/1‑70 (MQ=255)
                                    cATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGttttt                                                             <  1:412807/43‑2 (MQ=16)
                                    cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGa                                 <  1:2122378/71‑1 (MQ=255)
                                    cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGa                                 <  1:4763334/71‑1 (MQ=255)
                                     aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                                                >  1:3922039/1‑55 (MQ=255)
                                     aTAACCAGGAACACGAGGATGATCCCTTCGACCAGCGTTTTAACCACTTCGTGaa                                                >  1:1751603/1‑55 (MQ=255)
                                      tcaCCAGGAACACCAGCATAATGGCTTCGAACAGCGtttt                                                              <  1:3807416/38‑1 (MQ=16)
                                      tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAc                               <  1:4887055/71‑1 (MQ=255)
                                       aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGtttt                                                              >  1:620188/1‑39 (MQ=255)
                                       aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                            >  1:3942716/1‑57 (MQ=255)
                                       aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                            >  1:5133547/1‑57 (MQ=255)
                                       aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                            >  1:3212469/1‑57 (MQ=255)
                                       aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                            >  1:1121544/1‑57 (MQ=255)
                                       aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATaga                                            >  1:573596/1‑57 (MQ=255)
                                        aacAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                            >  1:4258127/3‑40 (MQ=255)
                                        aCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                            >  1:3716009/1‑40 (MQ=255)
                                        aCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                            >  1:4075661/1‑40 (MQ=255)
                                        aCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                            >  1:2845680/1‑40 (MQ=255)
                                          cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                                             <  1:4088301/37‑1 (MQ=255)
                                          cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                                             <  1:163417/37‑1 (MQ=255)
                                           aGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                          <  1:5056854/71‑1 (MQ=255)
                                           aGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                          <  1:5323838/71‑1 (MQ=255)
                                            ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                         >  1:2177368/1‑71 (MQ=255)
                                            ggAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtg                         <  1:3748263/71‑1 (MQ=255)
                                             gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                      <  1:1019905/57‑1 (MQ=255)
                                             gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                      <  1:326925/57‑1 (MQ=255)
                                             gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAt                                         >  1:2514016/1‑54 (MQ=255)
                                             gAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATag                                             <  1:4688178/50‑1 (MQ=255)
                                             gAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtgg                        >  1:4428081/1‑71 (MQ=255)
                                              aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt                       <  1:2953439/71‑1 (MQ=255)
                                              aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                         >  1:3123258/1‑69 (MQ=255)
                                              aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                         >  1:4267877/1‑69 (MQ=255)
                                              aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATATTCACGAACGGCgtg                         >  1:3374507/1‑69 (MQ=255)
                                               acacGAGGATGATCGCTTCGACCAGCGTTTTAACCa                                                         <  1:4256686/36‑1 (MQ=255)
                                               acacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggtg                      >  1:4649374/1‑71 (MQ=255)
                                                 acGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGgtgt                     <  1:1719512/70‑1 (MQ=255)
                                                  cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                      >  1:1144417/1‑52 (MQ=255)
                                                  cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTc                                     >  1:4707099/1‑53 (MQ=255)
                                                  cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                                    >  1:1715554/1‑54 (MQ=255)
                                                    aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                                    >  1:1313830/1‑52 (MQ=255)
                                                    aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGAGTGGTGTCGTa                 >  1:2278972/1‑71 (MQ=255)
                                                        tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg                              <  1:2902771/54‑1 (MQ=255)
                                                        tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                          >  1:4893259/1‑58 (MQ=255)
                                                        tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATggg             >  1:3992616/1‑71 (MQ=255)
                                                         gatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                                  <  1:2354159/49‑1 (MQ=255)
                                                          atggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAggg             >  1:3219919/4‑69 (MQ=12)
                                                          atggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAggg             >  1:3060770/4‑69 (MQ=12)
                                                          atggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTa           <  1:1614961/68‑1 (MQ=16)
                                                          atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACTGCGTGGTGTCGTATGGGTa           <  1:2824688/71‑1 (MQ=255)
                                                          atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTa           <  1:3785131/71‑1 (MQ=255)
                                                           tggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAggg             >  1:132925/3‑68 (MQ=16)
                                                             gCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCgtg                         >  1:614230/1‑54 (MQ=255)
                                                             gCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTa           >  1:4444151/1‑68 (MQ=16)
                                                               ttCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACaa      >  1:3547085/1‑71 (MQ=255)
                                                                tCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAt     >  1:3574996/1‑71 (MQ=255)
                                                                 cGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACaa      <  1:2911483/69‑1 (MQ=255)
                                                                   aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt  <  1:4342279/71‑1 (MQ=255)
                                                                   aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt  <  1:2827703/71‑1 (MQ=255)
                                                                   aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAttt   >  1:5557755/1‑70 (MQ=255)
                                                                    |                                                                     
AATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTT  >  minE/311772‑311909

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