Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,584,914	A→G	A616A (GCT→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,914	0	A	G	100.0%	107.5 / NA	31	A616A (GCT→GCC)	hsdR	type I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/14);  total (17/14)

CGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACG  >  NC_000913/4584779‑4585046                                                                                                                                        |                                                                                                                                     cgcgTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGa                                                                                                                                  >  2:64732/1‑140 (MQ=255)      tGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGAttttt                                                                                                                             >  1:161215/1‑140 (MQ=255)                     gAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCAt                                                                                                                      <  1:179196/132‑1 (MQ=255)                     gAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCAt                                                                                                                      >  2:179196/1‑132 (MQ=255)                      aaTAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGAt                                                                                                             >  2:91872/1‑140 (MQ=255)                           gCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCgaga                                                                                                        >  2:387816/1‑140 (MQ=255)                                tCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCg                                                                                                   >  2:401459/1‑140 (MQ=255)                                tCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCg                                                                                                   >  2:304407/1‑140 (MQ=255)                                tCCTGGTCGATCAGAAAACCGTCGATAACCGCGGAACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCCGGGGGGCGGTGAGGGCGATTTTTACCGCATCGAAGCGATCCAGAATGCg                                                                                                   >  2:156481/1‑140 (MQ=255)                                 ccTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGa                                                                                                  <  1:57194/140‑1 (MQ=255)                                   tGGTCGATCAGAAAACCGTCGATAACCGCGGTACTGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACg                                                                                                >  1:64052/1‑140 (MQ=255)                                     gTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGt                                                                                              <  1:387816/140‑1 (MQ=255)                                        gATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAgg                                                                                           <  1:156481/140‑1 (MQ=255)                                           cAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCag                                                                                        <  1:64732/140‑1 (MQ=255)                                           cAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCag                                                                                        <  1:304407/140‑1 (MQ=255)                                                    gTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCag                                                                                        >  1:284699/1‑131 (MQ=255)                                                    gTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCag                                                                                        <  2:284699/131‑1 (MQ=255)                                                     tCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATc                                                                              >  1:300991/1‑140 (MQ=255)                                                     tCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATc                                                                              >  1:322455/1‑140 (MQ=255)                                                        aTAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAg                                                                           <  2:322455/140‑1 (MQ=255)                                                                            ataACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGtt                                                       <  1:415897/140‑1 (MQ=255)                                                                            ataACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGtt                                                       <  2:362909/140‑1 (MQ=255)                                                                                   tAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCTACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTc                                                >  2:72644/1‑140 (MQ=255)                                                                                   tAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTc                                                >  2:128258/1‑140 (MQ=255)                                                                                   tAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTc                                                >  1:105690/1‑140 (MQ=255)                                                                                      aCCGGCTCTCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGt                                             <  2:294657/140‑1 (MQ=255)                                                                                        cGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCt                                           <  2:64052/140‑1 (MQ=255)                                                                                                    aaTCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCga                               <  2:105690/140‑1 (MQ=255)                                                                                                                    aGCGCCGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATg               >  2:360327/1‑140 (MQ=255)                                                                                                                         cGGGGTGGCGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCtt          >  2:15799/1‑140 (MQ=255)                                                                                                                                 cGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCaa    >  1:356979/1‑138 (MQ=255)                                                                                                                                 cGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCaa    <  2:356979/138‑1 (MQ=255)                                                                                                                                 cGGTGAGGGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAacg  <  1:401459/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     CGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATCCAGCTGGCTGCGGAACTGCAGTTCGCCTTCGGTCTGCTCTTTATCGAGAATATAGCCGCGATGCGCTTCGTCAACG  >  NC_000913/4584779‑4585046

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