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

TCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGCTTT  >  NZ_CP009273/360197‑360374
                                                                                         |                                                                                        
tCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                                                                          >  2:194802/1‑90 (MQ=38)
   ggAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggg                                                                                       >  2:64196/1‑90 (MQ=255)
   ggAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggg                                                                                       >  1:103372/1‑90 (MQ=255)
      gagCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaat                                                                                    >  1:124306/1‑90 (MQ=255)
      gagCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcc                                                                                    <  2:275907‑M2/90‑7 (MQ=255)
        gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctg                                                                                  <  2:276785‑M2/90‑9 (MQ=255)
        gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                                                                                  >  2:173687/1‑90 (MQ=255)
        gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                                                                                  >  1:27318/1‑90 (MQ=255)
        gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta                                                                                  >  1:230821/1‑90 (MQ=255)
             cACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagt                                                                             <  2:304356‑M2/90‑14 (MQ=255)
                cGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGa                                                                          <  1:196566/90‑1 (MQ=255)
                       caacaGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcggg                                                                   >  1:206519‑M2/1‑67 (MQ=255)
                         acaGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAg                                                                 >  2:295237/1‑90 (MQ=255)
                                 aaCGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcctgccc                                                         <  1:204005‑M2/90‑34 (MQ=255)
                                    gAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcctgcccgcc                                                      >  1:138661‑M2/1‑54 (MQ=255)
                                     aaaGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcctgcccgcca                                                     >  2:204005‑M2/1‑53 (MQ=255)
                                           ccAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcctgcccgccaccctcc                                               >  1:221057‑M2/1‑47 (MQ=255)
                                                        tGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggggtcct                                                             <  1:19111‑M2/63‑30 (MQ=255)
                                                        tGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcct                                                             >  2:19111‑M2/1‑34 (MQ=255)
                                                         gAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTaa                                 <  2:87949/90‑1 (MQ=255)
                                                            ccTTTCGTTTTATTTGATGCCTGGGGAATTaattcctggcagtttatggcgggcgtcctgcccgccaccctccgggccgttgcttcgcaa                              >  1:275907‑M2/1‑30 (MQ=255)
                                                                    tttATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGaa                      <  2:35370/90‑1 (MQ=255)
                                                                                     gaattgggaattAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGc     <  2:27318/90‑1 (MQ=255)
                                                                                     gaattgggaattAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGc     <  2:230821/90‑1 (MQ=255)
                                                                                     gAATTaattcctggcagtttatggcgggcgtcctgcccgccaccctccgggccgttgcttcgcaacgttcaaatccgctcccggcggatt     >  1:90630‑M2/1‑5 (MQ=255)
                                                                                        ttaattcctggcagtttatggcgggcgtcctgcccgccaccctccgggccgttgcttcgcaacgttcaaatccgctcccggcggatttgt  >  1:284135‑M2/1‑2 (MQ=255)
                                                                                         |                                                                                        
TCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTATTTTTATAGGTTAATGTCATGGGAAGATCCAGCCGACCAGGCTTT  >  NZ_CP009273/360197‑360374

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