New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3756668 =NA (NA)3 (0.060) 3/146 3.9 20.8% coding (1126/4134 nt) rhsA RHS element protein RhsA
?NZ_CP009273 3759111 = 13 (0.240)coding (3569/4134 nt) rhsA RHS element protein RhsA
Rejected: Coverage evenness skew score above cutoff.

GCCCTGCGTGTGCAGCACTTCACGGCGGTCCAGGCTGTCGGTGATGGTGATGCGGTCTTTCTCATACTGATACGTGTAGCTTAAgcctgccggg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3756751‑3756668
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAG  >  NZ_CP009273/3759111‑3759202
                                                                                                                                                                                
GCCCTGCGTGTGCAGCCCTTCACGGCGGTCCAGGCTGTCGGTGATGGTGATGCGGTCTTTCTCATACTGATACGTGTAGCTTAAGCCTGC                                                                                        <  2:901262/90‑1
GCCCTGCGTGTGCAGCACTTCACGGCGGTCCAGGCTGTCGGTGATGGTGATGCGGTCTTTCTCATACTGATACGTGTAGCTTAAGCCTGC                                                                                        <  2:1574035/90‑1
GCCCTGCGTGTGCAGCACTTCACGGCGGTCCAGGCTGTCGGTGATGGTGATGCGGTCTTTCTCATACTGATACGTGTAGCTTAAGCCTGC                                                                                        <  1:663637/90‑1
GCCCTGCGTGTGCAGCACTTCACGGCGGTCCAGGCTGTCGGGGATGGTGATGCGGTCTTTCTCATACTGATTCGTGTAGCTTAAGCCTGC                                                                                        >  2:720710/1‑90
  cagGCGTGTGCAGCACTTCACGGCGGTCCAGGCTGTCGGTGATGGTGATGCGGTCTTTCTCATACTGATACGTGTAGCTTAAGCCTGCCG                                                                                      <  1:1164289/87‑1
     GCGTGTGCAGCACTTCACGGCGGTCCAGGCTGTCGGGGATGGGGGTGCGGGCTTTTTCATACTGGTACGTGTAGCTTAAGCCTGCCGGGG                                                                                   >  1:1083509/1‑90
                                     TCGGTGATGGTGATGCGGTCTTTCTCATACTGATACGTGTAGCTTAAGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCT                                                         <  2:1299457/84‑1
                                     TCGGTGATGGTGATGCGGTCTTTCTCATACTGATACGTGTAGCTTAAGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCT                                                         >  1:1299457/1‑84
                                                                                    GCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTC    <  1:409529/90‑1
                                                                                    GCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTC    >  1:1072569/1‑90
                                                                                    GCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGC                                      <  1:417158/56‑1
                                                                                    GCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGC                                      >  2:417158/1‑56
                                                                                     CCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCA   >  1:87002/1‑90
                                                                                     CCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCA   <  2:88662/90‑1
                                                                                      CTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAG  >  1:819657/1‑90
                                                                                      CTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTAT                                                     >  1:840385/1‑39
                                                                                      CTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTAT                                                     <  2:840385/39‑1
                                                                                       TGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCC                                           >  2:1265223/1‑48
                                                                                       TGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCC                                           <  1:1265223/48‑1
                                                                                                                                                                                
GCCCTGCGTGTGCAGCACTTCACGGCGGTCCAGGCTGTCGGTGATGGTGATGCGGTCTTTCTCATACTGATACGTGTAGCTTAAgcctgccggg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3756751‑3756668
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAG  >  NZ_CP009273/3759111‑3759202

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