Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,584,295 A→G 13.8% G994G (GGT→GGC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,584,2950AG13.8% 210.9 / 22.3 87G994G (GGT→GGCmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base A (39/36);  new base G (2/10);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 2.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

AAATTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATT  >  minE/2584225‑2584361
                                                                      |                                                                  
aaaTTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGaa                                                                    >  1:2218682/1‑71 (MQ=255)
 aaTTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAAc                                                                   >  1:2807942/1‑71 (MQ=255)
  aTTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAAcc                                                                  <  1:956389/71‑1 (MQ=255)
       cAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCg                                                                 >  1:4007632/1‑67 (MQ=255)
         gCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCTCCAGAACCGGCa                                                              >  1:1079461/1‑68 (MQ=255)
         gCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAAc                                                                   >  1:1213764/1‑63 (MQ=255)
         gCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCa                                                              >  1:968101/1‑68 (MQ=255)
          cACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCAt                                                          <  1:2392213/71‑1 (MQ=255)
          cACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCAt                                                          <  1:2206919/71‑1 (MQ=255)
           aCTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCAt                                                          <  1:8923/70‑1 (MQ=255)
             tGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGa                                                        >  1:3745073/1‑70 (MQ=255)
              gTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGAc                                                       >  1:2874533/1‑70 (MQ=255)
              gTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACt                                                      >  1:3526394/1‑71 (MQ=255)
               ttGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACt                                                      >  1:4011988/1‑70 (MQ=255)
               ttGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTg                                                     <  1:5018544/71‑1 (MQ=255)
                 gCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGa                                                    <  1:1577447/70‑1 (MQ=255)
                 gCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGa                                                    <  1:3806177/70‑1 (MQ=255)
                 gCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGa                                                    <  1:1727403/70‑1 (MQ=255)
                  cGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGAt                                                   <  1:1188193/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:2202186/1‑69 (MQ=255)
                      aaaCATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAAc                                                >  1:3675681/1‑69 (MQ=255)
                       aaCATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAg                                             >  1:1900021/1‑71 (MQ=255)
                        aCATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGc                                            <  1:2005434/71‑1 (MQ=255)
                         cATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGc                                            <  1:3667268/70‑1 (MQ=255)
                         cATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCg                                           >  1:340559/1‑71 (MQ=255)
                          aTCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAAc                                                                   >  1:636662/1‑46 (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)
                                         ccGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGaa                            >  1:531115/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:2445552/71‑1 (MQ=255)
                                          cGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATa                          <  1:5205726/71‑1 (MQ=255)
                                            gTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaaa                        <  1:1516042/71‑1 (MQ=255)
                                            gTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATa                          >  1:2285179/1‑69 (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:1625428/1‑69 (MQ=255)
                                                   ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                   >  1:2544561/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)
                                                    agcgTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGcga                   >  1:4694479/2‑66 (MQ=9)
                                                     gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa                         >  1:1059270/1‑61 (MQ=255)
                                                     gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa                         >  1:3908394/1‑61 (MQ=255)
                                                     gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa                         >  1:5171755/1‑61 (MQ=255)
                                                      cgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagagag              >  1:5291444/1‑71 (MQ=255)
                                                       gTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                 <  1:4306235/67‑1 (MQ=255)
                                                        ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                 <  1:4521918/50‑1 (MQ=255)
                                                        ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                 <  1:4760626/50‑1 (MQ=255)
                                                        ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                 <  1:4764709/50‑1 (MQ=255)
                                                        ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGc                     <  1:4992602/62‑1 (MQ=255)
                                                        ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                 <  1:5564802/66‑1 (MQ=255)
                                                        ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                 <  1:562964/66‑1 (MQ=255)
                                                        ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                   <  1:4791385/64‑1 (MQ=255)
                                                        ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                   <  1:3214057/64‑1 (MQ=255)
                                                           tGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtcat         >  1:391086/1‑71 (MQ=255)
                                                            gcgcgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtcat         <  1:3554749/70‑1 (MQ=255)
                                                             cgcgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtcatc        <  1:1793153/70‑1 (MQ=255)
                                                             cgcgcCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtcatc        <  1:4469410/70‑1 (MQ=255)
                                                             cgcgcCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtcat         <  1:4900481/69‑1 (MQ=255)
                                                             cgcgcCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtcat         <  1:1735021/69‑1 (MQ=255)
                                                             cgcgcCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtcat         <  1:129252/69‑1 (MQ=255)
                                                             cgcgcCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtcat         <  1:4289323/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)
                                                                 ccAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAt   <  1:162097/71‑1 (MQ=255)
                                                                  cAGAACCGGCACCATGACTGATAACCAGAGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAtt  >  1:365942/1‑71 (MQ=255)
                                                                      |                                                                  
AAATTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATT  >  minE/2584225‑2584361

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: