New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1203246 =23 (0.720)3 (0.120) 3/134 2.3 13.2% coding (290/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 1205075 = 21 (0.820)coding (40/540 nt) BW25113_RS06005 phage tail protein
Rejected: Frequency below/above cutoff threshold.

TACTGATTGCTCAATACCACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCCCAAACCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/1203349‑1203246
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAC  >  NZ_CP009273/1205075‑1205159
                                                                                                                                                                                             
TACTGATTGCTCAATACCACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCC                                                                                                     >  2:25261/1‑90
 ACTGATTGCTCAATACCACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCT                                                                                                    <  1:834876/90‑1
  CTGATTGCTCAATACCACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTT                                                                                                   <  2:387937/90‑1
      TTGCTCAATACCACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCC                                                                                               >  2:645585/1‑90
      TTGCTCAATACCACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCC                                                                                               >  2:831267/1‑90
        GCTCAATACCACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCCCA                                                                                             >  2:298205/1‑90
                                  ATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAA                                                                   <  2:414494/90‑1
                                               CGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAA                                                        >  1:4901/1‑88
                                               CGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAA                                                        <  2:4901/88‑1
                                                                                          TTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATG                                                      >  1:162882/1‑47
                                                                                          TTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATG                                                      <  2:162882/47‑1
                                                                                            CTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGC         <  2:553543/90‑1
                                                                                            CTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGC         <  2:656218/90‑1
                                                                                              CCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTA       <  1:696215/90‑1
                                                                                               CCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAT      >  2:55999/1‑90
                                                                                                   ACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAC  <  2:675398/90‑1
                                                                                                                                                                                             
TACTGATTGCTCAATACCACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCCCAAACCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/1203349‑1203246
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAC  >  NZ_CP009273/1205075‑1205159

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