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

TGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC  >  NZ_CP009273/360022‑360193
                                                                                  |                                                                                         
tGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTg                                                                                    >  2:305604/1‑90 (MQ=255)
tGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTg                                                                                    >  2:85018/1‑90 (MQ=255)
  ccAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGc                                                                                  <  1:204182/90‑1 (MQ=255)
                        caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGcc                                                            <  1:83493/90‑1 (MQ=255)
                             tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCa                                                       <  1:143208/90‑1 (MQ=255)
                                          cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATCCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGt                                          <  2:359161/90‑1 (MQ=255)
                                                      gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTCATGGCGGGCGGCCTGCCCCCCACCCCCCCGGGCCTTGCTTCGCAACg                              >  1:136775/1‑90 (MQ=255)
                                                                   agaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAACGTTCAAATcccccc                 >  2:141741/1‑86 (MQ=17)
                                                                             tCCCCAATTCCTGGCAGTTCATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt       <  2:136775/90‑1 (MQ=255)
                                                                                  aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  1:247777/90‑1 (MQ=35)
                                                                                  aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  1:87466/90‑1 (MQ=35)
                                                                                  aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc  <  2:315464/90‑1 (MQ=35)
                                                                                  |                                                                                         
TGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCC  >  NZ_CP009273/360022‑360193

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