Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,583,863	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,863	0	C	T	100.0%	68.8 / NA	26	R967H (CGT→CAT)	hsdR	endonuclease R Type I restriction enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/12);  total (14/12)

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

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