Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 485,614 G→T N2K (AAC→AAA acrA ← multidrug efflux pump membrane fusion lipoprotein AcrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913485,6140GT98.6% 256.9 / ‑2.8 73coding (6/1194 nt)acrAmultidrug efflux pump membrane fusion lipoprotein AcrA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (55/17);  minor base . (0/1);  total (55/18)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-01

GTTAGGGCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTGTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAC  >  NC_000913/485552‑485674
                                                               |                                                            
gTTAGGGCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGtttttttt                                                            <  1:2384618/65‑1 (MQ=255)
 ttAGGGCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTc                                                           >  1:726457/1‑65 (MQ=255)
    gggCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCata                                                        <  1:1043241/65‑1 (MQ=255)
     ggCTAAGCTGCCTGAGGGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCatat                                                       >  1:607036/1‑65 (MQ=255)
     ggCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCatat                                                       >  1:335319/1‑65 (MQ=255)
     ggCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCatat                                                       >  1:3242143/1‑65 (MQ=255)
     ggCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCatat                                                       <  1:464574/65‑1 (MQ=255)
     ggCTAAGCTGCCTGAGAGCATCAGAACGACCACCAGAGGCGTAAACCCTCTGTTTTTTTTCatat                                                       >  1:3341834/1‑65 (MQ=255)
       cTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGt                                                     >  1:423717/1‑65 (MQ=255)
       cTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGt                                                     >  1:1862953/1‑65 (MQ=255)
         aaGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTCATATGTaaa                                                  <  1:3873217/65‑1 (MQ=255)
          aGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTaaa                                                  <  1:2965940/65‑1 (MQ=255)
          aGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTaaa                                                  <  1:3880955/65‑1 (MQ=255)
           gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTacac                                                 >  1:2403983/1‑65 (MQ=255)
           gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAc                                                 >  1:643300/1‑65 (MQ=255)
           gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAc                                                 >  1:1107044/1‑65 (MQ=255)
           gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAc                                                 >  1:1975004/1‑65 (MQ=255)
           gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAc                                                 >  1:1924692/1‑65 (MQ=255)
           gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAc                                                 >  1:3189791/1‑65 (MQ=255)
           gCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAc                                                 >  1:2877814/1‑65 (MQ=255)
            cTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAAcc                                                >  1:2875936/1‑65 (MQ=255)
               ccTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCg                                             <  1:1693897/65‑1 (MQ=255)
                cTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGa                                            >  1:2022171/1‑65 (MQ=255)
                 tGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTCAACCTCGAg                                           >  1:1842660/1‑65 (MQ=255)
                 tGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAg                                           >  1:3808725/1‑65 (MQ=255)
                   agagCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAgtg                                         >  1:1954030/1‑65 (MQ=255)
                     cgCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTc                                       >  1:1842609/2‑65 (MQ=255)
                     agCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTc                                       >  1:874504/1‑65 (MQ=255)
                     agCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTc                                       >  1:3060322/1‑65 (MQ=255)
                     agCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTc                                       >  1:3458874/1‑65 (MQ=255)
                     agCATCAGAACGACCGCCAGAGGCGGGAACCCTCTGTTTTTTTTTCATATATAAAACTCGCgttc                                        >  1:2147598/1‑63 (MQ=255)
                      gCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTcc                                      >  1:1968156/1‑65 (MQ=255)
                         tCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGAt                                   <  1:2601977/65‑1 (MQ=255)
                         tCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGAt                                   <  1:2143939/65‑1 (MQ=255)
                         tCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGAt                                   <  1:1040895/65‑1 (MQ=255)
                         tCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGAt                                   >  1:429448/1‑65 (MQ=255)
                         tCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGAt                                   <  1:907543/65‑1 (MQ=255)
                         tCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGAt                                   <  1:1268830/65‑1 (MQ=255)
                         tCAAAACGGCCCCCAGGGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCGCGAGTGTCCGAt                                   <  1:1506580/65‑1 (MQ=255)
                             aaCGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCa                               >  1:3534181/1‑65 (MQ=255)
                              aCGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTaaa                              >  1:2960067/1‑65 (MQ=255)
                              aCGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCaa                              >  1:3034545/1‑65 (MQ=255)
                              aCGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCaa                              >  1:2999805/1‑65 (MQ=255)
                              aCGACCGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCaa                              >  1:1469541/1‑65 (MQ=255)
                                   cGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtgg                         <  1:3910161/65‑1 (MQ=255)
                                   cGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtgg                         >  1:507807/1‑65 (MQ=255)
                                   cGCCAGAGGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtgg                         <  1:2758178/65‑1 (MQ=255)
                                        gagGCGTAAACCGTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGTTTTGAAATtggtgaat                    >  1:256540/1‑65 (MQ=255)
                                        gagGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaat                    >  1:2390620/1‑65 (MQ=255)
                                        gagGCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaat                    >  1:933400/1‑65 (MQ=255)
                                          ggCGTAACCCCTCTGTTTTTTTTCATATGTAAGCCTCGAGTGTCCGATTTCAAATtggtcaatgg                  <  1:1638066/65‑1 (MQ=255)
                                          ggCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatgg                  >  1:2625667/1‑65 (MQ=255)
                                          ggCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatgg                  >  1:1814128/1‑65 (MQ=255)
                                          ggCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatgg                  >  1:2521635/1‑65 (MQ=255)
                                          ggCGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatgg                  >  1:429310/1‑65 (MQ=255)
                                            cGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatggtc                >  1:15071/1‑65 (MQ=255)
                                            cGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatggtc                >  1:1614017/1‑65 (MQ=255)
                                            cGTAAACCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATtggtcaatggtc                >  1:3353630/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:3453903/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:3435389/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:814836/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:946678/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:3267379/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:1296673/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:1487975/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:1797404/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:1885318/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:2095348/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:248377/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:2614913/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:2743516/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:2982790/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:3328417/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:3142019/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGCCAATGGTCAAAAg           >  1:2615244/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATGTAAAACTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:3903316/1‑65 (MQ=255)
                                                 aCCCTCTGTTTTTTTTCATATATAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAg           >  1:1024592/1‑65 (MQ=255)
                                                  cccTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGt          >  1:1056171/1‑65 (MQ=255)
                                                  cccTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGt          >  1:1182416/1‑65 (MQ=255)
                                                  cccTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGt          >  1:610719/1‑65 (MQ=255)
                                                   ccTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGtt         >  1:2336320/1‑65 (MQ=255)
                                                   ccTCTGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGtt         >  1:1972696/1‑65 (MQ=255)
                                                    ctctGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTta        <  1:873758/65‑1 (MQ=255)
                                                     tctGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTtaa       <  1:1618779/65‑1 (MQ=255)
                                                     tctGTTTTTTTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTtaa       <  1:296369/65‑1 (MQ=255)
                                                           ttttttTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAc  >  1:1491943/1‑65 (MQ=255)
                                                           ttttttTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAc  <  1:2806495/65‑1 (MQ=255)
                                                           ttttttTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAc  <  1:1125461/65‑1 (MQ=255)
                                                           ttttttTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAc  <  1:2463162/65‑1 (MQ=255)
                                                           ttttttTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAc  <  1:3550724/65‑1 (MQ=255)
                                                               |                                                            
GTTAGGGCTAAGCTGCCTGAGAGCATCAGAACGACCGCCAGAGGCGTAAACCCTCTGTTTTTGTTCATATGTAAACCTCGAGTGTCCGATTTCAAATTGGTCAATGGTCAAAAGTTAATAAAC  >  NC_000913/485552‑485674

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