Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 481,228 C→A 10.5% intergenic (‑520/+26) tomB ← / ← acrB protein that modulates Hha toxicity/multidrug efflux pump RND permease AcrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913481,2280CA10.5% 110.2 / 7.7 48intergenic (‑520/+26)tomB/acrBprotein that modulates Hha toxicity/multidrug efflux pump RND permease AcrB
Reads supporting (aligned to +/- strand):  ref base C (21/22);  new base A (2/3);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-02

TTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCG  >  NC_000913/481092‑481365
                                                                                                                                        |                                                                                                                                         
ttAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGcc                                                                                                                                     <  1:781360/143‑1 (MQ=255)
  aaTATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtt                                                                                                                                   <  1:450959/143‑1 (MQ=255)
   atatAAATGAACTAACACTCAGTATTTTTTATGAATTTGATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGTGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTa                                                                                                                                  <  2:119586/143‑1 (MQ=255)
    tataAATTAACTAACAATCAGTATTCTCAATGAAGCTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtttt                                                                                                                                 <  1:769355/143‑3 (MQ=255)
         aTTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGAtt                                                                                                                            >  1:542952/1‑143 (MQ=255)
                aacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGtt                                                                                                                     <  1:746067/143‑1 (MQ=255)
                 acAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTg                                                                                                                    >  1:64038/1‑143 (MQ=255)
                           ttttttATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATACTGAGTTGGGGGTTCCATGACGCCTTCATGTACGTAGGTTATGCAGAAAAAAGGCCGCTTACGCGGCCTTCGTGATGACACGTTATAGCAatgtt                                                                                                          >  2:551819/1‑141 (MQ=255)
                                 aTGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGa                                                                                                    <  1:19291/143‑1 (MQ=255)
                                   gAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACa                                                                                                  <  2:147873/143‑1 (MQ=255)
                                           tCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCt                                                                                          <  1:776478/143‑1 (MQ=255)
                                             cGTGGTTACTACTGGTTTGCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtgt                                                                                        <  1:233447/143‑1 (MQ=255)
                                             cGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtgt                                                                                        <  1:800848/143‑1 (MQ=255)
                                               tGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGAGCGACAGGATGGCTGTGc                                                                                      >  1:320711/1‑143 (MQ=255)
                                                 gTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCta                                                                                    >  1:59087/1‑142 (MQ=255)
                                                    aaTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTCTCAATGATGATCGCCAGTCTGGCTCTGCCCGat                                                                                 >  2:258006/1‑143 (MQ=255)
                                                    aaTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGATCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTAGCAATGCTGGCGGACAGTATGGCTCTGCTCgct                                                                                 >  2:272413/1‑141 (MQ=255)
                                                       aCTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGTGTTGGTGATTACACGTTGTGTCAATGATGAACGGCAGTATGGGTGTGGTCGatagc                                                                              >  1:30645/1‑141 (MQ=255)
                                                         tGGTTTTCGTCTGAGATCTTGGGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCtt                                                                            <  2:385271/143‑1 (MQ=255)
                                                          ggTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATCTCTTc                                                                           >  2:673459/1‑143 (MQ=255)
                                                              ttCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGGTCGCAGGTTATGCAGAAAAAAGGCCGCGTACGCGGCCTACGTGATTACACGATGTATCAATGAAGATCGACAGTATGGCTGAGCCCGCTAAGTACATTc                                                                       >  2:725613/1‑143 (MQ=255)
                                                                 gTATGAGATCCTTGGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTg                                                                    <  2:538308/143‑1 (MQ=255)
                                                                  tATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGc                                                                   <  1:271497/143‑1 (MQ=255)
                                                                   aTGAGGTCCTGAGTTGGTGGTTCAATTACTCCTTAATGGTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCg                                                                  <  2:331715/143‑1 (MQ=255)
                                                                    tGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCgg                                                                 <  1:340141/143‑1 (MQ=255)
                                                                    tGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTGACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCgg                                                                 >  1:732085/1‑143 (MQ=255)
                                                                      agaTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCCTCATTCTTGCGGCt                                                               >  2:514692/1‑143 (MQ=255)
                                                                          ccTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAg                                                           >  2:90555/1‑143 (MQ=255)
                                                                                   gtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGACCGATTACGAGACTGTGGGTCTTTCACTTTTTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGGGGCTAAGCGCCGGGGaa                                                 >  2:171861/1‑143 (MQ=255)
                                                                                        tCAATTACTCCGTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAaccac                                             <  2:4955/143‑1 (MQ=255)
                                                                                         cAATTACTCCTTAATGTTCGGATGTTATGCATAAACAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTTCGGTTTAATCAGCCGGCAaccacc                                            <  2:689428/143‑1 (MQ=255)
                                                                                          aaTTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCAATATCTTCATTCTTGCGGCTACAGCGGCGGCGAAccagca                                           >  2:459567/1‑143 (MQ=255)
                                                                                                  cTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCGATGGTGATCGACAGTATGGCTGTGCTCGACATCTTCCTTCTTGCGGCTAAAGAGCAGGCGAGCCACCACGAAGAAt                                   >  2:311929/1‑143 (MQ=255)
                                                                                                    tAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAc                                 <  2:513448/143‑1 (MQ=255)
                                                                                                      gtgtgCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCg                               <  2:125305/142‑1 (MQ=255)
                                                                                                      aTGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCg                               <  1:546953/143‑1 (MQ=255)
                                                                                                         ttCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGAGAACCACCACAAAGAATACCGGaa                            >  1:451052/1‑143 (MQ=255)
                                                                                                         ttCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGCCTCAAGCGGCGGCGAACCACCACCAAGAATCCCGGaa                            >  2:374463/1‑143 (MQ=255)
                                                                                                               ggTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACgaaga                      <  1:197892/143‑1 (MQ=255)
                                                                                                                      cATAAAAGAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGc               <  2:226411/143‑1 (MQ=255)
                                                                                                                        tAAAAAAGGCCGCTTACGCGGCCTTTGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGGGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCCCAAAGAATCCCGGACCGACGGAGATCGCCa             >  2:407425/1‑143 (MQ=255)
                                                                                                                          actaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCAGCCAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGCAGAAGATTGCCAGt           >  1:463119/4‑143 (MQ=255)
                                                                                                                            taaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAGCCACCACAACGACTACCGCAACGCCGACGCTTGCCAGTAc         >  1:746195/2‑143 (MQ=255)
                                                                                                                                 ccGCTTACGCGGCCTTAGTGATGACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATCTCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCCCCACAAAGAATACCGGAACGAGGAAGACGGCCAGAACCgcag    >  2:217898/1‑140 (MQ=255)
                                                                                                                                   gcgtACGCGGCAGGAGTGATTCTGCGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCg  <  2:103633/140‑1 (MQ=255)
                                                                                                                                   gCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGCAGCACGTCGCCAGGACCGCTGCg  >  2:760768/1‑143 (MQ=255)
                                                                                                                                   gCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCg  <  1:90555/143‑1 (MQ=255)
                                                                                                                                   gCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCg  <  1:428460/143‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                         
TTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCG  >  NC_000913/481092‑481365

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