New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4585317143 (0.780)4 (0.020) 4/242 NT 2.9% coding (1445/3513 nt) hsdR endonuclease R Type I restriction enzyme
?NC_000913 = 4585344 137 (0.790)coding (1418/3513 nt) hsdR endonuclease R Type I restriction enzyme
Rejected: Frequency below/above cutoff threshold.

AGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4585164‑4585317
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atggcgaCCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTCTCTTGGCGAACAGGCGCTGGGCGCGTTTGAAGATACGCGTATTAACGGCGA  <  NC_000913/4585344‑4585197
                                                                                                                                                                                                                                                                                                              
AGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCAT                                                                                                                                                           <  4:236271/149‑1
  CCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCTCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGG                                                                                                                                                         >  4:301933/1‑149
  CCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGG                                                                                                                                                         >  4:275561/1‑149
   CCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGC                                                                                                                                                        <  4:1496918/149‑1
    CTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCG                                                                                                                                                       >  1:729512/1‑149
     TTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGA                                                                                                                                                      >  1:881786/1‑149
     TTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGA                                                                                                                                                      <  4:21914/149‑1
     TTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGA                                                                                                                                                      >  3:654064/1‑149
     TTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGA                                                                                                                                                      <  4:1092550/149‑1
     TTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGCAACATCATGGCGA                                                                                                                                                      >  3:618315/1‑149
                                         AATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGA                                                                                                                  <  2:1517560/149‑1
                                               CGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCC                                                                                                            <  3:808039/149‑1
                                                         AACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCA                                                                                                                     <  1:925763/130‑1
                                                         AACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCA                                                                                                                     >  2:925763/1‑130
                                                                                                                                                    TGGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTCTCTTGGCGAACAGGCGCTGGGCGCGTTTGAAGATACGCGTATTAAC       <  2:1519073/149‑1
                                                                                                                                                     GGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTCTCTTGGCGAACAGGCGCTGGGCGCGTTTGAAGATACGCGTATTAACG      >  2:1398208/1‑149
                                                                                                                                                     GGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTCTCTTGGCGAACAGGCGCTGGGCGCGTTTGAAGATACGCGTATTAACG      <  3:1337363/149‑1
                                                                                                                                                      GCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTCTCTTGGCGAACAGGCGCTGGGCGCGTTTGAAGATACGCGTATTAACGG     <  3:274202/149‑1
                                                                                                                                                       CGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTCTCTTGGCGAACAGGCGCTGGGCGCGTTTGAAGATACGCGTATTAACGGC    <  1:839614/149‑1
                                                                                                                                                        GACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTCTCTTGGCGAACAGGCGCTGGGCGCGTTTGAAGATACGCGTATTAACGGCG   >  2:1483476/1‑149
                                                                                                                                                         ACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTCTCTTGGCGAACAGGCGCTGGGCGCGTTTGAAGATACGCGTATTAACGGCGA  <  4:1547077/149‑1
                                                                                                                                                         ACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTCTCTTGGCGAACAGGCGCTGGGCGCGTTTGAAGATACGCGTATTAACGGCGA  <  3:311505/149‑1
                                                                                                                                                         ACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTCTCTTGGCGAACAGGCGCTGGGCGCGTTTGAAGATACGCGTATTAACGGCGA  <  2:438803/149‑1
                                                                                                                                                                                                                                                                                                              
AGCCCTTTAATGTCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4585164‑4585317
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atggcgaCCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTCTCTTGGCGAACAGGCGCTGGGCGCGTTTGAAGATACGCGTATTAACGGCGA  <  NC_000913/4585344‑4585197

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ 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.