New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 1391494NA (NA)33 (0.710) 19/158 0.8 91.8% intergenic (+145/‑128) BW25113_RS06980/ynaJ IS5‑like element IS5 family transposase/DUF2534 family protein
?NZ_CP009273 1972747 = 3 (0.060)intergenic (‑293/‑487) flhD/uspC flagellar transcriptional regulator FlhD/universal stress protein UspC

AGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACG  >  NZ_CP009273/1391417‑1391582
                                                                             |                                                                                        
agaAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCa                                                                              <  2:207725/90‑1 (MQ=38)
      aaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGt                                                                        <  1:981116/90‑1 (MQ=255)
        tCGGCTCAACCATGAAGAAATGAAATCACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCTCCTTCCCTAACTAATCAATGCGTTg                                                                      <  1:1218561/90‑1 (MQ=17)
        tCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTg                                                                      <  2:502898/90‑1 (MQ=255)
                   aTGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATcc                                                           <  1:1247267/90‑1 (MQ=255)
                   aTGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATcc                                                           <  1:450622/90‑1 (MQ=255)
                        gaaatgaaatGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTa                                                      <  1:1428163/90‑1 (MQ=255)
                           atgaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAaga                                                   >  2:1483453/1‑90 (MQ=255)
                                 tgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTATCAATGCGTTGATTGTAAATCCAGCTAAgaggtgagg                                            <  1:48683/90‑1 (MQ=255)
                                  gaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAgaggtgagg                                            <  2:577854/90‑1 (MQ=255)
                                        gTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCa                                      <  2:579920/90‑1 (MQ=255)
                                            gCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGc                                  >  2:436739/1‑90 (MQ=255)
                                              cGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAg                                <  2:605607/90‑1 (MQ=255)
                                                agaagaATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGAc                              <  2:1110777/90‑1 (MQ=255)
                                                agaagaATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGAc                              <  1:401754/90‑1 (MQ=255)
                                                  aagaaTTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACaa                            <  1:1483453/90‑1 (MQ=255)
                                                         tCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGaa                     >  1:481147/1‑90 (MQ=255)
                                                           cccGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAAt                   <  2:414148/90‑1 (MQ=255)
                                                                     cGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTAt         <  2:1178213/90‑1 (MQ=255)
                                                                           ttCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTAc   >  2:956190/1‑90 (MQ=255)
                                                                            tCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGTTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACg  >  2:706134/1‑90 (MQ=255)
                                                                             |                                                                                        
AGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACG  >  NZ_CP009273/1391417‑1391582

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