New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 360104 =NA (NA)29 (1.230) 21/172 0.1 NA intergenic (‑35/+251) lacZ/BW25113_RS23320 beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein
?NZ_CP009273 = 360286 NA (NA)intergenic (‑217/+69) lacZ/BW25113_RS23320 beta‑galactosidase/relaxase/mobilization nuclease domain‑containing protein

TTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC  >  NZ_CP009273/360032‑360193
                                                                        |                                                                                         
ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATg                                                                          >  1:569953/1‑90 (MQ=255)
       ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggcgggcg                                                                   >  2:538473/1‑90 (MQ=255)
        cAGCACCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGt                                                                  >  1:342235/1‑90 (MQ=255)
               accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGgcccgccc                                                           >  2:361648/1‑90 (MQ=255)
                   tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCa                                                       <  1:191341/90‑1 (MQ=255)
                   tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCa                                                       <  2:22706/90‑1 (MQ=255)
                       gTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggcggg                                                                     <  1:162884/72‑1 (MQ=255)
                       gTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATggcggg                                                                     >  2:162884/1‑72 (MQ=255)
                       gTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCCGccccccccccc                                                   >  2:513066/1‑89 (MQ=255)
                       gTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGccccccccccccccc                                                   >  1:606582/1‑89 (MQ=255)
                         gCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTcc                                                 <  2:569953/90‑1 (MQ=255)
                            ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGCccccccccccccccg                                              >  2:166326/1‑87 (MQ=38)
                              aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGcc                                            >  2:75818/1‑90 (MQ=37)
                                cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGt                                          <  2:342235/90‑1 (MQ=255)
                                cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGt                                          <  1:248857/90‑1 (MQ=255)
                                 tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCccgccgct                                         >  2:521019/1‑88 (MQ=37)
                                          aaGCTGTCAACCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                <  1:158755/90‑1 (MQ=255)
                                          aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                <  2:584128/90‑1 (MQ=255)
                                          aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                <  2:318799/90‑1 (MQ=255)
                                          aaGCTGCCAAACCTGAGAATTCCGCCCCCCAATCCCTGGCAGTTTATGGCGGGCGCCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCaa                                <  1:367508/90‑1 (MQ=17)
                                              tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCCCCCCCCCCCCCCCCCCGGCGGTGCTTTCCAACGtt                            >  1:51654/1‑90 (MQ=255)
                                                              tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc            <  2:201885/90‑1 (MQ=255)
                                                              tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc            <  1:143571/90‑1 (MQ=255)
                                                               ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCCGCCGCTGCTTCGCAACATTCAAAAcccccccccccg           >  1:98368/1‑86 (MQ=9)
                                                                   tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt       <  2:2284/90‑1 (MQ=255)
                                                                      ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt    <  2:641033/90‑1 (MQ=38)
                                                                        aaTTCCTGGCAGTTTTTGGCGGGCGTCCTTCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  2:374826/90‑1 (MQ=255)
                                                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  2:89526/90‑1 (MQ=25)
                                                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  2:404689/90‑1 (MQ=35)
                                                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  2:192701/90‑1 (MQ=35)
                                                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  2:606582/90‑1 (MQ=35)
                                                                        aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  1:222719/90‑1 (MQ=35)
                                                                        aaTTCCTGGCAGTTTATGGCGGCCGCCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  1:546087/90‑1 (MQ=11)
                                                                        aaTTCCGGGCAGTTTATGGCGGGCGTCCGGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  1:608572/90‑1 (MQ=11)
                                                                        |                                                                                         
TTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC  >  NZ_CP009273/360032‑360193

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