Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	469,429	A→G	K187E (AAG→GAG)	mdlA →	putative multidrug ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	469,429	0	A	G	100.0%	165.6 / NA	47	K187E (AAG→GAG)	mdlA	putative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/21);  total (26/21)

TGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACC  >  NC_000913/469295‑469560                                                                                                                                       |                                                                                                                                    tGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCt                                                                                                                                 >  1:472021/1‑139 (MQ=255)                             gatgTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAAt                                                                                                     >  1:1124190/1‑139 (MQ=255)                               tgTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCa                                                                                                             >  1:611233/1‑129 (MQ=255)                               tgTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCa                                                                                                             <  2:611233/129‑1 (MQ=255)                                   tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGc                                                                                               >  1:243797/1‑139 (MQ=255)                                   tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGc                                                                                               >  1:447974/1‑139 (MQ=255)                                     cGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcgctctc                                                                                                               <  1:794511/121‑6 (MQ=255)                                     cGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                               >  2:794511/1‑121 (MQ=255)                                          aTTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAggcggc                                                                                                                     <  2:661745/110‑1 (MQ=255)                                          aTTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAggcggc                                                                                                                     >  1:661745/1‑110 (MQ=255)                                               cTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                               <  1:14028/111‑1 (MQ=255)                                               cTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                               >  2:14028/1‑111 (MQ=255)                                                    aGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                              <  2:1109459/139‑1 (MQ=255)                                                    aGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                              <  1:53134/139‑1 (MQ=255)                                                              aTTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                              >  2:355057/1‑129 (MQ=255)                                                              aTTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                              <  1:355057/129‑1 (MQ=255)                                                               ttttCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCa                                                                   <  2:472021/139‑1 (MQ=255)                                                                    ccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATc                                                              >  2:1233808/1‑139 (MQ=255)                                                                     cTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAc                                                                                                  <  1:1145138/102‑1 (MQ=255)                                                                     cTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAc                                                                                                  >  2:1145138/1‑102 (MQ=255)                                                                              aTGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAggcgg                                                                                                                      >  2:123267/1‑73 (MQ=255)                                                                              aTGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAggcgg                                                                                                                      <  1:123267/73‑1 (MQ=255)                                                                                gCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAg                                                                                                                                   <  2:954577/58‑1 (MQ=255)                                                                                gCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAg                                                                                                                                   >  1:954577/1‑58 (MQ=255)                                                                                     tGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                              <  1:132875/106‑1 (MQ=255)                                                                                     tGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                              >  2:132875/1‑106 (MQ=255)                                                                                                  gATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCg                                <  1:1233808/139‑1 (MQ=255)                                                                                                    ttAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGc                              >  2:231066/1‑139 (MQ=255)                                                                                                    ttAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGc                              >  1:214712/1‑139 (MQ=255)                                                                                                       aGCGCAACGGCGATGCTTTACATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCtt                                                      >  1:220363/1‑112 (MQ=255)                                                                                                       aGCGCAACGGCGATGCTTTACATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCtt                                                      <  2:220363/112‑1 (MQ=255)                                                                                                         cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTa                                                                         <  2:724960/91‑1 (MQ=255)                                                                                                         cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTa                                                                         >  1:724960/1‑91 (MQ=255)                                                                                                            aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGc                                                                    >  2:359466/1‑93 (MQ=255)                                                                                                            aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGc                                                                    <  1:359466/93‑1 (MQ=255)                                                                                                            aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCtt                                                      <  1:11490/107‑1 (MQ=255)                                                                                                            aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCtt                                                      >  2:11490/1‑107 (MQ=255)                                                                                                              cGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCAc                                                                                             <  2:637177/66‑1 (MQ=255)                                                                                                              cGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCAc                                                                                             >  1:637177/1‑66 (MQ=255)                                                                                                                  gATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCgg                <  2:447974/139‑1 (MQ=255)                                                                                                                   aTGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCa                                                                                              >  2:532503/1‑60 (MQ=255)                                                                                                                   aTGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCa                                                                                              <  1:532503/60‑1 (MQ=255)                                                                                                                      cTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                              >  1:1370799/1‑73 (MQ=255)                                                                                                                      cTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                              <  2:1370799/73‑1 (MQ=255)                                                                                                                         tGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGa         >  2:522227/1‑139 (MQ=255)                                                                                                                               aaCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATAc   >  2:1086569/1‑139 (MQ=255)                                                                                                                                aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATAct  >  1:1118591/1‑138 (MQ=255)                                                                                                                                       |                                                                                                                                    TGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACC  >  NC_000913/469295‑469560

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