Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0131982,163,4370G→T100.0% 2.1 2N697K (AAC→AAALGG_02097type II restriction enzyme methyltransferase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/2):  ref base (0/0):  total (0/2)

CGATATGATCCACCCTTGTTGTAAGGAAACCATTT.CTTGATCGACCTAAGCGATTCTGCGGTGCTCTTTGTGTTAAAACTAATATTATGCACATTTACTTCA  >  NC_013198/2163366‑2163467
                                                                        |                              
cGATATGTCCCCCCCTGGTT.TAAAGAAACCCTTTTCTTAATGGACTTAAGCGATTC.GGGGGGCTCTTTGTTTTTAAACTAATATTATGCACATTTACTTCa  <  11:491447/101‑1
cGATATGTCCCCCCCTGGTT.TAAAGAAACCCTTTTCTTAATGGACTTAAGCGATTC.GGGGGGCTCTTTGTTTTTAAACTAATATTATGCACATTTACTTCa  <  26:491447/101‑1
                                                                        |                              
CGATATGATCCACCCTTGTTGTAAGGAAACCATTT.CTTGATCGACCTAAGCGATTCTGCGGTGCTCTTTGTGTTAAAACTAATATTATGCACATTTACTTCA  >  NC_013198/2163366‑2163467

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 18 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG