Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 481,400 C→T A1002T (GCA→ACA)  acrB ← multidrug efflux pump RND permease AcrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913481,4000CT100.0% 112.3 / NA 33A1002T (GCA→ACA) acrBmultidrug efflux pump RND permease AcrB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/11);  total (22/11)

AACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAAC  >  NC_000913/481338‑481461
                                                              |                                                             
aaCGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTgtgt                                                             <  1:3682227/65‑1 (MQ=255)
  cGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGGGTTc                                                           >  1:3826686/1‑65 (MQ=255)
      gaagaTTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTgcg                                                       >  1:3552522/1‑65 (MQ=255)
          aTTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGcc                                                   <  1:3092953/65‑1 (MQ=255)
             gCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGa                                                >  1:2824876/1‑65 (MQ=255)
             gCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGa                                                >  1:1582326/1‑65 (MQ=255)
             gCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGa                                                >  1:3578124/1‑65 (MQ=255)
                 gTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCa                                            >  1:2600485/1‑65 (MQ=255)
                   aCCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGc                                          >  1:2887995/1‑65 (MQ=255)
                   aCCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGc                                          >  1:3085799/1‑65 (MQ=255)
                    ccGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCa                                         <  1:2900995/65‑1 (MQ=255)
                    ccGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCa                                         <  1:2943605/65‑1 (MQ=255)
                    ccGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCa                                         >  1:1328000/1‑65 (MQ=255)
                     cGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGccc                                        >  1:2931473/1‑63 (MQ=255)
                      gTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCAcc                                       <  1:1024046/65‑1 (MQ=255)
                            gTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACt                                 <  1:1194697/65‑1 (MQ=255)
                               aCCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGAt                              >  1:84519/1‑65 (MQ=255)
                                  aTCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAAc                           >  1:2910116/1‑65 (MQ=255)
                                   tCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAAcc                          <  1:84848/65‑1 (MQ=255)
                                       gCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCg                      >  1:3499856/1‑65 (MQ=255)
                                         ccATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGc                    <  1:1084628/65‑1 (MQ=255)
                                         ccATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGc                    <  1:1186481/65‑1 (MQ=255)
                                         ccATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGc                    <  1:495138/65‑1 (MQ=255)
                                         ccATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGc                    <  1:1766288/65‑1 (MQ=255)
                                         ccATCACCCCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGc                    <  1:1718814/65‑1 (MQ=18)
                                               cacCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACg              >  1:2615030/1‑65 (MQ=255)
                                               cacCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACg              >  1:3257832/1‑65 (MQ=255)
                                                 ccGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGcc            >  1:553164/1‑65 (MQ=255)
                                                 ccGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGcc            >  1:739310/1‑65 (MQ=255)
                                                   ggTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAGCGCAGa          >  1:3899268/1‑65 (MQ=21)
                                                   ggTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGa          >  1:2395501/1‑65 (MQ=255)
                                                     tACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAgg        <  1:3612426/65‑1 (MQ=255)
                                                     tACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAgg        <  1:3056127/65‑1 (MQ=255)
                                                      aCCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGa       >  1:1810091/1‑65 (MQ=255)
                                                      aCCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGa       >  1:489589/1‑65 (MQ=255)
                                                       ccTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGAt      <  1:3967239/65‑1 (MQ=255)
                                                        cTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATa     >  1:61875/1‑65 (MQ=255)
                                                        cTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATa     >  1:2670375/1‑65 (MQ=255)
                                                          aCTGTGTTCTGCGCGCCGGAACCAGCACCCGTACTGATAACCAGCGGCATAACGCCGAGGATaaa   >  1:2916763/1‑65 (MQ=255)
                                                          aCTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATaaa   >  1:3643364/1‑65 (MQ=255)
                                                          aCTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATaaa   >  1:795322/1‑65 (MQ=255)
                                                          aCTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATaaa   >  1:2843033/1‑65 (MQ=255)
                                                          aCTGTGTTCTGCGCGCCGGAACCAGCACCAGAACTGATAACCAGCGGCATAGCGCGGAGGACaaa   >  1:288507/1‑65 (MQ=21)
                                                           cTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAAc  <  1:2626644/65‑1 (MQ=255)
                                                              |                                                             
AACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAAC  >  NC_000913/481338‑481461

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