New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 360104 =NA (NA)12 (0.840) 10/170 0.3 42.9% intergenic (‑35/+251) lacZ/BW25113_RS23320 beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein
?NZ_CP009273 = 360286 16 (1.120)intergenic (‑217/+69) lacZ/BW25113_RS23320 beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein

TCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGC  >  NZ_CP009273/360197‑360371
                                                                                         |                                                                                     
tCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                                                                       <  2:24584/90‑1 (MQ=38)
   ggAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaat                                                                                    <  2:47970‑M2/90‑4 (MQ=255)
   ggAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaat                                                                                    <  2:322762‑M2/90‑4 (MQ=255)
        gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctg                                                                               <  1:203110‑M2/90‑9 (MQ=255)
        gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                                                                               >  1:264817/1‑90 (MQ=255)
        gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                                                                               >  1:63347/1‑90 (MQ=255)
          gTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggc                                                                             <  1:272388‑M2/90‑11 (MQ=255)
                cGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGa                                                                       <  1:273041/90‑1 (MQ=255)
                cGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGa                                                                       <  2:161040/90‑1 (MQ=255)
                          caGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAgg                                                             <  2:264817/90‑1 (MQ=255)
                            gATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcc                                                           >  2:173069‑M2/1‑62 (MQ=255)
                            gATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGc                                                           >  1:303222/1‑90 (MQ=255)
                            gATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGc                                                           <  1:322762/90‑1 (MQ=255)
                              tAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCt                                                         <  1:366538/90‑1 (MQ=255)
                                    gAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcctgcccgcc                                                   >  1:122998‑M2/1‑54 (MQ=255)
                                    gAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAt                                                   <  2:240317/90‑1 (MQ=255)
                                        ggCCCAGTCTTTCGACTGAGCTTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGCCGAAAGGGCCTCGTGATACGc                                               <  1:363744/90‑1 (MQ=255)
                                          ccccGTATTTCGCATGAGCCATTCTTTTTATTTGATGCCCCGGGAACTGGGAATTAATACTTGAAGACGAAAGGGCCCCGTACTAAGCCt                                             <  2:269321/90‑1 (MQ=255)
                                                   tCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTtata                                    >  2:113904/1‑90 (MQ=255)
                                                        tGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctgg                                                                              <  1:322269‑M2/43‑10 (MQ=255)
                                                        tGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctgg                                                                              >  2:322269‑M2/1‑34 (MQ=255)
                                                        tGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTa                               >  1:120476/1‑90 (MQ=255)
                                                            ccTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcctgcccgccaccctccgggccgttgcttcgcaa                           >  1:207358‑M2/1‑30 (MQ=255)
                                                                   ttttATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGa                    >  2:321075/1‑90 (MQ=255)
                                                                            aTGCCTGGGGAATTaattcctggcagtttatggcgggcgtcctgcccgccac                                                 >  2:90845‑M2/1‑14 (MQ=255)
                                                                                     gaattgggaattAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGc  <  2:124099/90‑1 (MQ=255)
                                                                                     gaattgggaattAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGc  <  2:26009/90‑1 (MQ=255)
                                                                                     gAATTaattcctggcagtttatggcgggcgtcctgcccgccaccctccgggccgttgcttcgcaacgttcaaatccgctcccggcggatt  >  1:350436‑M2/1‑5 (MQ=255)
                                                                                     gAATTaattcctggcagtttatggcgggcgtcctgcccgccaccctccgggccgttgcttcgcaacgttcaaatccgctcccggcggatt  >  1:263809‑M2/1‑5 (MQ=255)
                                                                                         |                                                                                     
TCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGC  >  NZ_CP009273/360197‑360371

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