New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 360104 =NA (NA)9 (0.650) 7/164 0.6 47.4% intergenic (‑35/+251) lacZ/BW25113_RS23320 beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein
?NZ_CP009273 = 360286 10 (0.720)intergenic (‑217/+69) lacZ/BW25113_RS23320 beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein

GGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGCTTT  >  NZ_CP009273/360200‑360374
                                                                                      |                                                                                        
ggAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaat                                                                                       <  1:30582‑M2/90‑4 (MQ=255)
ggAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaat                                                                                       <  2:261947‑M2/90‑4 (MQ=255)
     gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctg                                                                                  <  1:249497‑M2/90‑9 (MQ=255)
     gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctg                                                                                  <  1:50596‑M2/90‑9 (MQ=255)
        ttCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggca                                                                               <  2:218219‑M2/90‑12 (MQ=255)
                  aacaacAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACg                                                                     >  2:44541/1‑90 (MQ=255)
                      acaGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAg                                                                 <  1:165344/90‑1 (MQ=255)
                         gATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGc                                                              <  1:218219/90‑1 (MQ=255)
                              aaCGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGt                                                         <  1:333385/90‑1 (MQ=255)
                                  aaaGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcctgcccgcca                                                     >  1:118811‑M2/1‑53 (MQ=255)
                                               ttCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTta                                         >  1:386596/1‑89 (MQ=255)
                                               ttCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTta                                         <  2:386596/89‑1 (MQ=255)
                                                     tGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcggggggcctgccccccccccccccggccgttgcttt                                  <  2:76127‑M2/90‑57 (MQ=255)
                                                            ttCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCAt                           <  1:44541/90‑1 (MQ=255)
                                                                ttttATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGa                       >  2:56342/1‑90 (MQ=255)
                                                                             cTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGAc          >  1:261201/1‑90 (MQ=255)
                                                                                  gaattgggaattAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGc     <  1:277622/90‑1 (MQ=255)
                                                                                  gAATTaattcctggcagtttatggcgggcgtcctgcccgccaccctccgggccgttgcttcgcaacgttcaaatccgctcccggcggatt     >  2:254078‑M2/1‑5 (MQ=255)
                                                                                     ttaattcctggcagtttatggcgggcgtcctgcccgccaccctccgggccgttgcttcgcaacgttcaaatccgctcccggcggatttgt  >  2:95725‑M2/1‑2 (MQ=255)
                                                                                      |                                                                                        
GGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGCTTT  >  NZ_CP009273/360200‑360374

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