Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,581,278	A→G	V62A (GTT→GCT)	hsdS ←	type I restriction enzyme EcoKI specificity subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,581,278	0	A	G	100.0%	38.4 / NA	13	V62A (GTT→GCT)	hsdS	type I restriction enzyme EcoKI specificity subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/5);  total (8/5)

AGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATC  >  NC_000913/4581116‑4581509                                                                                                                                                                   |                                                                                                                                                                                                                                        aGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTTCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAGCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTAt                                                                                                                                                 <  2:245975/251‑1 (MQ=255)         aGGACGTAATAAACCGAAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGCCATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAGCAACAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGct                                                                                                                                         >  2:169151/1‑251 (MQ=255)                          aaaaTGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTCCAAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGAAATTTTTTGACTTTCTTTAACAAGATTTTTAGGAGCAAAAACCAAGTCCGTAGTATCAAACTAGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCACTTTTTTATACGTTAct                                                                                                                        >  1:174259/1‑251 (MQ=255)                                                             tGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAGCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGc                                                                                     >  2:38426/1‑251 (MQ=255)                                                                             cAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTATGAGCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCCCGGATAAGAGTCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTAGACGTTACTCCTCGGATTAGAGTTGTGACCGTATATACTGGTGCGATACCCCACCCCTcc                                                                     <  2:174258/251‑1 (MQ=255)                                                                              aaCTACGGATTTGCTGCCTGATGACAATGCAATAACAATATCTTCAGGATATAGTTTTTGACTTTCTTTAACAAGATTTGTAGGAGCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCTCGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTTGTTTTTATACCTTACTGCTGGGAATAGATTTGGGACCGTGGATACTGGGGCGATCCCCCAACCCTCCg                                                                    <  2:283283/251‑1 (MQ=255)                                                                              aaCTACGGATTTGCTCCCTGATGACATTGCCATAAAAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAGCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCg                                                                    >  1:107323/1‑251 (MQ=255)                                                                                      aTTTGCTCCCTGATGACACTGTAATAACAATATCTTCAGTAGATATTTTTTGACTTTTTTAAACAAGATTTTTAGGAGCACAAACCAAATCCATAGTATCAAACATGCCATTCTGCATAATGTTCACCCGGATAACAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATAGGTTACTCCTCGGATTATAGTTGTGACCGTATTTACTAGGGCGATACCCCACCCCTCCGGCAATTTc                                                            >  2:206604/1‑251 (MQ=255)                                                                                           cTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAGCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGc                                                       >  1:133677/1‑251 (MQ=255)                                                                                                                       cTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAGCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGctt                           >  2:203317/1‑251 (MQ=255)                                                                                                                         tcatgagaTATTTTTTGACATTCTTTAACAAGATTTTTAGGAGCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCCCCAAACGc                             <  2:203720/243‑1 (MQ=255)                                                                                                                                      tttGACTTTCTTTAACAAGATTTTTAGGAGCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGa             <  1:156835/250‑1 (MQ=255)                                                                                                                                                tttAACAAGATTTTTAGGAGCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATc  >  1:162435/1‑251 (MQ=255)                                                                                                                                                                   |                                                                                                                                                                                                                                        AGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATC  >  NC_000913/4581116‑4581509

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