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

GTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTG  >  NC_000913/483582‑483865
                                                                                                                                                |                                                                                                                                            
gTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACcg                                                                                                                                       >  1:453130/1‑151 (MQ=255)
     caccacCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACcggcgg                                                                                                                                   <  2:104067/150‑1 (MQ=255)
     caccacCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCggcggc                                                                                                                                  >  1:343411/1‑151 (MQ=255)
            cAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTaccac                                                                                                                           <  2:784229/151‑1 (MQ=255)
                 cAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGc                                                                                                                      >  2:603196/1‑151 (MQ=255)
                         gcgACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACccgc                                                                                                             >  1:765018/1‑151 (MQ=255)
                         gcgACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTAGACTTCACCGGCGGCGAACCACCGAGCTGACccgc                                                                                                             >  1:313869/1‑151 (MQ=255)
                          cgACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACccgcc                                                                                                            >  1:287391/1‑151 (MQ=255)
                           gACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACtcgc                                                                                                             >  1:539458/1‑149 (MQ=255)
                                              cgGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTgcgc                                                                                        >  1:585859/1‑151 (MQ=255)
                                                    ccATCCTGTTTCTCTTGCAGCAGGATTTTGCCAACCTTTTCAGTAGAGGTCGGGGTCGTCGGAGAAAAAATAGGGGCTTTATGGTGTTGGCATCCCCCGGCGGCTCCCCCCGGAGCGGCCCCGCCGCAACCTGGGCGTTCTGCGCtt                                                                                      <  1:234002/147‑1 (MQ=255)
                                                      aTCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATgg                                                                                >  1:182263/1‑151 (MQ=255)
                                                       tCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGc                                                                               >  1:700536/1‑151 (MQ=255)
                                                        ccTGTTCCATTTACAGCAGTTTTTTGCCAACCTCTTCAGTAGGGGTAGGCGGCGTCTGAAAAAAAATAGGGGCTTAAAGTTGTGGGCCTTCCCCGGCGGCGTCCCCCCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCg                                                                              <  2:438177/151‑1 (MQ=255)
                                                               cACTTTCAGCGGGATTTTCCCGATCTTTTCAGTAGAGGTCAGCCGCGTCGGAGAAAAAATGGGGGCTTTAAGCTGTTGCCCTCCCCCGGCGGCTTCCCCCCGGGCTGCCCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAAt                                                                         <  2:181375/149‑1 (MQ=255)
                                                                aCTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGAc                                                                      >  2:828845/1‑151 (MQ=255)
                                                                     cagcagGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCa                                                                  >  2:594287/1‑150 (MQ=255)
                                                                       ggagatttttGCCGACCTCTTCAGGAGAGGCCAGCCGCGTCTGAAAAAAAAGAGGGGCATTAGGCTGTTGGCCTCCCCCGGCGCCGTCCCCCGGGGCTCACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCaa                                                                 <  2:267491/143‑1 (MQ=255)
                                                                              tttGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt                                                         >  2:8617/1‑150 (MQ=255)
                                                                                tGCCAACCTCTGCAGTAGAGGTCAGCCGCGTCTGAAAAATAAGGGGGGCGTAAAGTTGTGGCCCTCCCCGGGGGGCCTCCCCCGGAGCTGCCCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAg                                                       <  1:109432/150‑1 (MQ=255)
                                                                                      aCTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAAt                                                >  2:461147/1‑151 (MQ=255)
                                                                                        tctTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAAttt                                              >  2:712672/1‑151 (MQ=255)
                                                                                         ctTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTg                                             >  1:631798/1‑151 (MQ=255)
                                                                                                 ggggtcgaccgcgTTGGAGCAAAAATAGGGGCGTTAGGCTGTTGCCCTTCCCCGGCGGCGTCCCACCGAGTTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCt                                      <  2:74134/141‑1 (MQ=255)
                                                                                                       aGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCgg                               >  1:420059/1‑151 (MQ=255)
                                                                                                           gcgTCGGAACAAAAATGGGGTCTTAAAGCTGTTGCCCTCCCCCGGCGCCTTCCCACGAACCGACCCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCggg                              <  2:264716/148‑1 (MQ=255)
                                                                                                                  aGCAATAATAGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATcc                       >  1:61356/1‑148 (MQ=255)
                                                                                                                        aataGAGGCGTTAAGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGc               >  2:395438/1‑150 (MQ=255)
                                                                                                                        aagagagGCGTTAGGTGGTTGGCCTCCCCGGGCGGCGCCCCACCGGGCTGCCCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCa              <  1:595804/148‑1 (MQ=255)
                                                                                                                               gCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGt        >  2:443159/1‑151 (MQ=255)
                                                                                                                                    aaGCTGTTGGCCTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACtgtg  >  2:92964/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GTTCTCACCACCCAGCTCAATCTTCGCGACGTCACGCAGCAGCACGCGGGAACCATCCTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTG  >  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: