Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	470,620	C→A	R584S (CGC→AGC)	mdlA →	ABC transporter family protein MdlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	470,620	0	C	A	100.0%	113.2 / NA	33	R584S (CGC→AGC)	mdlA	ABC transporter family protein MdlA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (20/13);  total (20/13)

GTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAA  >  NC_000913/470494‑470752                                                                                                                               |                                                                                                                                     gtgATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGt                                                                                                                         >  2:557820/1‑140 (MQ=255)        agcaCGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCg                                                                                                                  <  1:245739/140‑1 (MQ=255)          caCGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGTCGTCGATGCGTa                                                                                                                >  1:637280/1‑140 (MQ=255)                    tCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCa                                                                                                      >  1:547943/1‑140 (MQ=255)                     cGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCaa                                                                                                     <  1:557820/140‑1 (MQ=255)                           gCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTgg                                                                                               >  1:113780/1‑140 (MQ=255)                            cGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGc                                                                                              >  1:285065/1‑140 (MQ=255)                                         atgTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAgcg                                                                                 >  1:648187/1‑140 (MQ=255)                                              cTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGt                                                                            >  1:108774/1‑140 (MQ=255)                                                ggCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTa                                                                          >  1:198597/1‑140 (MQ=255)                                                 gCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAg                                                                         >  2:247340/1‑140 (MQ=255)                                                   acaacaAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCg                                                                       >  1:281500/1‑140 (MQ=255)                                                      acaaAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTAc                                                                    >  1:27181/1‑140 (MQ=255)                                                           gCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTc                                                               >  2:648556/1‑140 (MQ=255)                                                           gCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTc                                                               >  1:217576/1‑140 (MQ=255)                                                              gCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGcc                                                            <  2:27181/140‑1 (MQ=255)                                                                      cgcgATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGTCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGTCTGTTAGCGTACGGTTCGCCGTGGCGTa                                                    <  2:637280/140‑1 (MQ=255)                                                                        cgATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTaaa                                                  <  1:475408/140‑1 (MQ=255)                                                                                     tATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTg                                     <  2:316928/140‑1 (MQ=255)                                                                                      aTCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                    <  2:108774/140‑1 (MQ=255)                                                                                       tCAACAACTGGAGGCGGCGCTTGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCg                                   >  1:314209/1‑140 (MQ=255)                                                                                       tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCg                                   >  2:27812/1‑140 (MQ=255)                                                                                         aacaacTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGTGGt                                 >  1:296291/1‑140 (MQ=255)                                                                                         aacaacTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGt                                 <  2:463714/140‑1 (MQ=255)                                                                                             acTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtg                             <  2:285065/140‑1 (MQ=255)                                                                                             acTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtg                             <  2:547943/140‑1 (MQ=255)                                                                                             acTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtg                             <  1:27812/140‑1 (MQ=255)                                                                                              cTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTg                                     <  2:505728/131‑1 (MQ=255)                                                                                              cTGGAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTg                                     >  1:505728/1‑131 (MQ=255)                                                                                                 gAGGCGGCGCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtgatga                         >  1:637554/1‑140 (MQ=255)                                                                                                     cggcgCTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGAtgtg                     >  1:369721/1‑140 (MQ=255)                                                                                                        cgcTCGACGACGCTCCGGAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGt                  >  1:635139/1‑140 (MQ=255)                                                                                                                        ggAAAATAGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGaa  <  2:281500/140‑1 (MQ=255)                                                                                                                               |                                                                                                                                     GTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAA  >  NC_000913/470494‑470752

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑