Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,582,054 C→T C332Y (TGT→TAT)  hsdM ← type I restriction enzyme EcoKI methylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,582,0540CT100.0% 11.8 / NA 7C332Y (TGT→TAT) hsdMtype I restriction enzyme EcoKI methylase subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (2/5);  total (2/5)

GCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCG  >  NC_000913/4581860‑4582290
                                                                                                                                                                                                  |                                                                                                                                                                                                                                            
gCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGGCATACCCCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGATACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACg                                                                                                                                                                                      <  1:102145/251‑1 (MQ=255)
                                                                                                    cGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGATACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAgg                                                                                  >  1:1107/1‑251 (MQ=255)
                                                                                                                                   ttGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGATACTTATCCATCAGGTCACGACGAATGTCgg                                                                                                                                                                                                                <  1:122573/94‑1 (MQ=255)
                                                                                                                                                      cGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGATACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATAGTCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCggg                                                                            <  1:157388/207‑1 (MQ=255)
                                                                                                                                                                  cGGTCGGCAGACGCAGAATGGTGTGCAGATGATACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGcc                    >  1:160494/1‑251 (MQ=255)
                                                                                                                                                                                     ggTGTGCAGATTATACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCCCCACCCCCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTTTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCg  <  1:90756/250‑1 (MQ=255)
                                                                                                                                                                                      gtggccAGATGATCCTTATCCATCAGGTCACGAAGAATTTCGGTGCCTTCACCTCCTGCAAAGAGCACGTTATCCGCCACCACCCCCGCCGCACGACCGCCGGGATGCATCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGTTGAACAAAGGTGCGTGTAATGTTGGTGCCTGCGGCGCTGc                                                    <  1:213301/194‑1 (MQ=255)
                                                                                                                                                                                                  |                                                                                                                                                                                                                                            
GCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAAAGGTGCGGGTAATGTTGGTGCCTGCGGCGCTGCCAAACGGCGGGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCG  >  NC_000913/4581860‑4582290

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: