New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727510 (0.000)35 (1.520) 17/164 0.3 100% 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

GTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTG  >  NZ_CP009273/2283508‑2283661
                                                                                    |                                                                     
gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAAc                                                                  <  2:551018/90‑1 (MQ=37)
gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAAc                                                                  <  2:433506/90‑1 (MQ=37)
            cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacc                                                      <  1:42704/90‑1 (MQ=255)
            cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacc                                                      <  1:650800/90‑1 (MQ=255)
                      aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGata                                            >  1:692882/1‑90 (MQ=255)
                        ccTGTTCTATGGCTCCAGATGACAAAAATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                           >  1:218313/1‑73 (MQ=255)
                        ccTGTTCTATGGCTCCAGATGACAAAAATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                           <  2:218313/73‑1 (MQ=255)
                         cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacc                                                      >  2:193579/1‑77 (MQ=255)
                         cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacc                                                      <  1:193579/77‑1 (MQ=255)
                             tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGa                                     <  2:178176/90‑1 (MQ=255)
                                  ggCTCCAGATGACAAACATGACCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGAc                                <  2:557212/90‑1 (MQ=255)
                                     tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCAt                                                                     >  2:591223/1‑50 (MQ=35)
                                     tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCAt                                                                     <  1:591223/50‑1 (MQ=35)
                                     tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaaa                             >  2:201773/1‑90 (MQ=255)
                                     tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaaa                             >  2:512485/1‑90 (MQ=255)
                                            gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAAc                      <  2:309390/90‑1 (MQ=255)
                                            gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAAc                      <  1:450942/90‑1 (MQ=255)
                                                  cATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAAcc                <  2:455240/90‑1 (MQ=255)
                                                        ctcATATCAGGGACTTGTTCGCACCTTCCATAACGc                                                                >  1:77055/1‑36 (MQ=255)
                                                        ctcATATCAGGGACTTGTTCGCACCTTCCATAACGc                                                                <  2:77055/36‑1 (MQ=255)
                                                                aGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTg  <  2:500992/90‑1 (MQ=255)
                                                                                    |                                                                     
GTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTG  >  NZ_CP009273/2283508‑2283661

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