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 13 [0] [0] 13 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)13 (0.780) 13/292 0.2 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

ATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTT  >  NC_000913/4586462‑4586895
                                                                                                                                                                                                                                             |                                                                                                                                                                                                    
ggtgaacgagcgcagccaacgcagaggcagcctgaaggatgaagtgtatacgtgtcaggccagctcgtcccagatataatcgctgaatttgcccagcagggtatcgagattatcgtcaaaggttctttgcagcatcgccttcccgccgcgacggtggaagttgccggttttgaagacatcgtcgtcgagcaccactttctctttcagcgcctgcgctaaacgatcgagccagcttaaGGCGATGGCATa                                                                                                                                                                                           >  1:102248‑M2/238‑249 (MQ=255)
                                                     gtcaggccagctcgtcccagatataagcgctgaatttgcccagcagggtatcgagatgatcttcaaaggttctgtgcagcatcgccttcccgccgcgacggtggaagttgccggttttgaagacatcgtcgtcgatcaccactttctctttcagcgcctgcgctaaacgatcgagccagcttaaGGc                                                                                                                                                                                                    <  1:247438‑M2/3‑1 (MQ=255)
                                                                      cagatataatcgctgaatttgcccagcagggtatcgagattatcgtcaaaggttctttgcagcatcgccttcccgccgcgacggtggaagctcccggttttgaagacatcgtcgtcgagcaccaatttctctttcagcgcctgcgctaaacgatcgagccagcttaaGGCGATGGCATAAGTGAAGTCGTTGACGCACTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCtgtg                                                                                                                   >  1:407590‑M2/168‑251 (MQ=255)
                                                                                                                                             cccgccgcgacggtggaagttgccggttttgaagacatcgtcgtcgagcaccactttctctttcagcgcctgcgctaaacgatcgagccagcttaaGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatc                                                                                 >  1:207023‑M2/97‑214 (MQ=255)
                                                                                                                                                gccgcgacggtggaagttgccggttttgaagacatcgtcggcgagcaccactttctctttcagcgcctgcgctaaacgatcgagccagcttaaGGCGATGGCATAAGTGAAGGCGTTGACGCCCTTCAGTAATTCAAAATTGGATTTATTCATCATTGTTAGTAATCCATTGCTGTGCGGGCCTGTTCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGAc                                              <  1:279182‑M2/153‑1 (MQ=255)
                                                                                                                                                           ggaagttgccggttttgaagacatcgtcgtcgagcacccctttctctttcagcgcctgcgctaaacgatcgagccagcttaaGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGACTTCAAAATTGGCTTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGcaac                                                 >  1:67773‑M2/83‑231 (MQ=255)
                                                                                                                                                                                    gtcgtcgagcaccactttctctttcagcgcctgcgctaaacgatcgagccagcttaaGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATGTAAGGCCCATAACATCTCATCTTAGc                                                                       <  1:143007‑M2/128‑1 (MQ=255)
                                                                                                                                                                                       gtcgagcaccactttctctttcagcgcctgcgctaaacgatcgagccagcttaaGGCGATGGCATAAGTGAAGTAGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGCCCAAATATTTAATGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCACTGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTAttt  >  1:300586‑M2/55‑251 (MQ=255)
                                                                                                                                                                                                       tctttcagcgcctgcgctaaacgatcgagccagcttaaGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGNNTTATtcnnc                                                                                                                                          >  2:25385‑M2/39‑99 (MQ=255)
                                                                                                                                                                                                       tctttcagcgcctgcgctaaacgatcgagccagcttaaGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTAc                                                              <  1:25385‑M2/137‑1 (MQ=255)
                                                                                                                                                                                                                            cgatcggnccagcttaaGGCGATGGCATAAGTGAAGTCGNTNACGCCCTTCAGGAATGCAAAATTGGATTTATTCATCATTGTTATTAATNCAGTGCTGTGNGGGNNTGTCNANANATTTAAGGCCCATAACATCTCATCTTAGCNTTCTGTACCTTTCCGGGCAATGANNACNNNCANNGCAAc                               <  2:317222‑M2/168‑1 (MQ=255)
                                                                                                                                                                                                                            cgatcgagccagcttaaGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGACTGTCCAAATATTTAAGGCCCATAACATCTCCTCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAAc                               >  1:317222‑M2/18‑185 (MQ=255)
                                                                                                                                                                                                                                 gagccagcttaaGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCAc                                           >  1:3899‑M2/13‑168 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                    
ATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTT  >  NC_000913/4586462‑4586895

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