Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,246,296	C→T	G42S (GGC→AGC)	malE ←	maltose ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,246,296	0	C	T	100.0%	72.5 / NA	22	G42S (GGC→AGC)	malE	maltose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/10);  total (12/10)

CGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCCTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATAATGCGAGGATGCGTGCACCTGTTTTTATTTTCATAATCTATGGTCCTTGTTGGTGAAGTGCTCGTGAAAACACCTAAACGGACTCTAGTTTCTTTATACGGCAACCTCTTTCCATCCTCCTTGCCCCTACGCCCCACCGTCGCTTTGTGTGATCT  >  NC_000913/4246049‑4246540                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     cGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCTGGGTGATTTCAGCCAACAGGCCAGATGGAGCGTAGCCACCAAAGCGGTCGTGTGCCCAGAATATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGTAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCtttt                                                                                                                                                                                                                                                   <  2:227317/251‑1 (MQ=255)      ggCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGc                                                                                                                                                                                                                                              <  1:293161/251‑1 (MQ=255)       gCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGcc                                                                                                                                                                                                                                             <  1:214519/251‑1 (MQ=255)                      ggATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTAc                                                                                                                                                                                                                              >  2:219417/1‑251 (MQ=255)                          aCAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAg                                                                                                                                                                                                                           <  1:287197/250‑1 (MQ=255)                                         aCGCTTTGTCCGGGGAGCTTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGCCGTGTGCTCAGAAGATAATGACAGGGCCATAGCCCGTTGCCGCAACCTGTGGGAATTACTCTTCCAGTTTATCCGGAAGCTCAACGGTGCCTTTAATTCTGATATCTTTCTCGAATTTCTTAACGACTTCAGCGAGACCGCTATAACTTTTATCGCCGTTAATCCAGATTACCAGTTTACCCTCTTCGAtt                                                                                                                                                                                                           >  1:408329/1‑251 (MQ=255)                                                  ccGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCg                                                                                                                                                                                                     >  1:233150/1‑248 (MQ=255)                                                             tCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCgg                                                                                                                                                                                        >  1:52902/1‑250 (MQ=255)                                                                gCCAACAGGACAGATTGGGCGTAGCCACCAAAGCGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCTCACCCTGGTGGAATTTCTCTTCCAGTTTATCGTGATGCTCATCGGTGACTTTAATTCCTGTATCTTTCTCGAATTTCTTACCGACTTCAGCGATACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGaaaa                                                                                                                                                                                    <  2:309697/251‑1 (MQ=255)                                                                              ttGAGCGTAGCCACCAAAGCGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATcgtcgtg                                                                                                                                                                       >  2:329511/1‑249 (MQ=255)                                                                                       gCCACCAAAGCGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATa                                                                                                                                                             <  1:97719/251‑1 (MQ=255)                                                                                                 cGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATAATGCGAGGAt                                                                                                                                                   <  2:102938/251‑1 (MQ=255)                                                                                                          gCCCAGAAGATAATGTCAGGTCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATAATGCGAGGATGCGTGCAc                                                                                                                                           <  1:144269/250‑1 (MQ=255)                                                                                                                        gTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATAATGCGAGGATGCGTGCACCTGTTTTTATTTTc                                                                                                                             >  2:73515/1‑250 (MQ=255)                                                                                                                                   cGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATTCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATAATGCGAGGATGCGTGCACCTGTTTTTATTATCATAATCTATgg                                                                                                                  <  1:210403/250‑1 (MQ=255)                                                                                                                                                   cTGTGGGAATTTCTCTTCCAGTTTATCCGGATTCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTAGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATAATGCGAGGATGCGTGCACCTGTTTTTATTTTCATAATCTATGGTCCTTGTTGGTGAAGTg                                                                                                 >  2:144483/1‑251 (MQ=255)                                                                                                                                                            tttCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTATTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATAATGCGAGGATGCGTGCACCTGTTTTTATTTTCATAATCTATGGTCCTTGTTGGTGAAGTGCTCGTGaa                                                                                         >  2:239651/1‑250 (MQ=255)                                                                                                                                                                                       aaCGGCGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGTGGCGCAAAACATCATCGTCGTTAATGCGGATAATGCGAGGATGCGTGCACCTGTTTTTATTTTCATACTCTATGCTCCTTGTTGTTGAAGTGCTCGTGAAAACACCTAAACGGACTCTAGTTTCttt                                                              >  2:52857/1‑250 (MQ=255)                                                                                                                                                                                           gTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATAATGCGAGGATGCGTGCACCTGTTTTTATTTTCATAATCTATGGTCCTTGTTGGTGAAGTGCTCGTGAAAACACCTAAACGGACTCTAGTTTCTTTATACg                                                         <  1:265691/251‑1 (MQ=255)                                                                                                                                                                                                                  ctcGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATAATGCGAGGATGCGTGCACCTGTTTTTATTTTCATAATCTATGGTCCTTGTTGGTGAAGTGCTCGTGAAAACACCTAAACGGACTCTAGTTTCTTTATACGGCAACCTCTTTCCATCCTCCTTg                                  >  2:125414/1‑251 (MQ=255)                                                                                                                                                                                                                                cGACTTCAGCGAGACCGTTATAGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATAATGCGAGGATGCGTGCACCTGTTTTTATTTTCATAATCTATGGTCCTTGTTGGTGAAGTGCTCGTGAAAACACCTAAACGGACTCTAGTTTCTTTATACGGCAACCTCTTTCCATCCTCCTTGCCCCTACGCCCCAc                    >  1:82283/1‑251 (MQ=255)                                                                                                                                                                                                                                                gTTATAGCTTTTATCGCCGTTAATCGATATTACCAGTTTACCGTCTTCGATTTGGGCGAGAGCCGAGGCTGAAAACATCATCGTCGTTAAGGCGGATAATGCGAGGATGCCTGCACCTTTTTTTAATTTCATAAACGATGGTCCTTGTTGGCGAAGTGCACGTTAAAACACCTAAACGGACTCTAGTTTCTTTATTAGTCAACCTCGCTCCTTCCTCCTTGCCCCCACGCTCCACCTTCGcgaggtctgat    <  2:405615/251‑11 (MQ=255)                                                                                                                                                                                                                                                  tataGCTTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATAATGCGAGGATGCGTGCACCTGTTTTTATTTTCATAATCTATGGTCCTTGTTGGTGAAGTGCTCGTGAAAACACCTAAACGGACTCTAGTTTCTTTATACGGCAACCTCTTTCCATCCTCCTTGCCCCTACGCCCCACCGTCGCTTTGTGTGAtct  >  2:228165/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     CGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTCTCGAATTTCTTACCGACTTCAGCGAGACCGTTATAGCCTTTATCGCCGTTAATCCAGATTACCAGTTTACCTTCTTCGATTTTGGCGAGAGCCGAGGCGGAAAACATCATCGTCGTTAATGCGGATAATGCGAGGATGCGTGCACCTGTTTTTATTTTCATAATCTATGGTCCTTGTTGGTGAAGTGCTCGTGAAAACACCTAAACGGACTCTAGTTTCTTTATACGGCAACCTCTTTCCATCCTCCTTGCCCCTACGCCCCACCGTCGCTTTGTGTGATCT  >  NC_000913/4246049‑4246540

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