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

GGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4583195‑4583431
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTT  >  NC_000913/4586699‑4586895
                                                                                                                                                                                                                                                                                                                                                                                                                                                  
GGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAGGCGATGGCATA                                                                                                                                                                                           >  1:102248/1‑249
                                                     GTCAGGCCAGCTCGTCCCAGATATAAGCGCTGAATTTGCCCAGCAGGGTATCGAGATGATCTTCAAAGGTTCTGTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGATCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAGGC                                                                                                                                                                                                    <  1:247438/187‑1
                                                                      CAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGCTCCCGGTTTTGAAGACATCGTCGTCGAGCACCAATTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAGGCGATGGCATAAGTGAAGTCGTTGACGCACTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTG                                                                                                                   >  1:407590/1‑251
                                                                                                                                             CCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATC                                                                                 >  1:207023/1‑214
                                                                                                                                                GCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGGCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAGGCGATGGCATAAGTGAAGGCGTTGACGCCCTTCAGTAATTCAAAATTGGATTTATTCATCATTGTTAGTAATCCATTGCTGTGCGGGCCTGTTCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGAC                                              <  1:279182/246‑1
                                                                                                                                                           GGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCCCTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGACTTCAAAATTGGCTTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAAc                                                 >  1:67773/1‑231
                                                                                                                                                                                    GTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATGTAAGGCCCATAACATCTCATCTTAGC                                                                       <  1:143007/185‑1
                                                                                                                                                                                       GTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAGGCGATGGCATAAGTGAAGTAGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGCCCAAATATTTAATGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCACTGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTT  >  1:300586/1‑251
                                                                                                                                                                                                       TCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGNNTTATTCNNC                                                                                                                                          >  2:25385/1‑99
                                                                                                                                                                                                       TCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTAC                                                              <  1:25385/175‑1
                                                                                                                                                                                                                            CGATCGGNCCAGCTTAAGGCGATGGCATAAGTGAAGTCGNTNACGCCCTTCAGGAATGCAAAATTGGATTTATTCATCATTGTTATTAATNCAGTGCTGTGNGGGNNTGTCNANANATTTAAGGCCCATAACATCTCATCTTAGCNTTCTGTACCTTTCCGGGCAATGANNACNNNCANNGCAAC                               <  2:317222/185‑1
                                                                                                                                                                                                                            CGATCGAGCCAGCTTAAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGACTGTCCAAATATTTAAGGCCCATAACATCTCCTCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAAC                               >  1:317222/1‑185
                                                                                                                                                                                                                                 GAGCCAGCTTAAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCAC                                           >  1:3899/1‑168
                                                                                                                                                                                                                                                                                                                                                                                                                                                  
GGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4583195‑4583431
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCACAGCAACTGACTCATTTCTAACGTGTTCGTCTATTT  >  NC_000913/4586699‑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: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.