Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,584,768	A→C	V665G (GTT→GGT)	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,584,768	0	A	C	100.0%	118.6 / NA	33	V665G (GTT→GGT)	hsdR	type I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (17/16);  total (17/16)

GATCGTCTTCCAGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTA  >  NC_000913/4584691‑4584845                                                                              |                                                                              gATCGTCTTCCAGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAcccc                                                                            >  2:871108/1‑81 (MQ=255)       cTTCCAGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCcctcct                                                                      >  1:871115/1‑81 (MQ=255)          ccAGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCcctcct                                                                      <  1:871059/78‑1 (MQ=255)          ccAGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTgc                                                                    <  1:871058/80‑1 (MQ=255)          ccAGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTg                                                                     <  2:871057/79‑1 (MQ=255)           cAGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTgcgc                                                                  >  1:871096/1‑81 (MQ=255)            aGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGcgcg                                                                 >  1:871098/1‑81 (MQ=255)            aGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGcgcg                                                                 <  2:871060/81‑1 (MQ=255)            aGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGcgcg                                                                 >  1:871113/1‑81 (MQ=255)              ggTGTCATTGATCACTTCTCCCTGCGGGCTGGTGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGtt                                                               <  2:871061/81‑1 (MQ=255)                tgtCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGc                                                             <  1:871062/81‑1 (MQ=255)                        aTCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGtgat                                                      <  2:871063/80‑1 (MQ=255)                         tCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGtgatga                                                    <  2:871064/81‑1 (MQ=255)                           aCTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATc                                                  >  2:871110/1‑81 (MQ=255)                            cTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCt                                                 <  1:871065/81‑1 (MQ=255)                               ctcCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTg                                                <  1:871066/79‑1 (MQ=255)                               ctcCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGaa                                              >  2:871109/1‑81 (MQ=255)                                 cccTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGaa                                              >  1:871107/1‑79 (MQ=255)                                      cGGGCTGATTCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCg                                        >  2:871132/1‑80 (MQ=255)                                        ggCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGAt                                     >  1:871104/1‑81 (MQ=255)                                         gCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATc                                    >  1:871125/1‑81 (MQ=255)                                         gCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATc                                    <  2:871067/81‑1 (MQ=255)                                                cgcTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTc                              <  2:871068/80‑1 (MQ=255)                                                 gcTCTACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGa                            >  1:871126/1‑81 (MQ=255)                                                    ctACCTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGAt                           >  2:871112/1‑79 (MQ=255)                                                       ccTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGa                       >  2:871155/1‑80 (MQ=255)                                                       ccTGCTCGCCTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGa                       <  1:871069/80‑1 (MQ=255)                                                               ccTTTGGAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCg              <  1:871070/81‑1 (MQ=255)                                                                    ggAGAGATAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAAc         >  2:871119/1‑81 (MQ=255)                                                                     gagagaTAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAAcc        <  1:871071/81‑1 (MQ=255)                                                                      agagaTAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACcg       >  2:871128/1‑81 (MQ=255)                                                                        agaTAACCCCCCTCCTTCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGAGAACcgcg     >  2:871111/1‑81 (MQ=255)                                                                          aTAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGcggt   <  1:871072/81‑1 (MQ=255)                                                                           tAACCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGcggta  >  2:871127/1‑81 (MQ=255)                                                                              |                                                                              GATCGTCTTCCAGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTA  >  NC_000913/4584691‑4584845

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