Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 482,533 G→A T624I (ACC→ATC)  acrB ← multidrug efflux pump RND permease AcrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913482,5330GA87.5% 16.7 / ‑3.4 8T624I (ACC→ATC) acrBmultidrug efflux pump RND permease AcrB
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (1/6);  total (2/6)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

CAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGCGAACACCGACTCAACGTTGT  >  NC_000913/482471‑482593
                                                              |                                                            
cAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGATATTc                                                         <  1:1991722/68‑1 (MQ=255)
  aCTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGATATTCTg                                                       >  1:992085/1‑68 (MQ=255)
         tttCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGATATTCTGACCACGt                                                <  1:2569333/68‑1 (MQ=255)
                  ccGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGATATTCTGACCACGTCCCGCAaag                                       <  1:1326447/68‑1 (MQ=255)
                   cGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGATATTCTGACCACGTCCCGCAaagc                                      <  1:2574416/68‑1 (MQ=255)
                   cGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGATATTCTGACCACGTCCCGCAaagc                                      <  1:2639424/68‑1 (MQ=255)
                   cGGACGATCGGCCCAGGCCTTCAAGGAAACGAACGCAATACCGATATTCTGACCACGTCCCGCAaagc                                      <  1:2144742/68‑1 (MQ=255)
                                   tCCTTCAAGAAAACGAACGCAATACCGATATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGc                      <  1:1780286/68‑1 (MQ=255)
                                    tcTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGCg                     >  1:1375080/2‑68 (MQ=255)
                                                       aaTGCCGATATGCCGCCCCCGCCCCGCAAAGCCGAAGCCGTTAACGGCGAACAGCGACTCAACgttgt  <  1:2922102/68‑1 (MQ=21)
                                                              |                                                            
CAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGCGAACACCGACTCAACGTTGT  >  NC_000913/482471‑482593

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: