Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,811,495 G→A L23L (CTG→CTA emrA → multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,811,4950GA89.8% 120.2 / 6.5 49L23L (CTG→CTAemrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (22/22);  total (24/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.87e-01

GCTGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGCTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGACGCATACGTGGCAGGGAATCAAATTCAAAT  >  NC_000913/2811365‑2811617
                                                                                                                                  |                                                                                                                          
gCTGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTActcct                                                                                                                       <  2:338110/136‑1 (MQ=255)
gCTGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTActcct                                                                                                                       >  1:338110/1‑136 (MQ=255)
  tGGCCTTTTTTGCAAGCAAGTCGGCTCAACCGATGAGTTAAGAAAATCCTTGAGAACCATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTtc                                                                                                                  <  1:248761/139‑1 (MQ=255)
   ggCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTtct                                                                                                                 <  1:449373/139‑1 (MQ=255)
    gCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGCTCCTCCTtc                                                                                                                  <  1:848360/137‑1 (MQ=255)
    gCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGCTCCTCCTtc                                                                                                                  >  2:848360/1‑137 (MQ=255)
          ttGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGCTCCTCCTTCTCACCTTg                                                                                                          <  2:910460/139‑1 (MQ=255)
                               ccGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTa                                                                                     >  1:14785/1‑139 (MQ=255)
                                      gTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAg                                                                              <  2:14785/139‑1 (MQ=255)
                                               aTCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACGTTGCTCTTTATAATTATTGCCGTa                                                                                     >  1:267177/1‑123 (MQ=255)
                                               aTCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTCCTCACCTTGCTCTTTATAATTATTGCCGTa                                                                                     <  2:267177/123‑1 (MQ=255)
                                                   tGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGtt                                                                 >  1:247188/1‑139 (MQ=255)
                                                               tGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTg                                                                                                          >  2:446917/1‑86 (MQ=255)
                                                               tGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTg                                                                                                          <  1:446917/86‑1 (MQ=255)
                                                                gAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTg                                                                                                          >  2:279708/1‑85 (MQ=255)
                                                                gAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTg                                                                                                          <  1:279708/85‑1 (MQ=255)
                                                                  gcgcAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCa                                                  >  1:212407/1‑139 (MQ=255)
                                                                    gcAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCAcc                                                                                                             >  1:215688/1‑78 (MQ=255)
                                                                    gcAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCAcc                                                                                                             <  2:215688/78‑1 (MQ=255)
                                                                     cAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTg                                                                                                          <  1:1071304/80‑1 (MQ=255)
                                                                     cAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTg                                                                                                          >  2:1071304/1‑80 (MQ=255)
                                                                     cAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTgg                                                            >  1:333345/1‑126 (MQ=255)
                                                                     cAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTgg                                                            <  2:333345/126‑1 (MQ=255)
                                                                       aaTGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGct                                                                                                        <  2:128397/80‑1 (MQ=255)
                                                                       aaTGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGct                                                                                                        >  1:128397/1‑80 (MQ=255)
                                                                         tGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCgaa                                           <  1:61140/139‑1 (MQ=255)
                                                                           cGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCgaaga                                         >  2:426353/1‑139 (MQ=255)
                                                                             gagaCTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGCTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGa                                                                                 <  2:791885/97‑1 (MQ=255)
                                                                             gagaCTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGCTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGa                                                                                 >  1:791885/1‑97 (MQ=255)
                                                                                aCTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCg                                    <  2:247188/139‑1 (MQ=255)
                                                                                    aaaCCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGa                                >  1:321440/1‑139 (MQ=255)
                                                                                        cccGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGACGCa                            <  2:212407/139‑1 (MQ=255)
                                                                                          cGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTg                                                                    >  2:938562/1‑97 (MQ=255)
                                                                                          cGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTg                                                                    <  1:938562/97‑1 (MQ=255)
                                                                                          cGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGACGCATa                          >  2:975536/1‑139 (MQ=255)
                                                                                            cagcaACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCgaag                                          <  1:834810/121‑1 (MQ=255)
                                                                                            cagcaACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCgaag                                          >  2:834810/1‑121 (MQ=255)
                                                                                            cagcaACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGACGCATACg                        <  1:975536/139‑1 (MQ=255)
                                                                                                           aagaGCGGCAAACGTAAGCGTCTACTCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGAc                               >  2:417666/1‑116 (MQ=255)
                                                                                                           aagaGCGGCAAACGTAAGCGTCTACTCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGAc                               >  2:573635/1‑116 (MQ=255)
                                                                                                           aagaGCGGCAAACGTAAGCGTCTACTCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGAc                               <  1:573635/116‑1 (MQ=255)
                                                                                                           aagaGCGGCAAACGTAAGCGTCTACTCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGAc                               <  1:417666/116‑1 (MQ=255)
                                                                                                            agagCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGACGCATACGTGGCAgg                 <  1:670944/130‑1 (MQ=255)
                                                                                                            agagCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGACGCATACGTGGCAgg                 >  2:670944/1‑130 (MQ=255)
                                                                                                             gagCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGACGCATACGTGGCAGGGAAtcaaatt       >  2:876073/1‑139 (MQ=255)
                                                                                                               gCGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGACGCATACGTGGCAGGGAAtcaaattca     >  1:104790/1‑139 (MQ=255)
                                                                                                                cGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTtat                                                                                                  >  2:548517/1‑45 (MQ=255)
                                                                                                                cGGCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTtat                                                                                                  <  1:548517/45‑1 (MQ=255)
                                                                                                                  gCAAACGTAAGCGTCTACTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGACGCATACGTGGCAGGGAAtcaaattcaaat  <  1:261451/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                          
GCTGGCCTTTTTGGCAAGCAGGTCGGCTCAGCCGATGAGTTAAGAAGATCGTGGAGAACAATATGAGCGCAAATGCGGAGACTCAAACCCCGCAGCAACCGGTAAAGAAGAGCGGCAAACGTAAGCGTCTGCTCCTCCTTCTCACCTTGCTCTTTATAATTATTGCCGTAGCGATAGGGATTTATTGGTTTTTGGTACTGCGTCACTTCGAAGAAACCGATGACGCATACGTGGCAGGGAATCAAATTCAAAT  >  NC_000913/2811365‑2811617

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: