Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 470,643 G→A *591* (TAG→TAA
S3N (AGT→AAT) 
mdlA →
mdlB →
putative multidrug ABC transporter ATPase
putative multidrug ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913470,6430GA92.9% 79.6 / ‑0.1 28coding (1773/1773 nt)
coding (8/1782 nt)
mdlA
mdlB
putative multidrug ABC transporter ATPase
putative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/8);  minor base . (1/1);  total (19/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGC  >  NC_000913/470518‑470765
                                                                                                                             |                                                                                                                          
cAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACtg                                                                                                               >  2:702195/1‑139 (MQ=255)
   cGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTgg                                                                                                            >  2:479995/1‑139 (MQ=255)
      ggCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCg                                                                                                         >  1:52984/1‑139 (MQ=255)
        cAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGAc                                                                                                       >  1:234160/1‑139 (MQ=255)
        cAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGAc                                                                                                       <  1:702195/139‑1 (MQ=255)
        cAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTAGCCAACTGTGGCCGAct                                                                                                      >  2:663577/1‑139 (MQ=255)
             atgatgTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCa                                                                                                  >  2:35060/1‑139 (MQ=255)
                   tgCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGcc                                                                                            <  1:664483/139‑1 (MQ=255)
                          acaacaAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGc                                                                                     >  2:217603/1‑139 (MQ=255)
                                  gCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAAtt                                                                                                                          >  2:25544/1‑94 (MQ=255)
                                                   atGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTCAGCGTACGGTTCGCCGTGGCGTAAACCg                                                            <  1:477425/139‑1 (MQ=255)
                                                              tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTg                                                  <  1:253066/138‑1 (MQ=255)
                                                              tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTg                                                  >  2:253066/1‑138 (MQ=255)
                                                              tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                 >  2:78653/1‑139 (MQ=255)
                                                              tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                 >  2:192955/1‑139 (MQ=255)
                                                              tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                 >  1:387953/1‑139 (MQ=255)
                                                                    acTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtg                                          <  1:663577/139‑1 (MQ=255)
                                                                     cTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGc                                                                                                           >  2:227236/1‑74 (MQ=255)
                                                                     cTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGc                                                                                                           <  1:227236/74‑1 (MQ=255)
                                                                        gAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtgatg                                       >  1:33971/1‑139 (MQ=255)
                                                                          ggcggcGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTgatgat                                     <  1:653724/139‑1 (MQ=255)
                                                                          ggcggcGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTgatgat                                     >  2:447649/1‑139 (MQ=255)
                                                                          ggcggcGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTAGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGTGATTGCGGTCCTgatgat                                     >  2:290146/1‑139 (MQ=255)
                                                                             ggcgCTCGACGACGCTCCGGAAAATCGCGAGGAGGCGGTCGATGCGTCATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGAGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGAtgtg                                  >  1:490287/1‑139 (MQ=255)
                                                                                                      cgcgAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGt         >  1:569343/1‑139 (MQ=255)
                                                                                                             aggCCGTCGATGCGTAATTTTAGCCAGCTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGTTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCgc  <  2:490287/139‑1 (MQ=255)
                                                                                                               gCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTggg                                                       <  2:236160/84‑1 (MQ=255)
                                                                                                               gCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTggg                                                       >  1:236160/1‑84 (MQ=255)
                                                                                                                             |                                                                                                                          
CAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGC  >  NC_000913/470518‑470765

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