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

CCGCTTCGGGCGGACGCTGACGTACCGGCGTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGC  >  NZ_CP009273/519358‑519529
                                                                                  |                                                                                         
ccGCTTCGGGCGGACGCTGACGTACCGGCGTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGa                                                                                    >  2:344801/1‑90 (MQ=255)
  gCTTCGGGCGGACGCTGACGTACCGGCGTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGt                                                                                  >  2:581264/1‑90 (MQ=255)
                     gTACCGGCGTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCAc                                                               <  2:7120/90‑1 (MQ=255)
                             gTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGgcgc                                                       <  1:28741/90‑1 (MQ=255)
                               gAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAg                                                     >  2:49599/1‑90 (MQ=255)
                                ggccgtattcggtaccgggcagtgtgtcggggaacgctgaggctgagagaGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGc                                                    <  1:702204‑M1/40‑1 (MQ=255)
                                    gtattcggtaccgggcagtgtgtcggggaacgctgaggctgagagaGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGc                                                >  2:702204‑M1/47‑90 (MQ=255)
                                                gggcagtgtgtcggggaacgctgaggctgagagaGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGgcgc                                                       >  1:534888‑M1/35‑71 (MQ=255)
                                                gggcagtgtgtcggggaacgctgaggctgagagaGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGgcgc                                                       <  2:534888‑M1/37‑1 (MQ=255)
                                                ggCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGGGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCg                                    >  1:282706/1‑90 (MQ=255)
                                                     gggAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACct                               >  1:493435/1‑90 (MQ=255)
                                                      ggAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACctc                              <  2:253276/90‑1 (MQ=255)
                                                        aaaTCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACCTCtt                            <  2:727355/90‑1 (MQ=255)
                                                                        ggcgggTGCCGGGGGGGAGTTCCGTCGGGTGCTGCCCACGCAGGCGCAGCGGGGGGAAGAGGCCCGCCCCTCTTCGCTATCTTCTTCTGa            <  1:608015/86‑1 (MQ=25)
                                                                          acggTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACa          <  2:493435/90‑1 (MQ=255)
                                                                             gTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAag                                           <  1:836751/54‑1 (MQ=255)
                                                                             gTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAag                                           >  2:836751/1‑54 (MQ=255)
                                                                                  GGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGc  <  2:657187‑M1/90‑1 (MQ=255)
                                                                                  |                                                                                         
CCGCTTCGGGCGGACGCTGACGTACCGGCGTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACCTCTTCGCTATCTTCTTCTGACAGTTCCCGC  >  NZ_CP009273/519358‑519529

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