New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727511 (0.020)29 (0.720) 18/170 0.6 96.7% intergenic (‑297/‑483) flhD/uspC flagellar transcriptional regulator FlhD/universal stress protein UspC
?NZ_CP009273 = 2283592 NA (NA)intergenic (‑68/‑387) BW25113_RS11435/narP IS5‑like element IS5 family transposase/nitrate/nitrite response regulator protein NarP

AGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAA  >  NZ_CP009273/2283506‑2283678
                                                                                      |                                                                                      
aGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATa                                                                                     >  2:869713/1‑90 (MQ=33)
aGGGAAGTTGATGACTCATGATGACCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCCTATCAGGGACTTGTTCCCACCTTCCATa                                                                                     <  1:290436/90‑1 (MQ=11)
 ggTAAGTTGATGACTCATGATGAACCCTGTTCTATGCCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATaa                                                                                    <  1:880671/90‑1 (MQ=18)
    aaGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGc                                                                                 <  2:1007725/90‑1 (MQ=39)
      gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTg                                                                               >  2:863767/1‑90 (MQ=40)
         gatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                                            <  1:863767/90‑1 (MQ=255)
         gatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                                            <  1:733965/90‑1 (MQ=255)
         gatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                                            >  2:612515/1‑90 (MQ=255)
           tgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGcc                                                                          >  2:645126/1‑90 (MQ=255)
              cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacc                                                                       <  2:82753/90‑1 (MQ=255)
               tCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacca                                                                      <  2:567607/90‑1 (MQ=255)
               tCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacca                                                                      <  2:488543/90‑1 (MQ=255)
                        aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGata                                                             >  1:590218/1‑90 (MQ=255)
                        aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGata                                                             >  2:759163/1‑90 (MQ=255)
                           cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTg                                                          >  2:46178/1‑90 (MQ=255)
                             gTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCg                                                        >  1:151505/1‑90 (MQ=255)
                                  aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACg                                                   >  2:622160/1‑90 (MQ=255)
                                  aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACg                                                   >  1:85139/1‑90 (MQ=255)
                                  aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACg                                                   >  1:869263/1‑90 (MQ=255)
                                       tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacc                                                                       <  1:522886/65‑1 (MQ=255)
                                       tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacc                                                                       >  2:522886/1‑65 (MQ=255)
                                       tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaaa                                              >  2:71697/1‑90 (MQ=255)
                                       tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaaa                                              >  2:9472/1‑90 (MQ=255)
                                       tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaaa                                              >  1:782099/1‑90 (MQ=255)
                                       tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaaa                                              >  1:116998/1‑90 (MQ=255)
                                         cAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAaga                                            <  2:210409/90‑1 (MQ=255)
                                          aGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAagag                                           <  1:1027268/90‑1 (MQ=255)
                                            aTGACAAACATGATCTCATATCAGGGACTTGTACGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaa                                         <  2:659946/90‑1 (MQ=255)
                                              gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAAc                                       <  1:851981/90‑1 (MQ=255)
                                              gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAAc                                       <  2:151505/90‑1 (MQ=255)
                                              gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAAc                                       <  2:271637/90‑1 (MQ=255)
                                               aCAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGaa                                                                   >  2:928987/1‑61 (MQ=255)
                                               aCAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGaa                                                                   <  1:928987/61‑1 (MQ=255)
                                                       gATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGa                              >  1:982415/1‑90 (MQ=255)
                                                           tcATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGa                          <  1:303237/90‑1 (MQ=255)
                                                           tcATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGa                          <  2:782099/90‑1 (MQ=255)
                                                           tcATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGa                          <  1:71697/90‑1 (MQ=255)
                                                             atatCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATg                        >  2:677016/1‑90 (MQ=255)
                                                                  aGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTg                   <  2:85139/90‑1 (MQ=255)
                                                                   gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaa                                         >  2:613231/1‑67 (MQ=255)
                                                                   gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaa                                         <  1:613231/67‑1 (MQ=255)
                                                                   gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGa                  >  2:1020537/1‑90 (MQ=255)
                                                                   gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGa                  <  2:889861/90‑1 (MQ=255)
                                                                   gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGa                  >  1:112534/1‑90 (MQ=255)
                                                                   gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGa                  >  1:214400/1‑90 (MQ=255)
                                                                   gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGa                  >  1:487563/1‑90 (MQ=255)
                                                                      aCTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCt               >  1:176767/1‑90 (MQ=255)
                                                                      aCTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCt               >  2:269802/1‑90 (MQ=255)
                                                                            tCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATc         >  1:661917/1‑90 (MQ=255)
                                                                             cGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATcc        >  1:273458/1‑90 (MQ=255)
                                                                             cGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATcc        >  1:395696/1‑90 (MQ=255)
                                                                                  cTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGa   <  1:612515/90‑1 (MQ=255)
                                                                                   ttCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGaa  >  1:162941/1‑90 (MQ=255)
                                                                                   ttCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGaa  >  2:234688/1‑90 (MQ=255)
                                                                                      |                                                                                      
AGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAA  >  NZ_CP009273/2283506‑2283678

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