Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 483,728 Δ1 bp coding (676/3150 nt) acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913483,7260T.86.2% 95.3 / 6.4 31coding (678/3150 nt)acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base . (20/5);  minor base C (0/4);  total (21/10)
Fisher's exact test for biased strand distribution p-value = 5.31e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTG  >  NC_000913/483582‑483865
                                                                                                                                                |                                                                                                                                            
gTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACcg                                                                                                                                       >  1:453130/1‑151 (MQ=255)
     caccacCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACcggcgg                                                                                                                                   <  2:104067/150‑1 (MQ=255)
     caccacCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCggcggc                                                                                                                                  >  1:343411/1‑151 (MQ=255)
            cAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTaccac                                                                                                                           <  2:784229/151‑1 (MQ=255)
                 cAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGc                                                                                                                      >  2:603196/1‑151 (MQ=255)
                         gcgACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑Cccgc                                                                                                             >  1:765018/1‑151 (MQ=255)
                         gcgACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTAGAC‑TTCACCGGCGGCGAACCACCGAGCTGA‑Cccgc                                                                                                             >  1:313869/1‑151 (MQ=255)
                          cgACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑Cccgcc                                                                                                            >  1:287391/1‑151 (MQ=255)
                           gACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑Ctcgc                                                                                                             >  1:539458/1‑149 (MQ=255)
                                              cgGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTgcgc                                                                                        >  1:585859/1‑151 (MQ=255)
                                                    ccATCCTGTTTCTCTTGCAGCAGGATTTTGCCAACCTTTTCAGTAGAGGTCGGGGTCGTCGGAGAAAAAATAGGGGCTTTATGGTGTTGG‑CATCCCCCGGCGGCTCCCCCCGGAGCGGC‑CCCGCCGCAACCTGGGCGTTCTGCGCtt                                                                                      <  1:234002/147‑1 (MQ=255)
                                                      aTCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATgg                                                                                >  1:182263/1‑151 (MQ=255)
                                                       tCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGc                                                                               >  1:700536/1‑151 (MQ=255)
                                                        ccTGTTCCATTTACAGCAGTTTTTTGCCAACCTCTTCAGTAGGGGTAGGCGGCGTCTGAAAAAAAATAGGGGCTTAAAGTTGTGGGCC‑TTCCCCGGCGGCGTCCCCCCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCg                                                                              <  2:438177/151‑1 (MQ=255)
                                                               cACTTTCAGCGGGATTTTCCCGATCTTTTCAGTAGAGGTCAGCCGCGTCGGAGAAAAAATGGGGGCTTTAAGCTGTT‑GCCCTCCCCCGGCGGCTTCCCCCCGGGCTGC‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAAt                                                                         <  2:181375/149‑1 (MQ=255)
                                                                aCTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGAc                                                                      >  2:828845/1‑151 (MQ=255)
                                                                     cagcagGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCa                                                                  >  2:594287/1‑150 (MQ=255)
                                                                       ggagatttttGCCGACCTCTTCAGGAGAGGCCAGCCGCGTCTGAAAAAAAAGAGGGGCATTAGGCTGTTGGCC‑TCCCCCGGCGCCGTCCCCCGGGGCTCA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCaa                                                                 <  2:267491/143‑1 (MQ=255)
                                                                              tttGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt                                                         >  2:8617/1‑150 (MQ=255)
                                                                                tGCCAACCTCTGCAGTAGAGGTCAGCCGCGTCTGAAAAATAAGGGGGGCGTAAAGTTG‑TGGCCCTCCCCGGGGGGCCTCCCCCGGAGCTGC‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAg                                                       <  1:109432/150‑1 (MQ=255)
                                                                                      aCTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAAt                                                >  2:461147/1‑151 (MQ=255)
                                                                                        tctTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAAttt                                              >  2:712672/1‑151 (MQ=255)
                                                                                         ctTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTg                                             >  1:631798/1‑151 (MQ=255)
                                                                                                 ggggtcgaccgcgTTGGAGCAAAAATAGGGGCGTTAGGCTGTT‑GCCCTTCCCCGGCGGCGTCCCACCGAGTTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCt                                      <  2:74134/141‑1 (MQ=255)
                                                                                                       aGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCgg                               >  1:420059/1‑151 (MQ=255)
                                                                                                           gcgTCGGAACAAAAATGGGGTCTTAAAGCTGTT‑GCCCTCCCCCGGCGCCTTCCCA‑CGAACCGACCCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCggg                              <  2:264716/148‑1 (MQ=255)
                                                                                                                  aGCAATAATAGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATcc                       >  1:61356/1‑148 (MQ=255)
                                                                                                                        aataGAGGCGTTAAGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGc               >  2:395438/1‑150 (MQ=255)
                                                                                                                        aagagagGCGTTAGGTGGTTGGCC‑TCCCCGGGCGGCGCCCCACCGGGCTGC‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCa              <  1:595804/148‑1 (MQ=255)
                                                                                                                               gCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGt        >  2:443159/1‑151 (MQ=255)
                                                                                                                                    aaGCTGTTGGCC‑TTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACtgtg  >  2:92964/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGA‑CCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTG  >  NC_000913/483582‑483865

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑