New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4585317116 (1.240)6 (0.070) 5/164 NT 5.5% coding (1445/3513 nt) hsdR endonuclease R Type I restriction enzyme
?NC_000913 = 4585344 99 (1.150)coding (1418/3513 nt) hsdR endonuclease R Type I restriction enzyme
Rejected: Frequency below cutoff threshold.

TATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4585213‑4585317
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atggcgaCCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTC  <  NC_000913/4585344‑4585248
                                                                                                                                                                                                          
TATCTTCAAACGCGCCCAGCGCCTGTTCGTGAAGAGCACGGCCGTCGACAATGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCA                                                                                                         <  2:885861/99‑1
  TCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGG                                                                                                      >  2:737188/1‑100
   CTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGC                                                                                                     >  1:116257/1‑100
   CTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGC                                                                                                     >  2:332681/1‑100
   CTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGC                                                                                                     >  2:1717072/1‑100
    TTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGCAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCG                                                                                                    >  2:1185077/1‑100
    TTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCG                                                                                                    >  1:1530810/1‑100
    TTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCG                                                                                                    >  1:698051/1‑100
    TTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCG                                                                                                    >  1:409598/1‑100
                       TGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGACCGGTACCGGTAAAACCC                                                                                 >  2:2363511/1‑100
                                      CGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCA                                                                  <  1:937445/100‑1
                                       GGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCAT                                                                 <  1:2363524/100‑1
                                          GGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGAT                                                              <  2:2390808/100‑1
                                          GGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGAT                                                              <  1:155300/100‑1
                                            TCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGT                                                            <  1:166298/100‑1
                                                                                                   TGGCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCG     >  1:1728836/1‑100
                                                                                                     GCGACCGGTACCGGTAAAACCCGTACGGCAATCGCCACGATGTTCCGCCTGAACCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCG     >  2:1232298/1‑98
                                                                                                      CGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTC  >  1:155195/1‑100
                                                                                                      CGACCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCGCTTCCTTGTCGACCGCCGTTC  <  2:2462184/100‑1
                                                                                                                                                                                                          
TATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4585213‑4585317
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atggcgaCCGGTACCGGTAAAACCCGTACGGCAATCGCCATGATGTTCCGCCTGATCCAGTCCCAGCGTTTTAAACGCATTCTCTTCCTTGTCGACCGCCGTTC  <  NC_000913/4585344‑4585248

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