Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,583,657	G→A	L1035L (CTC→CTT)	hsdR ←	type I restriction enzyme EcoKI endonuclease component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,583,657	0	G	A	100.0%	111.0 / NA	31	L1035L (CTC→CTT)	hsdR	type I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/15);  total (16/15)

CATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGT  >  NC_000913/4583580‑4583727                                                                              |                                                                       cATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTaaga                                                                     <  2:951054/81‑1 (MQ=255)   tGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATc                                                                   <  2:951055/81‑1 (MQ=255)        tGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATcgcg                                                                <  1:951057/79‑1 (MQ=255)        tGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATcgcg                                                                <  1:951056/79‑1 (MQ=255)          cgcAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCggg                                                            <  1:951058/81‑1 (MQ=255)           gcAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGc                                                           >  1:951088/1‑81 (MQ=255)               ggAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGAtt                                                       <  2:951059/81‑1 (MQ=255)                 aaGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATtaa                                                     >  1:951114/1‑81 (MQ=255)                    aTTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATtaataa                                                  >  2:951095/1‑81 (MQ=255)                    aTTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATtaataa                                                  >  1:951098/1‑81 (MQ=255)                     ttCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAAc                                                 <  1:951060/81‑1 (MQ=255)                       ccTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAActg                                               >  2:951100/1‑81 (MQ=255)                            aaGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAActgcctg                                           >  2:951096/1‑80 (MQ=255)                             aGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAActgcctgc                                          >  2:951097/1‑80 (MQ=255)                             aGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCa                                         <  1:951061/81‑1 (MQ=255)                               tGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAAt                                       <  1:951062/81‑1 (MQ=255)                                   ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCg                                   >  1:951101/1‑81 (MQ=255)                                   ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCg                                   >  1:951090/1‑81 (MQ=255)                                      ggTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGtt                                <  2:951063/81‑1 (MQ=255)                                           aaCCACTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGcgcg                           >  2:951109/1‑81 (MQ=255)                                               aCTCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCg                       <  2:951064/81‑1 (MQ=255)                                                 tCCTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCggg                     >  2:951111/1‑81 (MQ=255)                                                   cTGTAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCgggg                    <  2:951065/80‑1 (MQ=255)                                                      tAGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACg                >  2:951112/1‑81 (MQ=255)                                                       aGCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGt               >  2:951129/1‑81 (MQ=255)                                                        gCTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGt               <  1:951066/80‑1 (MQ=255)                                                         cTCGACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTg             <  1:951067/81‑1 (MQ=255)                                                            gACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCa           <  1:951068/80‑1 (MQ=255)                                                            gACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCAc          >  1:951117/1‑81 (MQ=255)                                                            gACCAGCCCTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCAc          >  1:951107/1‑81 (MQ=255)                                                                   ccTTTACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAg   <  2:951069/81‑1 (MQ=255)                                                                     tttACGGGTAAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGt  >  2:951120/1‑80 (MQ=255)                                                                              |                                                                       CATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGT  >  NC_000913/4583580‑4583727

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