New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 539537432 (1.010)3 (0.100) 3/228 NT 8.8% coding (2268/2364 nt) hsdR type I restriction modification system endonuclease
?NC_002947 = 5395392 33 (1.130)coding (2286/2364 nt) hsdR type I restriction modification system endonuclease

GTAAAGAACGCGCCGAACAAGTGAAAAAGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5395217‑5395374
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctttaaccagTTCGATGGGTGCACCGAGCTGGTTAAAGGGATCGAGCTGGAGGATTTTGATGTCCTCGATGTGCTCAATGCCGGTGTCGGCGTATTTATCCAGCAGCGCTTCCAGCACCTGGCGCGCTGCACCGGAGTATTTGGCGAAGTAGTTAC  <  NC_002947/5395392‑5395247
                                                                                                                                                                                                                                                                                                                
GTAAAGAACGCGCCGAACAAGTGAAAAAGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCT                                                                                                                                                            <  2:446163/150‑1
   AAGAACGCGCCGAACAAGTGAAAAAGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACAGCCAGGACATCAAAATCCTCCAGCTCGCTCCATTTA                                                                                                                                                         >  2:779854/1‑150
     GAACGCGCCGAACAAGTGAAAAAGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAAC                                                                                                                                                       <  2:474733/150‑1
               AACAAGTGAAAAAGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAA                                                                                                                                                        >  1:470508/1‑139
               AACAAGTGAAAAAGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAA                                                                                                                                                        <  2:470508/139‑1
                                                                 CAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGTTCGATGGGTGCACCGAGCTGGTTAAAGGGA                                                                                                                     <  1:712027/124‑1
                                                                 CAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGTTCGATGGGTGCACCGAGCTGGTTAAAGGGA                                                                                                                     >  2:712027/1‑124
                                                                                      GATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAAC                                                                                                                                                       <  1:782319/69‑1
                                                                                      GATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAAC                                                                                                                                                       >  2:782319/1‑69
                                                                                                                                           GCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGGGATCGAGCTGGAGGATTTTGATGTCCTCGATGTGCTCAATGCCGGTGTCGGCGTATTTATCCAGCAGCGCTTCCAGCACCTGGCGCGCTGCACCGGAGTATT                 >  2:369609/1‑150
                                                                                                                                                     TTTAACCAGTTCGATGGGTGCACCGAGCTGGTTAAAGGGATCGAGCTGGAGGATTTTGATGTCCTCGATGTGCTCAATGCCGGTGTCGGCGTATTTATCCAGCAGCGCTTCCAGCACCTGGCGCGCTGCAC                          >  1:111148/1‑131
                                                                                                                                                     TTTAACCAGTTCGATGGGTGCACCGAGCTGGTTAAAGGGATCGAGCTGGAGGATTTTGATGTCCTCGATGTGCTCAATGCCGGTGTCGGCGTATTTATCCAGCAGCGCTTCCAGCACCTGGCGCGCTGCAC                          <  2:111148/131‑1
                                                                                                                                                          CCAGTTCGATGGGTGCACCGAGCTGGTTAAAGGGATCGAGCTGGAGGATTTTGATGTCCTCGATGTGCTCAATGCCGGTGTCGGCGTATTTATCCAGCAGCGCTTCCAGCACCTGGCGCGCTGCACCGGAGTATTTGGCGAAGTAGTTAC  <  2:683154/150‑1
                                                                                                                                                                                                                                                                                                                
GTAAAGAACGCGCCGAACAAGTGAAAAAGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5395217‑5395374
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctttaaccagTTCGATGGGTGCACCGAGCTGGTTAAAGGGATCGAGCTGGAGGATTTTGATGTCCTCGATGTGCTCAATGCCGGTGTCGGCGTATTTATCCAGCAGCGCTTCCAGCACCTGGCGCGCTGCACCGGAGTATTTGGCGAAGTAGTTAC  <  NC_002947/5395392‑5395247

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ 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.