New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 1391494NA (NA)14 (0.670) 10/170 0.6 93.4% intergenic (+145/‑128) BW25113_RS06980/ynaJ IS5‑like element IS5 family transposase/DUF2534 family protein
?NZ_CP009273 1972747 = 1 (0.050)intergenic (‑293/‑487) flhD/uspC flagellar transcriptional regulator FlhD/universal stress protein UspC

GAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCAT  >  NZ_CP009273/1391416‑1391570
                                                                              |                                                                            
gagaAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATc                                                                   >  2:101912/1‑90 (MQ=37)
         tCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTg                                                          <  1:119703/90‑1 (MQ=255)
                           aatgaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAag                                        >  2:285074/1‑90 (MQ=255)
                                   gaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAgaggtgagg                                <  2:135287/90‑1 (MQ=255)
                                       gAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGtttt                            >  2:456292/1‑90 (MQ=255)
                                             gCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTg                                                          <  1:111869/54‑1 (MQ=255)
                                             gCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTg                                                          >  2:111869/1‑54 (MQ=255)
                                                 agaagaATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGAc                  >  1:386264/1‑90 (MQ=255)
                                                  gaagaaTTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACa                 >  2:268209/1‑90 (MQ=255)
                                                    agaaTTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAAc               <  2:386264/90‑1 (MQ=255)
                                                         ttCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGa          <  2:350128/90‑1 (MQ=255)
                                                          tCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGaa         <  2:34479/90‑1 (MQ=255)
                                                           ccccGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACATCGGTGaaa        <  2:248019/90‑1 (MQ=255)
                                                                 ttATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCAt  <  1:69349/90‑1 (MQ=255)
                                                                              |                                                                            
GAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCAT  >  NZ_CP009273/1391416‑1391570

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