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%	156.0 / NA	43	K187E (AAG→GAG)	mdlA	putative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/21);  total (22/21)

TGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGAT  >  NC_000913/469298‑469548                                                                                                                                    |                                                                                                                        tGATGGGCTGCGCTGAGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGAGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGc                                                                                                                  >  1:674934/1‑139 (MQ=255)   aTGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGcac                                                                                                                >  2:115055/1‑139 (MQ=255)        cTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAggcg                                                                                                           >  1:606175/1‑139 (MQ=255)          gcgcTGTGTTGATTATGATGTGTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTGAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAggcggc                                                                                                         >  1:608291/1‑139 (MQ=255)                tgtTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                   <  2:495030/139‑1 (MQ=255)                tgtTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                   >  2:236568/1‑139 (MQ=255)                          gatgTCTACGCAAATTGGCTGGAAGTTGACCTTATTTTCCCTGTTGCCTTTGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGGACAGGCGGCGGTTTCCAGTCTTAAt                                                                                         >  1:842354/1‑139 (MQ=255)                          gatgTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                   <  1:408616/129‑1 (MQ=255)                          gatgTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                   >  2:408616/1‑129 (MQ=255)                          gatgTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAAt                                                                                         >  2:537639/1‑139 (MQ=255)                                tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcg                                                                                                        <  2:857928/118‑1 (MQ=255)                                tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcg                                                                                                        >  1:857928/1‑118 (MQ=255)                                tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAcc                                                                                     <  1:943817/137‑1 (MQ=255)                                tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAcc                                                                                     >  2:943817/1‑137 (MQ=255)                                tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGc                                                                                   >  2:409447/1‑139 (MQ=255)                                       aTTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAgg                                                                            >  1:793655/1‑139 (MQ=255)                                           gCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAg                                                                        <  1:115055/139‑1 (MQ=255)                                                cAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCa                                                                   <  1:537639/139‑1 (MQ=255)                                                 aGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGGACGCTTAGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                  <  1:556131/139‑1 (MQ=255)                                                 aGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                  >  2:696690/1‑139 (MQ=255)                                                      aCCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTa                                                             >  1:936913/1‑139 (MQ=255)                                                            ttttCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCa                                                       <  1:696690/139‑1 (MQ=255)                                                                cccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                   <  1:782340/91‑1 (MQ=255)                                                                cccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                   >  2:782340/1‑91 (MQ=255)                                                                cccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGAt                                                   >  2:250234/1‑139 (MQ=255)                                                                 ccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                   >  1:602485/1‑90 (MQ=255)                                                                 ccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                   <  2:602485/90‑1 (MQ=255)                                                                 ccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGATTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATc                                                  <  2:608291/139‑1 (MQ=255)                                                                  cTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTCTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCCGGAAAGCCTCACCAGTATCCGCATGATCa                                                 <  2:674934/139‑1 (MQ=255)                                                                         cGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCtt                                          <  2:793655/139‑1 (MQ=255)                                                                         cGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCtt                                          <  1:236568/139‑1 (MQ=255)                                                                         cGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCtt                                          <  1:409447/139‑1 (MQ=255)                                                                           aTGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTg                                        <  1:250234/139‑1 (MQ=255)                                                                                 gTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTgg                                  <  2:606175/139‑1 (MQ=255)                                                                                                      cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                  <  2:490887/86‑1 (MQ=255)                                                                                                      cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                  >  1:490887/1‑86 (MQ=255)                                                                                                         aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                  >  2:803350/1‑83 (MQ=255)                                                                                                         aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                  <  1:803350/83‑1 (MQ=255)                                                                                                         aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTg          >  2:816469/1‑139 (MQ=255)                                                                                                              cGATGCTTTGCATGAACACTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGa                                                                                       <  2:500387/57‑1 (MQ=255)                                                                                                              cGATGCTTTGCATGAACACTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGa                                                                                       >  1:500387/1‑57 (MQ=255)                                                                                                                             aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGAt  <  1:476662/127‑1 (MQ=255)                                                                                                                             aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGAt  >  2:476662/1‑127 (MQ=255)                                                                                                                                    |                                                                                                                        TGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGAT  >  NC_000913/469298‑469548

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