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

GGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCG  >  NZ_CP009273/519444‑519607
                                                                                        |                                                                           
GGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCctc                                                                            >  1:302164‑M1/1‑90 (MQ=255)
                                cAGCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGccc                                            <  1:177162/90‑1 (MQ=255)
                                  gCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCgg                                          <  1:756393/90‑1 (MQ=255)
                                  gCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCgg                                          >  2:756393/1‑90 (MQ=255)
                                    gTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGacacac                                                 >  1:203023/1‑81 (MQ=255)
                                    gTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGacacac                                                 <  2:203023/81‑1 (MQ=255)
                                            gcacgctgcgcctgcgtggtcagcaccagacggaactcccgcccTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATAc                                >  2:702204‑M2/45‑90 (MQ=255)
                                            gagGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATAc                                <  2:263890/90‑1 (MQ=255)
                                              ggCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGacacac                                                 >  1:825726/1‑71 (MQ=255)
                                              ggCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGacacac                                                 <  2:825726/71‑1 (MQ=255)
                                                gctgcgcctgcgtggtcagcaccagacggaactcccgcccTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGcc                            <  1:702204‑M2/50‑1 (MQ=255)
                                                   gcgcctgcgtggtcagcaccagacggaactcccgcccTCTCTCAGCCTCAGCGTTCCCCGACACACTGccc                                            >  1:534888‑M2/38‑71 (MQ=255)
                                                   gcgcctgcgtggtcagcaccagacggaactcccgcccTCTCTCAGCCTCAGCGTTCCCCGACACACTGccc                                            <  2:534888‑M2/34‑1 (MQ=255)
                                                   gCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCg                         >  1:505373/1‑90 (MQ=255)
                                                       ctctTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAg                     <  1:657187/90‑1 (MQ=255)
                                                                          gACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTcggcg  >  2:531435/1‑90 (MQ=255)
                                                                                        |                                                                           
GGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGCCCTCTCTCAGCCTCAGCGTTCCCCGACACACTGCCCGGTACCGAATACGGCCCCGACAGGGGTATCCGCCTTTCGGCG  >  NZ_CP009273/519444‑519607

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