New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4584195 =122 (1.130)15 (0.150) 7/164 NT 11.9% coding (2567/3513 nt) hsdR endonuclease R Type I restriction enzyme
?NC_000913 4584232 = 109 (1.100)coding (2530/3513 nt) hsdR endonuclease R Type I restriction enzyme
Rejected: Frequency below cutoff threshold.

AACAGATGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4584300‑4584195
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgacacAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACGCA  >  NC_000913/4584232‑4584330
                                                                                                                                                                                                             
AACAGATGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGT                                                                                                           >  1:1704061/1‑100
AACAGATGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGT                                                                                                           >  1:576740/1‑100
 ACAGATGAAAGGCCGCGCCACGCGCTTATGTCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTC                                                                                                          <  1:829463/100‑1
  CAGATGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCG                                                                                                         <  2:431990/100‑1
  CAGATGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCG                                                                                                         <  1:2017397/100‑1
    GATGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGAC                                                                                                       >  2:420716/1‑100
    GATGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGAC                                                                                                       <  2:341959/100‑1
    GATGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGAC                                                                                                       <  1:2870934/100‑1
      TGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGACTAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACAC                                                                                                     >  1:2140775/1‑100
      TGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACAC                                                                                                     >  1:1969858/1‑100
      TGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACAC                                                                                                     >  2:1568088/1‑100
      TGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACAC                                                                                                     >  1:2974817/1‑100
      TGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACAC                                                                                                     >  1:2974883/1‑100
           GGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTTACTTTTTCTATATCTAAAGCACGCTAGAGAGCGTCGACAacat                                                                                                  >  2:30776/1‑94
                                                         AGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGG                                                  >  2:440088/1‑100
                                                                TGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCC                                           >  1:1285876/1‑100
                                                                TGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCC                                           >  1:2345093/1‑100
                                                                TGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCC                                           >  1:2871431/1‑100
                                                                 GACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTTACGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCT                                          >  1:541975/1‑100
                                                                 GACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCT                                          >  1:2922747/1‑100
                                                                 GACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCT                                          >  2:2072841/1‑100
                                                                 GACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCT                                          >  2:2175175/1‑100
                                                                 GACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCT                                          >  2:2847718/1‑100
                                                                 GACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAAGCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCC                                           >  2:598736/1‑99
                                                                  ACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTT                                         >  2:2593206/1‑100
                                                                  ACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACATCCGGGCATAAGCGCGTGGCGCGGCCTT                                         >  1:1749977/1‑100
                                                                   CTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTT                                        >  1:1071747/1‑100
                                                                    TGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTC                                       >  2:2249503/1‑100
                                                                          GATATCTACAGCACGCTGGAGAGCGTCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGT                                 >  2:1488250/1‑100
                                                                                                      ACACAGTCAAAAATCTTAAAGCTGGTCTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTAC     >  1:850377/1‑100
                                                                                                        ACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACGC   <  1:2010368/100‑1
                                                                                                        ACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACG    >  2:342720/1‑99
                                                                                                         CAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACGCA  <  2:2176987/100‑1
                                                                                                         CAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACGCA  <  1:19642/100‑1
                                                                                                                                                                                                             
AACAGATGAAAGGCCGCGCCACGCGCTTATGCCCGGAGGTGAATAAAACCAGCTTTAAGATTTTTGACTGTGTCGATATCTACAGCACGCTGGAGAGCGTCGACAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4584300‑4584195
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgacacAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACGCA  >  NC_000913/4584232‑4584330

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