Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,584,306 G→T 16.2% H991N (CAT→AAT) ‡ mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,584,3060GT16.2% 341.3 / 32.7 117H991N (CAT→AAT) ‡mdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base G (46/52);  new base T (7/12);  total (53/64)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 117 total aligned reads displayed.

TGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGC  >  Exported/2584238‑2584373
                                                                    |                                                                   
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:3806177/70‑1 (MQ=255)
    gCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGa                                                                <  1:1577447/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:859642/70‑1 (MQ=255)
        cAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAAc                                                            >  1:2414109/1‑70 (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)
             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)
                            ccGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGaa                                        >  1:1165698/1‑70 (MQ=255)
                            ccGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGaa                                        >  1:531115/1‑70 (MQ=255)
                             cgatacccacaGCGTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAAt                                       >  1:1730401/4‑70 (MQ=9)
                             cGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGAc                                                                   >  1:5354845/1‑42 (MQ=255)
                             cGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAAt                                       >  1:1950990/1‑70 (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:4123503/70‑1 (MQ=9)
                               atacccacaGCGTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATaaa                                    <  1:852092/70‑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:4008447/71‑1 (MQ=255)
                                       cgcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGCCAGCACGCCGAGAATAAAGGCCagag                            <  1:768882/71‑1 (MQ=255)
                                       agcgTTCTGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGcga                               >  1:4694479/2‑66 (MQ=9)
                                        gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa                                     >  1:5171755/1‑61 (MQ=255)
                                        gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa                                     >  1:3908394/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:4521918/50‑1 (MQ=255)
                                           ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                             <  1:4764709/50‑1 (MQ=255)
                                           ttttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                             <  1:4760626/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: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)
                                                cgcgcCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtcat                     <  1:1735021/69‑1 (MQ=255)
                                                cgcgcCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtcat                     <  1:4900481/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:5537871/1‑71 (MQ=255)
                                                   gcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGa                >  1:503638/1‑71 (MQ=255)
                                                    ccAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAt               <  1:162097/71‑1 (MQ=255)
                                                     cAGAACCGGCACCATGACTGATAACCAGAGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAtt              >  1:365942/1‑71 (MQ=255)
                                                        tgccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGAt               <  1:4171865/65‑1 (MQ=12)
                                                         gccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCgaga                             >  1:1120242/2‑52 (MQ=255)
                                                         gccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCgaga                             >  1:2230835/2‑52 (MQ=255)
                                                         gccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCgaga                             >  1:4517874/2‑52 (MQ=255)
                                                          ccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAGGACg         <  1:1165271/71‑1 (MQ=9)
                                                          ccGGCAACGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAGGACg         >  1:735007/1‑71 (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:245938/71‑1 (MQ=255)
                                                            ggCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa       <  1:4030879/71‑1 (MQ=255)
                                                             gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAAt                                       <  1:2535918/38‑1 (MQ=255)
                                                             gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagaga                           <  1:1835747/50‑1 (MQ=255)
                                                             gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagaga                           <  1:5239277/50‑1 (MQ=255)
                                                             gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagaga                           <  1:3166582/50‑1 (MQ=255)
                                                             gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagaga                           <  1:4760639/50‑1 (MQ=255)
                                                             gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagaga                           <  1:4724158/50‑1 (MQ=255)
                                                             gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa       >  1:3605538/1‑70 (MQ=255)
                                                             gCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa       >  1:4600369/1‑70 (MQ=255)
                                                                ccATGACTGATAACCAGCGGCAGCACGCCGAGAATaaa                                    >  1:4262358/1‑38 (MQ=255)
                                                                ccATGACTGATAACCAGCGGCAGCACGCCGAGAATaaa                                    >  1:2335523/1‑38 (MQ=255)
                                                                ccATGACTGATAACCAGCGGCAGCACGCCGAGAATaaa                                    >  1:345050/1‑38 (MQ=255)
                                                                ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagtgag                          <  1:2822000/48‑1 (MQ=255)
                                                                ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                 <  1:4370925/57‑1 (MQ=255)
                                                                ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                 <  1:3692448/57‑1 (MQ=255)
                                                                 cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                 <  1:4535269/56‑1 (MQ=255)
                                                                 cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGa                <  1:5678795/57‑1 (MQ=255)
                                                                 cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGa                <  1:4750877/57‑1 (MQ=255)
                                                                 cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc  <  1:1370591/71‑1 (MQ=255)
                                                                  aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagaga                           >  1:1880683/1‑45 (MQ=255)
                                                                  aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                 >  1:66371/1‑55 (MQ=255)
                                                                  aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTg             >  1:4813082/1‑59 (MQ=255)
                                                                  aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGACAGGTCATCAGGATTGGGCGTaa      >  1:1121840/1‑66 (MQ=255)
                                                                    |                                                                   
TGTTGCGGCAAACATCCCGCCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGC  >  Exported/2584238‑2584373

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: