Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,583,432 Δ3,267 bp coding (64‑3330/3513 nt) hsdR ← type I restriction enzyme EcoKI endonuclease subunit

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4583432 4586698 3267 5 [0] [0] 5 hsdR type I restriction enzyme EcoKI endonuclease subunit

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 45834310 (0.000)5 (0.680) 4/326 0.3 100% coding (3331/3513 nt) hsdR type I restriction enzyme EcoKI endonuclease subunit
?NC_000913 4586699 = 0 (0.000)coding (63/3513 nt) hsdR type I restriction enzyme EcoKI endonuclease subunit

AGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTCCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGG  >  NC_000913/4583223‑4583615
                                                                                                                                                                                                                |                                                                                                                                                                                        
aGCCGGAAGGAGGAAGTGTATAGGTGTCAGGCCAGCTCGTCCCAGATAGAATCGCTGAATTTGCCCAGCAGGGTTTCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGGCGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTGTCTCTGTCAGCGCCTGCGCTAAACGATCGAGGCAGCTTAAggcgatggcataagtgaagtcgttgacgcccttcaggaattc                                                                                                                                                <  1:76095‑M1/251‑43 (MQ=255)
                         gtCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGAATAGAGATTATAGTCCATGGTTATTTGCAGCCTCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGTAGAGATCGTCGTCGCGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcataagtgaagtcgttgacgcccttcaggaattcaaaattggatttattcatccttgtt                                                                                                                       >  1:110537‑M1/1‑184 (MQ=255)
                                                                                                                               ggAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcataagtgaagtcgttgacgcccttcaggaattcaaaattggatttattcatcattgttattaatccattgctgtgcgggcctgtccaaatatttaagn                                                                               <  1:94867‑M1/189‑108 (MQ=255)
                                                                                                                               ggAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcataagtgaagtcgttgacgcccttcaggaattcaaaattgganttannnnncann                                                                                                                          >  2:94867‑M1/1‑82 (MQ=255)
                                                                                                                                               gAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcataagtgaagtcgttgacgcccttcaggaattcaaaattggatttattcatcattgttattaatccatttctttgctggcctgtccaaatatttaaggcccataacatctcatcttagctttctgtacctttccggtcaatgaccacggtcacagcaactgactcatttctaacg  <  1:152410‑M1/250‑185 (MQ=255)
                                                                                                                                                                                                                |                                                                                                                                                                                        
AGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTCCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGG  >  NC_000913/4583223‑4583615

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