Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,584,283 T→C 8.2% Q998Q (CAA→CAG mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,584,2830TC 8.2% 203.7 / 6.1 73Q998Q (CAA→CAGmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base T (41/26);  new base C (2/4);  total (43/30)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGAAGTAAATTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAG  >  Exported/2584217‑2584349
                                                                  |                                                                  
aaCGAAGTAAATTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTgc                                                                 >  1:5464394/1‑70 (MQ=255)
 aCGAAGTAAATTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTgcgc                                                               <  1:4191195/71‑1 (MQ=255)
 aCGAAGTAAATTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTgcgc                                                               <  1:313613/71‑1 (MQ=255)
    aaGTAAATTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTgcgcgc                                                             >  1:1203647/1‑70 (MQ=255)
        aaaTTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTg                                                                  >  1:2673245/1‑61 (MQ=255)
        aaaTTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGaa                                                        >  1:2218682/1‑71 (MQ=255)
         aaTTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTgcg                                                                >  1:3031471/1‑62 (MQ=255)
         aaTTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAAc                                                       >  1:2807942/1‑71 (MQ=255)
          aTTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAAcc                                                      <  1:956389/71‑1 (MQ=255)
               cAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCg                                                     >  1:4007632/1‑67 (MQ=255)
                 gCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGtttt                                                                   >  1:4456243/1‑51 (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:2206919/71‑1 (MQ=255)
                  cACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCAt                                              <  1:2392213/71‑1 (MQ=255)
                   aCTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCAt                                              <  1:8923/70‑1 (MQ=255)
                     tGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGa                                            >  1:3745073/1‑70 (MQ=255)
                      gTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCac                                                          >  1:4058827/1‑54 (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)
                        tGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGtttt                                                                   >  1:4769934/1‑44 (MQ=255)
                         gCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGa                                        <  1:3806177/70‑1 (MQ=255)
                         gCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGa                                        <  1:1727403/70‑1 (MQ=255)
                         gCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGa                                        <  1:1577447/70‑1 (MQ=255)
                          cGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGAt                                       <  1:1188193/70‑1 (MQ=255)
                           ggCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATaa                                     <  1:3434867/71‑1 (MQ=255)
                            gCAAACATCCCGCCCATCACGCCGGTACCTACCGCGttt                                                                    <  1:869946/39‑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)
                                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:3668171/1‑71 (MQ=255)
                                        cccATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCAc                        >  1:2568325/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:531115/1‑70 (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)
                                                    gTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATa              >  1:2285179/1‑69 (MQ=255)
                                                    atacccacaGCGTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATaaa            <  1:852092/70‑1 (MQ=9)
                                                    atacccacaGCGTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATaaa            <  1:4123503/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)
                                                            agcgTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGcga       >  1:4694479/2‑66 (MQ=9)
                                                             gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa             >  1:1059270/1‑61 (MQ=255)
                                                             gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa             >  1:5171755/1‑61 (MQ=255)
                                                             gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa             >  1:3908394/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: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)
                                                                  |                                                                  
AACGAAGTAAATTGCCAGCACTGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAG  >  Exported/2584217‑2584349

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: