Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,584,630	A→G	V711A (GTC→GCC)	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,630	0	A	G	100.0%	149.6 / NA	43	V711A (GTC→GCC)	hsdR	type I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/27);  total (16/27)

GCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCAAAATCCTGATCGTCTTCCAGGGT  >  NC_000913/4584552‑4584706                                                                               |                                                                             gCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCg                                                                            <  1:951288/81‑1 (MQ=255)   aTGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAggcggc                                                                          <  2:951289/81‑1 (MQ=255)      ggCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGcgcg                                                                       >  2:951353/1‑81 (MQ=255)                aCGCAGAAGACCAGCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGc                                                              <  1:951290/80‑1 (MQ=255)                    agaagaCCAGCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGa                                                         <  1:951291/81‑1 (MQ=255)                       agaCCAGCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCa                                                      >  2:951366/1‑81 (MQ=255)                        gaCCAGCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCAc                                                     >  1:951374/1‑81 (MQ=255)                         aCCAGCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCAc                                                     >  2:951324/1‑80 (MQ=255)                           cAGCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAg                                                  >  2:951364/1‑81 (MQ=255)                           cAGCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAg                                                  >  1:951355/1‑81 (MQ=255)                             gCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGc                                                <  1:951292/81‑1 (MQ=255)                              cGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGc                                                >  1:951393/1‑80 (MQ=255)                              cGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGc                                                >  1:951370/1‑80 (MQ=255)                               gTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCa                                              <  2:951293/81‑1 (MQ=255)                               gTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCa                                              >  1:951362/1‑81 (MQ=255)                                ttttttGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCAc                                             <  1:951294/81‑1 (MQ=255)                                ttttttGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCAc                                             <  2:951295/81‑1 (MQ=255)                                  ttttGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACgg                                           >  1:951384/1‑81 (MQ=255)                                     tGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGtt                                         >  2:951385/1‑80 (MQ=255)                                      gCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTa                                        <  2:951296/80‑1 (MQ=255)                                       cGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaaa                                      <  1:951298/81‑1 (MQ=255)                                       cGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa                                       <  2:951297/80‑1 (MQ=255)                                        gATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAg                                     <  1:951299/81‑1 (MQ=255)                                          tCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTc                                   >  2:951375/1‑81 (MQ=255)                                          tCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTc                                   <  2:951300/81‑1 (MQ=255)                                             ggTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGc                                <  2:951301/81‑1 (MQ=255)                                              gTCGGGTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCg                               <  2:951302/81‑1 (MQ=255)                                                 gggTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGAc                             <  1:951303/80‑1 (MQ=255)                                                   gTCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTc                          <  2:951304/81‑1 (MQ=255)                                                    tCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCa                         <  2:951306/81‑1 (MQ=255)                                                    tCAAGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCa                         <  1:951305/81‑1 (MQ=255)                                                      aaGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCaaa                       <  2:951307/81‑1 (MQ=255)                                                      aaGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCaaa                       <  1:951308/81‑1 (MQ=255)                                                       aGATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCaaaa                      >  1:951383/1‑81 (MQ=255)                                                        gATAATTGGTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCAAAAt                     <  1:951309/81‑1 (MQ=255)                                                                gTGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCAAAATCCTGATCg             <  1:951310/81‑1 (MQ=255)                                                                 tGAGTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCAAAATCCTGATCGt            <  2:951311/81‑1 (MQ=255)                                                                    gTTCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCAAAATCCTGATCGTCtt         <  2:951312/81‑1 (MQ=255)                                                                     ttCGTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCAAAATCCTGATCGTCTTc        <  1:951313/81‑1 (MQ=255)                                                                        gTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCAAAATCCTGATCGTCTTCCa      <  2:951314/80‑1 (MQ=255)                                                                        gTTACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCAAAATCCTGATCGTCTTCCAg     >  2:951380/1‑81 (MQ=255)                                                                         ttACAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCAAAATCCTGATCGTCTTCCAgg    >  1:951373/1‑81 (MQ=255)                                                                           aCAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCAAAATCCTGATCGTCTTCCAGGgt  >  2:951372/1‑81 (MQ=255)                                                                            cAGGCGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCAAAATCCTGATCGTCTTCCAGGgt  <  1:951315/80‑1 (MQ=255)                                                                               |                                                                             GCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCGGGTCAAGATAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAAAGTCGGCGACTTCAAAATCCTGATCGTCTTCCAGGGT  >  NC_000913/4584552‑4584706

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: ‑