New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4585553 =128 (1.040)2 (0.020) 5/166 NT 1.8% coding (1209/3513 nt) hsdR endonuclease R Type I restriction enzyme
?NC_000913 4585576 = 101 (0.880)coding (1186/3513 nt) hsdR endonuclease R Type I restriction enzyme
Rejected: Frequency below cutoff threshold.

AAACAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCG‑ACGTGCGTGATACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4585657‑4585553
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gataccAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCC  >  NC_000913/4585576‑4585674
                                                                                                                                                                                                             
AAACAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATCCCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCG‑ACGTGCGTG                                                                                                          >  2:1432345/1‑100
AAACAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCG‑ACGTGCGTG                                                                                                          >  2:1802821/1‑100
AAACAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCG‑ACGTGCGTG                                                                                                          <  2:1608341/100‑1
   CAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGCCCAAAAGCGGCATCTGGTATCGCG‑ACGTGCGTGATA                                                                                                       >  2:2622369/1‑100
   CAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAACAGCGGCATCTGGTATCGCG‑ACGTGCGTGATA                                                                                                       >  2:3281976/1‑100
   CAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCG‑ACGTGCGTGATA                                                                                                       >  2:2802676/1‑100
     ACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGTCCAAAAGCGGCATCTGGTATCGCG‑ACGTGCGTGATACC                                                                                                     <  1:2873624/100‑1
     ACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCG‑ACGTGCGTGATACC                                                                                                     <  2:2109585/100‑1
        GTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCACGCACGT‑CGCGATACCAGA                                                                                                  <  1:1605474/100‑1
         TTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCACGCACGT‑CGCGATACCAGAT                                                                                                 <  1:1155869/100‑1
           TAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCACGCACGT‑CGCGATACCAGATGC                                                                                               <  2:73586/100‑1
            AAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCACGCACGT‑CGCGATACCAGATG                                                                                                <  1:388130/98‑1
                                                                                  GGTATCGCG‑ACGTGCGTGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGT                        >  2:1935982/1‑100
                                                                                                    GATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTC       <  1:3111340/100‑1
                                                                                                     ATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCC      <  1:708786/100‑1
                                                                                                     ATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCC      <  2:575764/100‑1
                                                                                                      TACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCT     <  1:866051/100‑1
                                                                                                      TACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCT     <  2:462967/100‑1
                                                                                                       ACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTG    >  2:837930/1‑100
                                                                                                       ACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTG    <  1:1180487/100‑1
                                                                                                        CCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAGCCGTTGTTTGCCGCTGGTGTCCTG    <  1:2216400/99‑1
                                                                                                        CCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGC   <  2:2594580/100‑1
                                                                                                        CCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGC   <  2:2594367/100‑1
                                                                                                        CCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGC   >  1:3292278/1‑100
                                                                                                         CAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCC  >  2:318046/1‑100
                                                                                                         CAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCC  <  1:1263013/100‑1
                                                                                                                                                                                                             
AAACAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCG‑ACGTGCGTGATACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4585657‑4585553
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gataccAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCC  >  NC_000913/4585576‑4585674

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