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.

AAACAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCGACGTGCGTGATACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4585657‑4585553
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gataccAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCC  >  NC_000913/4585576‑4585674
                                                                                                                                                                                                             
AAACAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATCCCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCGACGTGCGTG                                                                                                          >  2:1432345/1‑100
AAACAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCGACGTGCGTG                                                                                                          >  2:1802821/1‑100
AAACAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCGACGTGCGTG                                                                                                          <  2:1608341/100‑1
   CAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGCCCAAAAGCGGCATCTGGTATCGCGACGTGCGTGATA                                                                                                       >  2:2622369/1‑100
   CAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAACAGCGGCATCTGGTATCGCGACGTGCGTGATA                                                                                                       >  2:3281976/1‑100
   CAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCGACGTGCGTGATA                                                                                                       >  2:2802676/1‑100
     ACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGTCCAAAAGCGGCATCTGGTATCGCGACGTGCGTGATACC                                                                                                     <  1:2873624/100‑1
     ACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCGACGTGCGTGATACC                                                                                                     <  2:2109585/100‑1
        GTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCACGCACGTCGCGATACCAGA                                                                                                  <  1:1605474/100‑1
         TTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCACGCACGTCGCGATACCAGAT                                                                                                 <  1:1155869/100‑1
           TAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCACGCACGTCGCGATACCAGATGC                                                                                               <  2:73586/100‑1
            AAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCACGCACGTCGCGATACCAGATG                                                                                                <  1:388130/98‑1
                                                                                  GGTATCGCGACGTGCGTGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGT                        >  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
                                                                                                                                                                                                             
AAACAACGGTTTAAAATCCCCTTCTGCTACTCGACCAACGGGCGCGAATACCGCGCAACAATGAAGACCAAAAGCGGCATCTGGTATCGCGACGTGCGTGATACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  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.