Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,583,862	C→T	R967H (CGT→CAT)	hsdR ←	endonuclease R Type I restriction enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,583,862	0	C	T	100.0%	64.6 / NA	24	R967H (CGT→CAT)	hsdR	endonuclease R Type I restriction enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/10);  total (14/10)

GAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTAT  >  NC_000913/4583724‑4584000                                                                                                                                            |                                                                                                                                           gAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGc                                                                                                                                        <  1:226860/143‑1 (MQ=255)    tCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCCTCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGATGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGAt                                                                                                                                     >  2:392442/1‑143 (MQ=255)                  gAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAATCCAGGCAGtt                                                                                                                       >  2:355887/1‑143 (MQ=255)                       ccTGCGGCGTGTCGGAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTa                                                                                                                  <  1:415175/143‑1 (MQ=255)                         tggggcgTGTGGTAATCCCCGTACAGCGATTTTACGCTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTaaa                                                                                                                <  2:293262/140‑1 (MQ=255)                                      aaTCACCGTGCGGCGTGTTTGCGCGCGGCACTGCATCGTCGTTATCGCGGAAGATCGGCGCATCATTCAGGTTGTGGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTc                                                                                                   <  2:461767/143‑1 (MQ=255)                                        tCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCa                                                                                                 <  1:513076/143‑1 (MQ=255)                                                            acTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGCTGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCa                                                                             >  1:415379/1‑143 (MQ=255)                                                              tCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCACGAAGCTTGTTAAAGCATTCGGCGCTCCAGTGCGGCCCTTTTTCTCGCAgg                                                                           >  2:483626/1‑143 (MQ=255)                                                                     ttCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGCgg                                                                    >  1:348721/1‑143 (MQ=255)                                                                         tcgtcgATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGaa                                                                <  2:122262/143‑1 (MQ=255)                                                                               atatCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGtt                                                          >  2:276606/1‑143 (MQ=255)                                                                               atatCGAGGAAAATCGGCTCATCATT‑ACGTTTTTTTCGTCCCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGtt                                                          <  1:217007/143‑1 (MQ=255)                                                                                       gAAGATCGGCGCATCATTCAGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCGGGCGCGAGGCCAAGCCGTCAAAGTTCa                                                  >  1:348145/1‑143 (MQ=255)                                                                                                        tccggTTGTTGACGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTgg                                 <  1:3565/140‑1 (MQ=255)                                                                                                          aGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCa                               <  1:363991/143‑1 (MQ=255)                                                                                                                  gATGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCCTAGCTCAt                       >  1:385717/1‑143 (MQ=255)                                                                                                                   aTGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGCGGCGAACCCGTTAAAGTTCACGCCCGGCGCGTCCTGGGATCGCTccgc                      >  2:301091/1‑140 (MQ=255)                                                                                                                   aTGT‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCAGAGCCCagc                      >  2:402882/1‑141 (MQ=255)                                                                                                                      t‑CCGTTTTCAGCTTTTCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATcccg                   >  1:374366/1‑141 (MQ=255)                                                                                                                        cctttttCAGTTGTGCCAGTTGGGCCATAGACCCGTGTAGTTTGATATAGACTTCGGCGCTCCAGTGCGGCCCTTTTCCCCGCCGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAga                  <  2:58166/140‑1 (MQ=255)                                                                                                                                    tttCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCGTAGCTCATCCAGACGACGCACCTGt      >  1:398737/1‑143 (MQ=255)                                                                                                                                     ttCCAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACACCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGta     >  2:88646/1‑142 (MQ=255)                                                                                                                                        cAGATGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCAAGTGCGGCCCTTTTTCCCGCAGGCGCGCGGCGAAGCCGTTAAAGATCACGCCCGCCGCTGCCGGGCATAGCTCGTCCAGACGACGCCCCTGgttct  >  2:454026/1‑141 (MQ=255)                                                                                                                                            |                                                                                                                                           GAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGT‑CCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTAT  >  NC_000913/4583724‑4584000

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