Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 311,710 A→T 10.6% F388Y (TTC→TAC)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported311,7100AT10.6% 336.3 / 26.8 122F388Y (TTC→TAC) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (57/52);  new base T (7/6);  total (64/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

Only 100 of 131 total aligned reads displayed.

TGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGT  >  Exported/311645‑311780
                                                                 |                                                                      
tGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                                   <  1:1071670/71‑1 (MQ=255)
 ggCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                                   >  1:958283/1‑70 (MQ=255)
  gcatCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCaa                                                                  >  1:1607906/4‑70 (MQ=9)
  gCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                                   >  1:3622721/1‑69 (MQ=255)
    gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                                               >  1:5396240/1‑71 (MQ=255)
    gTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAgg                                                               <  1:2435204/71‑1 (MQ=255)
     tCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggc                                                              <  1:4854868/71‑1 (MQ=255)
     tCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggc                                                              <  1:1420939/71‑1 (MQ=255)
        tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                                                                  <  1:2642586/64‑1 (MQ=255)
        tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                                                                  <  1:5550556/64‑1 (MQ=255)
        tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                           <  1:5510265/71‑1 (MQ=255)
        tCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcgg                                                            >  1:4525407/1‑70 (MQ=255)
        tCCACCAACAGGCAGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                                                                  <  1:4851666/64‑1 (MQ=255)
          caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                           >  1:65230/1‑69 (MQ=255)
          caccaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAg                                                                >  1:2174257/1‑64 (MQ=255)
           accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                                        <  1:2964428/71‑1 (MQ=255)
           accaACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                                        <  1:3421285/71‑1 (MQ=255)
             gagcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAagc                                                              <  1:600496/60‑1 (MQ=9)
             gagcAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGa                                                      <  1:3877372/68‑1 (MQ=12)
             caACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                      >  1:2970675/1‑71 (MQ=255)
              aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGc                                                                     >  1:2718452/1‑55 (MQ=255)
              aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                                                                 >  1:1893672/1‑59 (MQ=255)
              aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaaa                                                                 >  1:497563/1‑59 (MQ=255)
              aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                                   >  1:560383/1‑57 (MQ=255)
              aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCa                                                                   >  1:5318445/1‑57 (MQ=255)
              aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                                     >  1:5483622/1‑71 (MQ=255)
              aaCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGAc                                                     >  1:3244469/1‑71 (MQ=255)
               aCAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAg                                                        <  1:5028315/67‑1 (MQ=255)
                 aGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGc                                                  <  1:841289/71‑1 (MQ=255)
                  ggCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                                 <  1:4191323/71‑1 (MQ=255)
                   cccGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCaa                                                <  1:1437865/70‑1 (MQ=16)
                    ccGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                           <  1:3557547/59‑1 (MQ=255)
                     cGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                              <  1:1094316/71‑1 (MQ=255)
                      gATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                                                                    <  1:82209/48‑1 (MQ=255)
                      gATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                                                                    <  1:1626237/48‑1 (MQ=255)
                      gATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                                                                    <  1:5422774/48‑1 (MQ=255)
                      gATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc                                                                    <  1:3037034/48‑1 (MQ=255)
                       aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                           <  1:4838550/56‑1 (MQ=255)
                       aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                           <  1:3707842/56‑1 (MQ=255)
                       aTGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcggc                                                           <  1:219989/56‑1 (MQ=255)
                       aTGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGt                                            <  1:1344227/71‑1 (MQ=12)
                        tGGCGAGCACCATCCCTAACATTGTTATCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                                 <  1:4981848/65‑1 (MQ=255)
                         ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgaga                                                                        >  1:37014/1‑41 (MQ=255)
                         ggCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                               >  1:5263517/1‑66 (MQ=255)
                          gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                         <  1:3276581/71‑1 (MQ=255)
                          gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                         <  1:4340837/71‑1 (MQ=255)
                          gCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                         <  1:1676517/71‑1 (MQ=255)
                             aGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                                                                  <  1:3260133/43‑1 (MQ=255)
                             aGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCaa                                                                  <  1:646078/43‑1 (MQ=255)
                             aGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAg                                      >  1:2721745/1‑71 (MQ=255)
                             aGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGat                                                     >  1:4923939/1‑55 (MQ=12)
                             aGCACCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTccctaa                                      >  1:1227868/1‑68 (MQ=12)
                               cACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTcc                                          >  1:4627673/1‑65 (MQ=255)
                               cACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                                     >  1:1881803/1‑70 (MQ=255)
                                aCCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                      >  1:785974/1‑52 (MQ=255)
                                aCCATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaa                                   <  1:1302844/71‑1 (MQ=12)
                                 ccATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCa                                    >  1:2964362/1‑69 (MQ=255)
                                 ccATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAAt                                  <  1:4201267/71‑1 (MQ=255)
                                 ccATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaac                                  <  1:1962639/71‑2 (MQ=12)
                                  cATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                                     <  1:1733292/67‑1 (MQ=255)
                                  cATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTc                                           >  1:128890/1‑61 (MQ=12)
                                  cATCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaaca                                 >  1:4166017/1‑69 (MQ=12)
                                    tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAggcgg                                                            >  1:3876650/1‑42 (MQ=255)
                                    tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCaa                                                         >  1:2579569/1‑45 (MQ=255)
                                    tCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                                <  1:4952290/54‑1 (MQ=255)
                                    tCCCGAACATCGTTAGTGTGTTGATGGAGTAACCAAAAGCGGCGAGGATGGCAAACGTCCCTAACaaca                                 >  1:91188/1‑67 (MQ=12)
                                      ccGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                                     <  1:4266106/63‑1 (MQ=255)
                                        gAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                           >  1:4877978/1‑71 (MQ=255)
                                         aaCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                                     <  1:2224202/60‑1 (MQ=255)
                                          aCATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                          >  1:4653480/1‑70 (MQ=255)
                                           cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                      >  1:3814848/1‑41 (MQ=255)
                                           cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGa                                                      >  1:2130522/1‑41 (MQ=255)
                                           cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCa                                                 <  1:1352518/46‑1 (MQ=255)
                                           cATTGTTAGCGTGTTTATCGAGAAGCCAAAGGAGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCg                           >  1:1376820/1‑68 (MQ=255)
                                            aTTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTAc                       <  1:4654934/71‑1 (MQ=255)
                                               gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                         >  1:727980/1‑50 (MQ=255)
                                               gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                         >  1:1032302/1‑50 (MQ=255)
                                               gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                         >  1:4456264/1‑50 (MQ=255)
                                               gTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACgg                     >  1:1724965/1‑70 (MQ=255)
                                                    cGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCaa                                                <  1:800495/38‑1 (MQ=255)
                                                     gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                              >  1:4259160/1‑39 (MQ=255)
                                                     gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                              >  1:2100296/1‑39 (MQ=255)
                                                     gtgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAg                                              >  1:753075/1‑39 (MQ=255)
                                                      tgtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGc                                     <  1:1681563/47‑1 (MQ=255)
                                                       gtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgt                          <  1:4398025/57‑2 (MQ=255)
                                                       gtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                          <  1:3540056/57‑1 (MQ=255)
                                                       gtTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                          <  1:2718617/57‑1 (MQ=255)
                                                        tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                               >  1:5392864/1‑35 (MQ=255)
                                                        tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                               >  1:5407435/1‑35 (MQ=255)
                                                        tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                               >  1:2547344/1‑35 (MQ=255)
                                                        tttATCGAGAAGCCAAAGGCGGCAAGGACGGCaaa                                               >  1:1771258/1‑35 (MQ=255)
                                                        tttATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCgg                          <  1:141612/56‑1 (MQ=255)
                                                        tttATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCgg           <  1:494426/71‑1 (MQ=255)
                                                          tATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTccc                                         <  1:4050109/39‑1 (MQ=255)
                                                            tCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATc       <  1:2303522/71‑1 (MQ=255)
                                                               agaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCa      <  1:208985/69‑1 (MQ=255)
                                                                gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCATCGGTACGGCaa                  >  1:893220/1‑56 (MQ=255)
                                                                gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg   <  1:2329970/71‑1 (MQ=255)
                                                                gaAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACg   <  1:113161/71‑1 (MQ=255)
                                                                 aaGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAAc    >  1:2228641/1‑69 (MQ=255)
                                                                 |                                                                      
TGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAGGCGGCAAGGACGGCAAAGGTCCCGAGCAATACCACCGGTACGGCAATGGTCGGAATCAACGT  >  Exported/311645‑311780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: