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% 82.0 / NA 24A1002T (GCA→ACA) acrBmultidrug efflux pump RND permease AcrB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/8);  total (16/8)

AACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAA  >  NC_000913/481338‑481460
                                                              |                                                            
aaCGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTgtgt                                                            <  1:1526434/65‑1 (MQ=255)
          aTTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGCTCGGCGCGcc                                                  >  1:515251/1‑65 (MQ=25)
             gCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGa                                               >  1:3521828/1‑65 (MQ=255)
              ccAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGaa                                              >  1:1810165/1‑65 (MQ=255)
               cAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAAc                                             <  1:2144374/65‑1 (MQ=255)
                   aCCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGc                                         >  1:1890394/1‑65 (MQ=255)
                    ccGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCa                                        >  1:623781/1‑65 (MQ=255)
                    ccGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCa                                        <  1:580903/65‑1 (MQ=255)
                    ccGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCa                                        <  1:1669910/65‑1 (MQ=255)
                      gTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCAcc                                      <  1:1427135/65‑1 (MQ=255)
                         gCGGTCACCACCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGt                                   <  1:2585480/65‑1 (MQ=37)
                            gTCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACt                                <  1:535368/65‑1 (MQ=255)
                             tCACCATCCCGCCCATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTg                               >  1:3266016/1‑65 (MQ=255)
                                         ccATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGc                   <  1:733260/65‑1 (MQ=255)
                                         ccATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGc                   >  1:1750708/1‑65 (MQ=255)
                                         ccATTACACCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGc                   >  1:849817/1‑65 (MQ=255)
                                         ccATTACACCGGTACCTACTGTGTGCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGc                   <  1:2754518/65‑1 (MQ=255)
                                                acCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGc            >  1:2952103/1‑65 (MQ=255)
                                                acCGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGc            <  1:2194468/65‑1 (MQ=255)
                                                 ccGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCGCCAGTACTGATAACCAGCGGCATAACGcc           >  1:918331/1‑65 (MQ=25)
                                                 ccGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGcc           >  1:3172239/1‑65 (MQ=255)
                                                 ccGGTACCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGcc           >  1:3023718/1‑65 (MQ=255)
                                                      aCCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGa      >  1:170989/1‑65 (MQ=255)
                                                      aCCTACTGTGTTCTGCGCGCCGGAACCAGCACCAGCACTGATAACCAGCGGCATAACGCCGAGGa      >  1:2856752/1‑65 (MQ=25)
                                                        cTACTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACGAGCGGCATAACGCCGAGGATa    >  1:2037723/1‑65 (MQ=255)
                                                          aCTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATaaa  >  1:342765/1‑65 (MQ=255)
                                                          aCTGTGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATaaa  >  1:308873/1‑65 (MQ=255)
                                                              |                                                            
AACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAA  >  NC_000913/481338‑481460

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