Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,529 T→C 5.3% V448V (GTA→GTG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,5290TC 5.3% 278.0 / 5.6 95V448V (GTA→GTGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (52/38);  new base C (3/2);  total (55/40)
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

GAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCA  >  minE/311461‑311599
                                                                    |                                                                      
gagaACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCAc                                                                      >  1:2879698/1‑71 (MQ=9)
 agaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAcc                                                                     <  1:2908535/71‑1 (MQ=255)
     cTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATcg                                                                 >  1:3074943/1‑71 (MQ=255)
       gACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                                                                >  1:2217854/1‑70 (MQ=255)
       gACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                                                                >  1:2164353/1‑70 (MQ=255)
        gcGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGa                                                              <  1:57355/70‑1 (MQ=255)
        aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                                                                <  1:2122659/69‑1 (MQ=255)
        aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                                                                <  1:135907/69‑1 (MQ=255)
        aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                                                                <  1:3675328/69‑1 (MQ=255)
        aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                                                                <  1:2649993/69‑1 (MQ=255)
             agatagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATAcc                                                          >  1:203348/1‑70 (MQ=255)
               atagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                                                                <  1:440967/62‑1 (MQ=255)
               atagCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                                                                <  1:136807/62‑1 (MQ=255)
                 agCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGAc                                                       >  1:1004187/1‑69 (MQ=255)
                  gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGa                                                              >  1:2320443/1‑61 (MQ=255)
                  gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCag                                                    <  1:3208573/71‑1 (MQ=255)
                  gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCag                                                    <  1:1857996/71‑1 (MQ=255)
                  gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCag                                                    <  1:304715/71‑1 (MQ=255)
                    aCCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAcc                                                                     <  1:1521534/52‑1 (MQ=255)
                       aGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAcc                                                                     >  1:3144201/1‑49 (MQ=255)
                       aGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                                                                      >  1:81596/1‑48 (MQ=255)
                       aGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAgcgc                                               >  1:3139594/1‑71 (MQ=255)
                         tAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                                                                <  1:3356690/52‑1 (MQ=255)
                          aGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGccc                                             >  1:3841881/1‑70 (MQ=255)
                          aGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt                                            >  1:1238850/1‑71 (MQ=255)
                           gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgc                                                                <  1:3635217/50‑1 (MQ=255)
                           gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGa                                                              >  1:2131698/1‑52 (MQ=255)
                           gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCa                                                     >  1:551963/1‑61 (MQ=255)
                           gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGccc                                             >  1:3752865/1‑69 (MQ=255)
                           gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt                                            >  1:3172011/1‑70 (MQ=255)
                           gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt                                            >  1:3295307/1‑70 (MQ=255)
                           gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt                                            >  1:1065092/1‑70 (MQ=255)
                            aaCCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa                                          <  1:2121909/71‑1 (MQ=255)
                             aCCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa                                          <  1:2506293/70‑1 (MQ=255)
                             aCCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa                                          <  1:907978/70‑1 (MQ=255)
                             aCCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa                                          <  1:307462/70‑1 (MQ=255)
                              ccGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAgcgc                                               >  1:2925843/1‑64 (MQ=255)
                              ccGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAgcgc                                               >  1:3389224/1‑64 (MQ=255)
                              ccGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAgcgc                                               >  1:3394083/1‑64 (MQ=255)
                               cGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATc                                       <  1:2007500/71‑1 (MQ=255)
                                gCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                                                                      >  1:2684275/1‑39 (MQ=255)
                                gCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGaa                                         >  1:1976328/1‑68 (MQ=255)
                                  cAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCAt                                                                   >  1:859337/1‑40 (MQ=255)
                                   aaaGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt                                            >  1:917748/1‑62 (MQ=255)
                                   aaaGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa                                          <  1:1683033/64‑1 (MQ=255)
                                   aaaGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTccc                                   <  1:1776318/71‑1 (MQ=255)
                                    aagaagGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGAc                                                       >  1:2206760/1‑50 (MQ=255)
                                     agaagGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATAcc                                                          >  1:3513594/1‑46 (MQ=255)
                                     agaagGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATAcc                                                          >  1:2018142/1‑46 (MQ=255)
                                     agaagGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCa                                                     <  1:1153421/51‑1 (MQ=255)
                                     agaagGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGaa                                         <  1:2970877/63‑1 (MQ=255)
                                       aagGCCATCGGTACGAATACCGCCGACAGTACCATCGCGa                                                              >  1:718917/1‑40 (MQ=255)
                                         tgCCATCGGAATAAATACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCg                             <  1:1810850/70‑1 (MQ=255)
                                         ggCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCg                             <  1:3198546/71‑1 (MQ=255)
                                          gCCATCGGTACGAATACCGCCGACAGTACCATcgc                                                                <  1:3344247/35‑1 (MQ=255)
                                          gCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCACATCGa                            <  1:942695/71‑1 (MQ=255)
                                           ccATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCaga                                                   >  1:3116835/1‑47 (MQ=255)
                                           ccATCGGAATAAATACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGat                           >  1:3998371/1‑70 (MQ=255)
                                            cATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGccc                                             <  1:1901301/52‑1 (MQ=255)
                                            cATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt                                            >  1:282216/1‑53 (MQ=255)
                                             aTCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGccc                                             >  1:2053496/1‑51 (MQ=255)
                                             aTCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGccc                                             >  1:3682092/1‑51 (MQ=255)
                                             aTCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACtt                         <  1:2069669/71‑1 (MQ=255)
                                             aTCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACtt                         <  1:3627637/71‑1 (MQ=255)
                                             aTCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACt                          >  1:1312938/1‑70 (MQ=255)
                                                ggTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATc                                       >  1:2068817/1‑54 (MQ=255)
                                                ggTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACg                      >  1:2697424/1‑71 (MQ=255)
                                                ggTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACg                      >  1:3798475/1‑71 (MQ=255)
                                                ggTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACg                      >  1:3190938/1‑71 (MQ=255)
                                                   aCGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACggg                    >  1:2551911/1‑70 (MQ=255)
                                                   aCGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACg                      >  1:2472916/1‑68 (MQ=255)
                                                     gAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCg                             <  1:3716955/59‑1 (MQ=255)
                                                      aaTACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATc                                       >  1:1409621/1‑48 (MQ=255)
                                                      aaTACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGa                            >  1:2374307/1‑59 (MQ=255)
                                                      aaTACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACtt                         <  1:3926985/62‑1 (MQ=255)
                                                      aaTACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGc                >  1:3307058/1‑71 (MQ=255)
                                                      aaTACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGc                >  1:65617/1‑71 (MQ=255)
                                                       aTACCGCCGACAGTACCATCGCGATACCGCCCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGct               >  1:24267/1‑71 (MQ=255)
                                                        tACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa                                          <  1:131555/43‑1 (MQ=255)
                                                        tACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCCTCGACTTACGGGTAg                 >  1:3256523/1‑68 (MQ=255)
                                                        tACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGct               >  1:3960291/1‑70 (MQ=255)
                                                          ccgccgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAAt                                        >  1:370888/1‑43 (MQ=255)
                                                          ccgccgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGcttct            <  1:3643074/71‑1 (MQ=255)
                                                          ccgccgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGcttct            <  1:3223678/71‑1 (MQ=255)
                                                          ccgccgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGcttc             <  1:326535/70‑1 (MQ=255)
                                                            gccgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCttt          >  1:902879/1‑71 (MQ=255)
                                                             ccgACAGTACCATCGCGATATCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTa                        >  1:404071/1‑56 (MQ=255)
                                                             ccgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTa                        >  1:2314051/1‑56 (MQ=255)
                                                             ccgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTa                        >  1:3809441/1‑56 (MQ=255)
                                                             ccgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTa                        >  1:366257/1‑56 (MQ=255)
                                                              cgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCtttt         <  1:82532/70‑1 (MQ=255)
                                                              cgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCtttt         <  1:834924/70‑1 (MQ=255)
                                                              cgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCtttt         <  1:3426185/70‑1 (MQ=255)
                                                              cgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTg        <  1:2786649/71‑1 (MQ=255)
                                                                 cAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTggcg     >  1:114347/1‑71 (MQ=255)
                                                                    tACCATCGCGATACCGACCAGAGCGCCCTGAATCTACCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCa  <  1:3499010/71‑1 (MQ=255)
                                                                    |                                                                      
GAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCA  >  minE/311461‑311599

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: