Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,584,974	C→T	Q596Q (CAG→CAA)	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,974	0	C	T	100.0%	159.8 / NA	44	Q596Q (CAG→CAA)	hsdR	type I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/17);  total (27/17)

CGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGA  >  NC_000913/4584896‑4585050                                                                               |                                                                             cgcCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAgttg                                                                            <  1:951398/81‑1 (MQ=255)    cGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAgttggct                                                                         <  2:951399/81‑1 (MQ=255)    cGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAgttggct                                                                         <  2:951400/81‑1 (MQ=255)    cGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAgttggct                                                                         >  1:951478/1‑81 (MQ=255)      gggTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGc                                                                       >  1:951505/1‑81 (MQ=255)         tGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCgg                                                                     >  1:951460/1‑80 (MQ=255)         tGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGa                                                                    <  1:951402/81‑1 (MQ=255)         tGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGa                                                                    <  2:951401/81‑1 (MQ=255)              gTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTg                                                                >  1:951464/1‑80 (MQ=255)               tGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTg                                                                >  1:951403/1‑79 (MQ=255)                gagaGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAg                                                             >  2:951479/1‑81 (MQ=255)                 agagCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAg                                                             <  2:951404/80‑1 (MQ=255)                   agCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTc                                                          <  2:951406/81‑1 (MQ=255)                   agCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTc                                                          <  2:951405/81‑1 (MQ=255)                   agCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTc                                                          <  1:951408/81‑1 (MQ=255)                   agCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTc                                                          <  1:951407/81‑1 (MQ=255)                      gATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGcc                                                       >  1:951480/1‑81 (MQ=255)                      gATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGcc                                                       >  1:951485/1‑81 (MQ=255)                              ccGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCt                                               <  1:951409/81‑1 (MQ=255)                               cGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTg                                              >  2:951465/1‑81 (MQ=255)                                 cATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGct                                            >  2:951452/1‑81 (MQ=255)                                 cATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGct                                            >  1:951486/1‑81 (MQ=255)                                  aTCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGctc                                           <  2:951410/81‑1 (MQ=255)                                    cGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCtt                                         >  2:951483/1‑81 (MQ=255)                                       aGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTAt                                      <  2:951411/81‑1 (MQ=255)                                        gTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATc                                     >  1:951461/1‑81 (MQ=255)                                           aTCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCgag                                  <  2:951412/81‑1 (MQ=255)                                            tCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCgag                                  >  2:951475/1‑80 (MQ=255)                                            tCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCgag                                  >  1:951459/1‑80 (MQ=255)                                              gagaATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGaa                                >  1:951467/1‑80 (MQ=255)                                              gagaATGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGaa                                >  2:951469/1‑80 (MQ=255)                                                 aaTGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAata                              >  2:951482/1‑79 (MQ=255)                                                  aTGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAg                           >  1:951488/1‑81 (MQ=255)                                                   tGCGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGc                          >  1:951489/1‑81 (MQ=255)                                                     cGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCcg                        >  1:951474/1‑81 (MQ=255)                                                     cGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCcg                        >  1:951470/1‑81 (MQ=255)                                                     cGACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCcg                        <  2:951413/81‑1 (MQ=255)                                                      gACGGTAGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCcgc                       >  1:951476/1‑81 (MQ=255)                                                            aGGCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATgcg                 >  2:951462/1‑81 (MQ=255)                                                             ggCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATgcgc                >  1:951493/1‑81 (MQ=255)                                                              gCAGAGACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCt               <  2:951414/81‑1 (MQ=255)                                                                agagACGTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTc             <  1:951415/81‑1 (MQ=255)                                                                      gTAATCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAac       >  1:951503/1‑81 (MQ=255)                                                                        aaTCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAacga     >  2:951457/1‑81 (MQ=255)                                                                         aTCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAacgac    <  2:951417/81‑1 (MQ=255)                                                                         aTCCAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAacgac    <  1:951416/81‑1 (MQ=255)                                                                            cAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAAcgacga  <  2:951418/80‑1 (MQ=255)                                                                            cAGTTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAAcgacga  >  1:951487/1‑80 (MQ=255)                                                                               |                                                                             CGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACGACGA  >  NC_000913/4584896‑4585050

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