New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 72783557 (0.480)12 (0.110) 10/154 3.5 18.9% coding (2797/4194 nt) rhsC rhs element protein RhsC
?NZ_CP009273 = 729729 NA (NA)pseudogene (126/1506 nt) BW25113_RS03635 RHS repeat protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/727743‑727835
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccgtgacgGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCGGGTCGG  <  NZ_CP009273/729729‑729641
                                                                                                                                                                                      
CCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTG                                                                                              >  2:3212542/1‑90
 CACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGA                                                                                             >  2:2611315/1‑90
  ACAGACCCGGCAGGTAACCGCCTGCCCGCCCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGCGGCCGGAGAACCGTAGCGCCCGTGgc                                                                                            >  2:1256852/1‑88
  ACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAC                                                                                            >  2:688554/1‑90
  ACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAC                                                                                            >  2:3345291/1‑90
  ACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAC                                                                                            >  2:3270545/1‑90
  ACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAC                                                                                            <  1:511179/90‑1
  ACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGAC                                                                                            >  2:2632835/1‑90
   CAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACG                                                                                           >  1:1240930/1‑90
   CAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACG                                                                                           <  2:2648879/90‑1
                                             CGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGGTCATAC                                                                                    <  1:1174543/55‑1
                                             CGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGGTCATAC                                                                                    >  2:1174543/1‑55
                                                               TGTGGCCGGATAACCGTATCGCCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAG                               >  2:3152332/1‑90
                                                               TGTGGCCGGATAACCGTATCGCCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAG                               >  1:2253067/1‑90
                                                                GTGGCCGGATAACCGTATCGCCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGG                              >  2:198279/1‑90
                                                                 TGGCCGGATAACCGTATCGCCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGG                             >  1:1862624/1‑90
                                                                  GGCCGGATAACCGTATCGCCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGT                            >  2:1801443/1‑90
                                                                      GGATAACCGTATCGCCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTG                        >  1:1762564/1‑90
                                                                        ATAACCGTATCGCCCGTGACGGTCATACCGGTAAAGATAGGGCGCGTCACGGGCGATCCGGTTATCCGGCCACATGCTGAGGGGGCTGTC                      >  2:970064/1‑90
                                                                                 TCGCCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAG             >  1:485046/1‑90
                                                                                 TCGCCCGTGACGCTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAG             >  2:1789211/1‑90
                                                                                   GCCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCT           >  1:2992251/1‑90
                                                                                      CGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCG        >  2:629544/1‑90
                                                                                      CGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGG                                                            <  2:2463461/38‑1
                                                                                      CGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGG                                                            >  1:2463461/1‑38
                                                                                         GACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCGGGT     >  2:140591/1‑90
                                                                                         GACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCGGGT     >  1:340482/1‑90
                                                                                         GACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCGGGT     <  2:1145750/90‑1
                                                                                         GACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCGGGT     <  2:2148300/90‑1
                                                                                         GACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCGGGT     >  2:214263/1‑90
                                                                                          ACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCGGGTC    <  2:2869408/90‑1
                                                                                          ACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCGGGTC    <  1:1811228/90‑1
                                                                                          ACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCGGGTC    <  2:361951/90‑1
                                                                                          ACGGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCGGGTC    <  2:2054725/90‑1
                                                                                            GGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCGGGTCGG  >  2:1030397/1‑90
                                                                                                                                                                                      
CCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/727743‑727835
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccgtgacgGTCATACCGGTAAAGATAGTGCGCGTCACGGGCGATACGGTTATCCGGCCACATGCTGAGGGTGCTGTCCGGGTGCAGCTCCGGGTCGG  <  NZ_CP009273/729729‑729641

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