New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 730478 =11 (0.210)3 (0.060) 3/146 3.9 23.7% pseudogene (875/1506 nt) BW25113_RS03635 RHS repeat protein
?NZ_CP009273 3613677 = NA (NA)coding (1126/4236 nt) rhsB rhs element protein RhsB
Rejected: Coverage evenness skew score above cutoff.

ATCCTGAGTGATATACCTCCCCTGCAGCGGGTCATAATAGCGGTGGCGGTTGTAATACAGGCCGGACTCCTCATCATACTGCTGCCCCGGCAGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/730572‑730478
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cccggcaggcTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGC  >  NZ_CP009273/3613677‑3613760
                                                                                                                                                                                   
ATCCTGAGTGATATACCTCCCCTGCAGCGGGTCATAATAGCGGTGGCGGTTGTAATACAGGCCGGACTCCTCATCATACTGCTGCCCCGG                                                                                           >  1:1482600/1‑90
   CTGAGTGATATACCTCCCCTGCAGCGGGTCATAATAGCGGTGGCGGTTGTAATACAGGCCGGACTCCTCATCATACTGATGCCCCGGCAG                                                                                        <  2:460292/90‑1
                                       GCGGTGGCGGTTGTAATACAGGCCGGACTCCTCATCATACTGCTGCCCCGGCAGGC                                                                                      >  1:417158/1‑56
                                       GCGGTGGCGGTTGTAATACAGGCCGGACTCCTCATCATACTGCTGCCCCGGCAGGC                                                                                      <  2:417158/56‑1
                                            GGCGGTTGTAATACAGGCCGGACTCCTCATCATACTGCTGCCCCGGCA                                                                                         >  1:1265223/1‑48
                                            GGCGGTTGTAATACAGGCCGGACTCCTCATCATACTGCTGCCCCGGCA                                                                                         <  2:1265223/48‑1
                                                      ATACAGGCCGGACTCCTCATCATACTGCTGCCCCGGCAG                                                                                        <  1:840385/39‑1
                                                      ATACAGGCCGGACTCCTCATCATACTGCTGCCCCGGCAG                                                                                        >  2:840385/1‑39
                                                          AGGCCGGACTCCTCATCATACTGCTGCCCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGA                                       <  1:1299457/84‑1
                                                          AGGCCGGACTCCTCATCATACTGCTGCCCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGA                                       >  2:1299457/1‑84
                                                                                    CCCCGGCAGGCTTAAGCTACACGTACCAGTATGAAAAAGCCCGCACCCCCATCCCCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGC       <  1:1083509/90‑1
                                                                                       CGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCctg    >  1:1164289/1‑87
                                                                                         GCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGC  >  1:663637/1‑90
                                                                                         GCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGC  >  2:1574035/1‑90
                                                                                         GCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGGGCTGCACACGCAGGGC  >  2:901262/1‑90
                                                                                         GCAGGCTTAAGCTACACGAATCAGTATGAGAAAGACCGCATCACCATCCCCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGC  <  2:720710/90‑1
                                                                                                                                                                                   
ATCCTGAGTGATATACCTCCCCTGCAGCGGGTCATAATAGCGGTGGCGGTTGTAATACAGGCCGGACTCCTCATCATACTGCTGCCCCGGCAGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/730572‑730478
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cccggcaggcTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGC  >  NZ_CP009273/3613677‑3613760

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