New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3614374 =10 (0.370)5 (0.200) 5/154 1.6 23.2% coding (1823/4236 nt) rhsB rhs element protein RhsB
?NZ_CP009273 3614426 = 24 (0.970)coding (1875/4236 nt) rhsB rhs element protein RhsB

CAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTAT  >  NZ_CP009273/3614289‑3614455
                                                                                     |                                                                                 
gtggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgtcagcccgccctgcgtggtacggcc                                                                               >  1:655971‑M1/86‑88 (MQ=255)
            ttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgtcagcccgccctgcgtggtaCGGACGGCAGCAGAAAc                                                                   >  1:338010‑M1/74‑90 (MQ=255)
              cactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgtcagcccgccctgcgtggtaCGGACGGCAGCAGAAACgg                                                                 >  1:119370‑M1/72‑90 (MQ=255)
                                                          ctgcgcgtcagcccgccctgcgtggtaCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAgg                     <  2:76866‑M1/63‑1 (MQ=255)
                                                                    gcccgccctgcgtggtaCGGACGGCAGCAGAAACGGGACACAGTACGATGCGGGGGGAAAGGCCGTCCGTGCCACGCAGGGCGGGCTGAc           <  1:212804‑M1/73‑1 (MQ=255)
                                                                        gccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACgcgc       <  1:122101/90‑1 (MQ=35)
                                                                        gccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACgcgc       <  1:350677/90‑1 (MQ=33)
                                                                        gccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACgcgc       <  1:439275/90‑1 (MQ=35)
                                                                        gccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACgcgc       <  2:579868/90‑1 (MQ=35)
                                                                        gccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACgcgc       <  2:586358/90‑1 (MQ=35)
                                                                         ccctgcgtggtaCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGGCCGGACCACGCAGGGCGGGCTGACGCGCa      <  1:384583‑M1/78‑1 (MQ=255)
                                                                          cgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAg     <  1:728884/90‑1 (MQ=35)
                                                                          cgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAg     <  2:152260/90‑1 (MQ=35)
                                                                             cATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTAt  <  1:215885/90‑1 (MQ=35)
                                                                                  ccccGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCa      >  1:527003/1‑81 (MQ=33)
                                                                                  ccccGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCa      <  2:527003/81‑1 (MQ=35)
                                                                                     |                                                                                 
CAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTAT  >  NZ_CP009273/3614289‑3614455

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