New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 519440 =13 (0.450)4 (0.150) 4/154 1.8 26.0% coding (724/4281 nt) rhsD rhs element protein RhsD
?NZ_CP009273 519532 = 11 (0.420)coding (816/4281 nt) rhsD rhs element protein RhsD

GAGGGCGGGAACTGTCAGAAGAAGATAGCGAAGAGGTGCGGGCCTCTTCCGCACGCTGCGCCTGCGTGGTCAGCACCAGACGGAACTCCCGCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/519533‑519440
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcccgcccTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCC  >  NZ_CP009273/519532‑519581
                                                                                                                                                
GAGGGCGGGAACTGTCAGAAGAAGATAGCGAAGAGGTGCGGGCCTCTTCCGCACGCTGCGCCTGCGTGGTCAGCACCAGACGGAACTCCC                                                        >  1:302164/1‑90
  GGGCGGGAACTGTCAGAAGAAGATAGCGAAGAGGTGCGGGCCTCTTCCGCACGCTGCGCCTGCGTGGTCAGCACCAGACGGAACTCCCGC                                                      >  1:741004/1‑90
  GGGCGGGAACTGTCAGAAGAAGATAGCGAAGAGGTGCGGGCCTCTTCCGCACGCTGCGCCTGCGTGGTCAGCACCAGACGGAACTCC                                                         <  1:832111/87‑1
  GGGCGGGAACTGTCAGAAGAAGATAGCGAAGAGGTGCGGGCCTCTTCCGCACGCTGCGCCTGCGTGGTCAGCACCAGACGGAACTCC                                                         >  2:832111/1‑87
  GGGCGGGAACTGTCAGAAGAAGATAGCGAAGAGGTGCGGGCCTCTTCCGCACGCTGCGCCTGCGTGGTCAGCACCAGACGAAACTCCCGC                                                      >  2:2220/1‑90
    GCGGGAACTGTCAGAAGAAGATAGCGAAGAGGTGCGGGCCTCTTCCGCACGCTGCGCCTGCGTGGTCAGCACCAGACGGAACTCCCGCCC                                                    <  2:657187/90‑1
                                                  GCACGCTGCGCCTGCGTGGTCAGCACCAGACGGAACTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATAC      >  2:702204/1‑90
                                                      GCTGCGCCTGCGTGGTCAGCACCAGACGGAACTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCC  <  1:702204/90‑1
                                                         GCGCCTGCGTGGTCAGCACCAGACGGAACTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCC                  >  1:534888/1‑71
                                                         GCGCCTGCGTGGTCAGCACCAGACGGAACTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCC                  <  2:534888/71‑1
                                                                                                                                                
GAGGGCGGGAACTGTCAGAAGAAGATAGCGAAGAGGTGCGGGCCTCTTCCGCACGCTGCGCCTGCGTGGTCAGCACCAGACGGAACTCCCGCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/519533‑519440
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcccgcccTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCC  >  NZ_CP009273/519532‑519581

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