Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,584,304 A→G 16.9% H991H (CAT→CAC) ‡ mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,584,3040AG16.9% 264.8 / 39.1 112H991H (CAT→CAC) ‡mdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base A (39/54);  new base G (7/12);  total (46/66)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

CACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGC  >  minE/2584235‑2584373
                                                                     |                                                                     
cACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCAt                                                                      <  1:2392213/71‑1 (MQ=255)
cACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCAt                                                                      <  1:2206919/71‑1 (MQ=255)
 aCTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCAt                                                                      <  1:8923/70‑1 (MQ=255)
    gTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACt                                                                  >  1:3526394/1‑71 (MQ=255)
     ttGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACt                                                                  >  1:4011988/1‑70 (MQ=255)
       gCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGa                                                                <  1:3806177/70‑1 (MQ=255)
       gCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGa                                                                <  1:1727403/70‑1 (MQ=255)
       gCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGa                                                                <  1:1577447/70‑1 (MQ=255)
         ggCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATaa                                                             <  1:3434867/71‑1 (MQ=255)
          gCAAACATCCCACCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAAc                                                            <  1:1517813/71‑1 (MQ=255)
           cAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAAc                                                            >  1:2414109/1‑70 (MQ=255)
           cAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAAc                                                            <  1:859642/70‑1 (MQ=255)
           cAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAAc                                                            >  1:3157592/1‑70 (MQ=255)
            aaaCATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAAc                                                            >  1:3675681/1‑69 (MQ=255)
            aaaCATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAAc                                                            >  1:2202186/1‑69 (MQ=255)
             aaCATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAg                                                         >  1:1900021/1‑71 (MQ=255)
               cATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGc                                                        <  1:3667268/70‑1 (MQ=255)
                aTCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGa                                                                    >  1:2902715/1‑57 (MQ=255)
                      cccATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCAc                                                >  1:2568325/1‑71 (MQ=255)
                      cccATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCAc                                                >  1:3668171/1‑71 (MQ=255)
                       ccATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACg                                               <  1:765931/71‑1 (MQ=255)
                        cATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGc                                              >  1:5622025/1‑71 (MQ=255)
                          tCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGAt                                                               <  1:238504/52‑1 (MQ=255)
                            aCGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                             >  1:5471326/1‑68 (MQ=255)
                               ccGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGaa                                        >  1:1165698/1‑70 (MQ=255)
                                cgatacccacaGCGTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAAt                                       >  1:1730401/4‑70 (MQ=9)
                                cgatacccacaGCGTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATa                                      <  1:2237258/68‑1 (MQ=9)
                                cGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGAc                                                                   >  1:5354845/1‑42 (MQ=255)
                                cGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAAt                                       >  1:1950990/1‑70 (MQ=255)
                                cGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATa                                      <  1:5205726/71‑1 (MQ=255)
                                cGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATa                                      <  1:2445552/71‑1 (MQ=255)
                                  gTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaaa                                    <  1:1516042/71‑1 (MQ=255)
                                  atacccacaGCGTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATaaa                                    <  1:4123503/70‑1 (MQ=9)
                                  atacccacaGCGTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATaaa                                    <  1:852092/70‑1 (MQ=9)
                                   tacccacaGCGTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAg                                   <  1:4715207/71‑1 (MQ=9)
                                     cctaccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAg                                                         <  1:5339496/47‑1 (MQ=255)
                                     cctaccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGc                                 >  1:3371844/1‑71 (MQ=255)
                                        accGCGTTTTGCGCGCCAGAACCGGCACCATGACTGa                                                                >  1:4548272/1‑37 (MQ=255)
                                        accGCGTTTTGAGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCAc                                                >  1:384676/1‑53 (MQ=255)
                                         ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                               >  1:2544561/1‑69 (MQ=255)
                                         ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                               >  1:1625428/1‑69 (MQ=255)
                                          cgcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagag                            <  1:1895414/71‑1 (MQ=255)
                                          cgcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagag                            <  1:4008447/71‑1 (MQ=255)
                                          cgcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGCCAGCACGCCGAGAATAAAGGCCagag                            <  1:768882/71‑1 (MQ=255)
                                           gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa                                     >  1:3908394/1‑61 (MQ=255)
                                           gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa                                     >  1:5171755/1‑61 (MQ=255)
                                           gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa                                     >  1:1059270/1‑61 (MQ=255)
                                            cgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagagag                          >  1:5291444/1‑71 (MQ=255)
                                             gTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                             <  1:4306235/67‑1 (MQ=255)
                                              ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                             <  1:4764709/50‑1 (MQ=255)
                                              ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                             <  1:4760626/50‑1 (MQ=255)
                                              ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                             <  1:4521918/50‑1 (MQ=255)
                                              ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGc                                 <  1:4992602/62‑1 (MQ=255)
                                              ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                             <  1:562964/66‑1 (MQ=255)
                                              ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                               <  1:4791385/64‑1 (MQ=255)
                                              ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                               <  1:3214057/64‑1 (MQ=255)
                                                  gcgcgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtcat                     <  1:3554749/70‑1 (MQ=255)
                                                   cgcgcCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtcat                     <  1:4289323/69‑1 (MQ=255)
                                                   cgcgcCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtcat                     <  1:129252/69‑1 (MQ=255)
                                                   cgcgcCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtcat                     <  1:1735021/69‑1 (MQ=255)
                                                   cgcgcCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtcat                     <  1:4900481/69‑1 (MQ=255)
                                                   cgcgcCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGgcat                     <  1:4282731/69‑1 (MQ=255)
                                                     cgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                 <  1:3756472/71‑1 (MQ=255)
                                                      gcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGa                >  1:503638/1‑71 (MQ=255)
                                                      gcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGa                >  1:5537871/1‑71 (MQ=255)
                                                        cAGAACCGGCACCATGACTGATAACCAGAGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAtt              >  1:365942/1‑71 (MQ=255)
                                                            gccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCgaga                             >  1:1120242/2‑52 (MQ=255)
                                                            gccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCgaga                             >  1:4517874/2‑52 (MQ=255)
                                                            gccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCgaga                             >  1:2230835/2‑52 (MQ=255)
                                                             ccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAGGACg         <  1:1165271/71‑1 (MQ=9)
                                                              cGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAGGACg         >  1:20626/1‑70 (MQ=9)
                                                               ggCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa       <  1:4758177/71‑1 (MQ=255)
                                                               ggCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa       <  1:5418138/71‑1 (MQ=255)
                                                               ggCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa       <  1:4030879/71‑1 (MQ=255)
                                                                gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAAt                                       <  1:2535918/38‑1 (MQ=255)
                                                                gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagaga                           <  1:4724158/50‑1 (MQ=255)
                                                                gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagaga                           <  1:1835747/50‑1 (MQ=255)
                                                                gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagaga                           <  1:3166582/50‑1 (MQ=255)
                                                                gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagaga                           <  1:5239277/50‑1 (MQ=255)
                                                                gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa       >  1:4600369/1‑70 (MQ=255)
                                                                   ccATGACTGATAACCAGCGGCAGCACGCCGAGAATaaa                                    >  1:345050/1‑38 (MQ=255)
                                                                   ccATGACTGATAACCAGCGGCAGCACGCCGAGAATaaa                                    >  1:688751/1‑38 (MQ=255)
                                                                   ccATGACTGATAACCAGCGGCAGCACGCCGAGAATaaa                                    >  1:4262358/1‑38 (MQ=255)
                                                                   ccATGACTGATAACCAGCGGCAGCACGCCGAGAATaaa                                    >  1:2335523/1‑38 (MQ=255)
                                                                   ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagtgag                          <  1:2822000/48‑1 (MQ=255)
                                                                   ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                 <  1:3692448/57‑1 (MQ=255)
                                                                   ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                 <  1:4370925/57‑1 (MQ=255)
                                                                   ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACg   <  1:3864851/71‑1 (MQ=255)
                                                                    cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                 <  1:4535269/56‑1 (MQ=255)
                                                                    cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGa                <  1:4750877/57‑1 (MQ=255)
                                                                    cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGa                <  1:5678795/57‑1 (MQ=255)
                                                                    cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc  <  1:2464338/71‑1 (MQ=255)
                                                                    cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc  <  1:1370591/71‑1 (MQ=255)
                                                                     aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagaga                           >  1:1880683/1‑45 (MQ=255)
                                                                     aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                               >  1:1698327/1‑41 (MQ=255)
                                                                     aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtc                       <  1:4150350/49‑1 (MQ=255)
                                                                     aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGt                        >  1:3468221/1‑48 (MQ=255)
                                                                     aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                 >  1:66371/1‑55 (MQ=255)
                                                                     aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTg             >  1:4813082/1‑59 (MQ=255)
                                                                     aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGACAGGTCATCAGGATTGGGCGTaa      >  1:1121840/1‑66 (MQ=255)
                                                                     |                                                                     
CACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGC  >  minE/2584235‑2584373

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: