New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 727836 =NA (NA)3 (0.080) 3/154 2.8 17.2% coding (2798/4194 nt) rhsC rhs element protein RhsC
?NZ_CP009273 727859 = 16 (0.400)coding (2821/4194 nt) rhsC rhs element protein RhsC
Rejected: Frequency below/above cutoff threshold.

ATCATCCGTGCGGATAACCCCTTCCGGGATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGcgtcacgg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/727921‑727836
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGC  >  NZ_CP009273/727859‑727955
                                                                                                                                                                                       
ATCATCCGTGCGGATAACCCCTTCCGGGATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTC                                                                                               <  1:994495/90‑1
 TCATCCGTGCGGATAACCCCTTCCGGGATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCA                                                                                              >  1:228547/1‑90
   ATCCGTGCGGATAACCCCTTCCGGGATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACG                                                                                            <  1:1060897/90‑1
   ATCCGTGCGGATAACCCCTTCCGGGATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACG                                                                                            <  2:693547/90‑1
   ATCCGTGCGGATAACCCCTTCCGGGATGAGGTCGGTTTTCCCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACG                                                                                            <  1:148062/90‑1
    TCCGTGCGGATAACCCCTTCCGGGATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGG                                                                                           <  1:295804/90‑1
    TCCGTGCGGATAACCCCTTCCGGGATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGG                                                                                           <  2:653424/90‑1
    TCCGTGCGGATAACCCCTTCCGGGATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGG                                                                                           >  2:596602/1‑90
       GTGCGGATAACCCCTTCCGGGATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGG                                                                                           >  2:855405/1‑87
       GTGCGGATAACCCCTTCCGGGATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGG                                                                                           <  1:855405/87‑1
          CGGATAACCCCTTCCGGGATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGGCAGGCT                                                                                     <  1:955049/90‑1
                           GATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGCGTCACGGCAGGCTGACAGAGAAAACCGACC                                                                    <  1:393493/90‑1
                                     GTTTTCCGTGTCCGCCTGCCGGGACGGTCATCCCGGTAAAGATAGGGCGCGTCACGGCAGGCTGACAGAGAAACCCGACCTCATCCCGGA                                                          <  1:494519/90‑1
                                                                                        TCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAG       >  2:594915/1‑90
                                                                                         CACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGT      >  2:1063443/1‑90
                                                                                          ACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTC     <  1:286463/90‑1
                                                                                             GCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGAGGAGCGCACCCACCGGTACCATTACGACCGTCAGC  >  1:284240/1‑90
                                                                                                                                                                                       
ATCATCCGTGCGGATAACCCCTTCCGGGATGAGGTCGGTTTTCTCTGTCAGCCTGCCGTGACGGTCATACCGGTAAAGATAGTGCGcgtcacgg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/727921‑727836
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGC  >  NZ_CP009273/727859‑727955

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