New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 376748 =NA (NA)4 (0.200) 4/152 1.4 NA intergenic (‑5/‑59) BW25113_RS01855/BW25113_RS01870 transferase/IS3‑like element IS2 family transposase
?NZ_CP009273 376768 = NA (NA)intergenic (‑25/‑39) BW25113_RS01855/BW25113_RS01870 transferase/IS3‑like element IS2 family transposase

ATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGA  >  NZ_CP009273/376664‑376802
                                                                                    |                                                      
aTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAc                                                   >  1:188125/1‑90 (MQ=255)
            ttGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGa                                       <  2:280151/90‑1 (MQ=255)
              gCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                     <  2:73082/90‑1 (MQ=255)
               cTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAg                                    <  1:330325/90‑1 (MQ=255)
                   ttattGTGCCAGGACGAATAAACGCATAATTACTGGACTTGCCCCCTGAATCTCCAGACAACCAAAAACACTTAAATAAGTGATAGTCtt                                >  1:317633/1‑90 (MQ=21)
                        gtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                        <  1:584905/61‑1 (MQ=255)
                        gtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                        >  2:584905/1‑61 (MQ=255)
                        gtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATAc                           >  1:151730/1‑90 (MQ=255)
                             aGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                      >  1:394082/1‑90 (MQ=255)
                                  gAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                 >  1:415938/1‑90 (MQ=255)
                                        aCGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGt           >  2:217087/1‑90 (MQ=255)
                                              aaTTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg     <  2:1963/90‑1 (MQ=35)
                                                 tACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgaga  <  2:373561/90‑1 (MQ=31)
                                                                                    |                                                      
ATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGA  >  NZ_CP009273/376664‑376802

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