New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 360104 =NA (NA)14 (0.730) 10/166 0.6 48.3% intergenic (‑35/+251) lacZ/BW25113_RS23320 beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein
?NZ_CP009273 = 360286 15 (0.780)intergenic (‑217/+69) lacZ/BW25113_RS23320 beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein

GGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGCTTT  >  NZ_CP009273/360200‑360374
                                                                                      |                                                                                        
ggAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggg                                                                                       >  1:165916/1‑90 (MQ=255)
ggAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGACCCTTTCGTTTTATTTGATGCCTGgggaattggg                                                                                       >  1:360345/1‑90 (MQ=255)
     gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                                                                          >  1:200735/1‑82 (MQ=38)
     gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                                                                          <  2:200735/82‑1 (MQ=38)
     gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctg                                                                                  <  2:311678‑M2/90‑9 (MQ=255)
     gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctg                                                                                  <  2:369048‑M2/90‑9 (MQ=255)
     gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                                                                                  >  2:352893/1‑90 (MQ=255)
       gTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggc                                                                                >  1:511592‑M2/1‑80 (MQ=255)
          cACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagt                                                                             <  1:263258‑M2/90‑14 (MQ=255)
          cACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagt                                                                             <  2:507094‑M2/90‑14 (MQ=255)
                   acaacaGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgg                                                                    >  1:258207‑M2/1‑68 (MQ=255)
                      acaGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAg                                                                 <  1:487550/90‑1 (MQ=255)
                          aTAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcct                                                             >  2:263258‑M2/1‑61 (MQ=255)
                                        ccAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTa                                               >  2:309457/1‑90 (MQ=255)
                                               ttCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTtat                                        <  2:478002/90‑1 (MQ=255)
                                                     tGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggggtgctgccccccccccccccggtcgttgtttt                                  <  1:510728‑M2/90‑57 (MQ=255)
                                                     tGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTa                                  >  2:523877/1‑90 (MQ=255)
                                                     tGAGCCTTTCGTTTTATGTGATGCCTGGGGAATTaattcctggcagtttatggcgggggtcctgccccccccccccccggccgttgcttt                                  <  2:12294‑M2/90‑57 (MQ=255)
                                                         ccTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcctgcccgccaccctccgggccgttgcttcgcaa                              >  2:58185‑M2/1‑30 (MQ=255)
                                                         ccTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcctgcccgccaccctccgggccgttgcttcgcaa                              >  2:239097‑M2/1‑30 (MQ=255)
                                                                ttttATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGa                       >  2:275850/1‑90 (MQ=255)
                                                                ttttATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGa                       <  1:507094/90‑1 (MQ=255)
                                                                         aTGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGc              >  1:517692/1‑90 (MQ=255)
                                                                             cTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGAc          >  2:332241/1‑90 (MQ=255)
                                                                              tGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGAcc         >  2:295544/1‑90 (MQ=255)
                                                                               ggggAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCa        <  1:369048/90‑1 (MQ=255)
                                                                                  gaattgggaattAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGc     <  1:265778/90‑1 (MQ=255)
                                                                                    attgggaattAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGCtt   >  1:268554/1‑90 (MQ=255)
                                                                                    attaattcctggcagtttatggcgggcgtcctgcccgccac                                                    >  1:132252‑M2/1‑3 (MQ=255)
                                                                                    attaattcctggcagtttatggcgggcgtcctgccccccac                                                    <  2:132252‑M2/41‑39 (MQ=255)
                                                                                     ttaattcctggcagtttatggcgggcgtcctgcccgccaccctccgggccgttgcttcgcaacgttcaaatccgctcccggcggatttgt  >  2:442374‑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 < 20 ≤ ATCG/ATCG < 33 ≤ 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.