Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,582,590	A→G	R153R (CGT→CGC)	hsdM ←	type I restriction enzyme EcoKI methyltransferase component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,582,590	0	A	G	100.0%	73.5 / NA	21	R153R (CGT→CGC)	hsdM	type I restriction enzyme EcoKI methyltransferase component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/9);  total (12/9)

GTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATC  >  NC_000913/4582455‑4582718                                                                                                                                        |                                                                                                                                  gtcgtcCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGTTTCAATAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGcgc‑g                                                                                                                              <  2:287185/140‑1 (MQ=255) gtcgtcCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGcgc‑g                                                                                                                              >  1:24843/1‑140 (MQ=255)      ccAGATCTTTGGTTTGCGACTAAACATGGCGGTCGGCTTAAATAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTg                                                                                                                         <  2:13695/140‑1 (MQ=255)               tGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGcgc‑g                                                                                                                              <  1:17415/126‑1 (MQ=255)               tGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGcgc‑g                                                                                                                              >  2:17415/1‑126 (MQ=255)                      cGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCAc                                                                                                         >  1:221029/1‑140 (MQ=255)                                  gCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGttt                                                                                             <  2:257245/140‑1 (MQ=255)                                   cGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTc                                                                                            >  2:349481/1‑140 (MQ=255)                                           ttCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGt                                                                                    >  1:27556/1‑140 (MQ=255)                                                 aaaaaaGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGCGGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGttcttc                                                                               <  2:24843/139‑1 (MQ=255)                                                          cgccgTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGcc                                                                     >  2:89416/1‑140 (MQ=255)                                                                               cTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGt                                                >  1:87051/1‑140 (MQ=255)                                                                                     cacTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATcgc                                          >  2:303064/1‑140 (MQ=255)                                                                                       cTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATcgcgc                                        >  1:4325/1‑140 (MQ=255)                                                                                        ttCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATcgcgcg                                       >  1:22599/1‑140 (MQ=255)                                                                                            cgcgGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACtt                                   >  1:239232/1‑140 (MQ=255)                                                                                                     ggTTTCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgcgc                          <  1:97686/140‑1 (MQ=255)                                                                                                         tCAGCAGATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCg                      <  1:303064/140‑1 (MQ=255)                                                                                                          cagcagATGAATAATGGTTTTAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGt                     <  2:22599/140‑1 (MQ=255)                                                                                                                           ttttAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGaa    <  2:87051/140‑1 (MQ=255)                                                                                                                             ttAATCAGCGGGCGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATc  >  1:138988/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                  GTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGC‑GGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATC  >  NC_000913/4582455‑4582718

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