New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727510 (0.000)31 (1.500) 20/166 0.1 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

CGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972667‑1972840
                                                                                    |                                                                                         
cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                                                                                      >  2:78177‑M1/1‑85 (MQ=255)
  tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                                                                                    >  2:541624‑M1/1‑83 (MQ=255)
      tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAcgaaggtgcga                                                                                <  1:541624‑M1/90‑12 (MQ=255)
        cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaac                                                                              >  1:595630‑M1/1‑77 (MQ=255)
               aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                       <  2:165778‑M1/90‑21 (MQ=255)
                    aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                  >  2:98659‑M1/1‑65 (MQ=255)
                    aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                  <  2:8511‑M1/90‑26 (MQ=255)
                    aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                  >  1:489837‑M1/1‑65 (MQ=255)
                      ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatg                                                                >  2:265794‑M1/1‑63 (MQ=255)
                        aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgag                                                              >  2:4528‑M1/1‑61 (MQ=255)
                         ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga                                                             >  2:145137‑M1/1‑60 (MQ=255)
                         ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga                                                             <  1:351229‑M1/90‑31 (MQ=255)
                                 gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgttt                                                     >  1:474755‑M1/1‑52 (MQ=255)
                                           ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                           <  1:145137‑M1/90‑49 (MQ=255)
                                              atGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagcc                                        >  1:391957‑M1/1‑39 (MQ=255)
                                                        cAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagg                              <  1:494575‑M1/90‑62 (MQ=255)
                                                            tCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttc                          <  2:489837‑M1/90‑66 (MQ=255)
                                                               cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                  <  1:564471‑M1/47‑26 (MQ=255)
                                                               cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                  >  2:564471‑M1/1‑22 (MQ=255)
                                                                 caAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcat                     >  2:341103‑M1/1‑20 (MQ=255)
                                                                   aaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatga                   >  2:409269‑M1/1‑18 (MQ=255)
                                                                      cacTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtc                >  2:422345‑M1/1‑15 (MQ=255)
                                                                                 ttAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                                      <  2:326834‑M1/57‑54 (MQ=255)
                                                                                 ttAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                                      >  1:326834‑M1/1‑4 (MQ=255)
                                                                                  tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct    <  1:265577‑M1/90‑88 (MQ=255)
                                                                                  tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct    <  1:130655‑M1/90‑88 (MQ=255)
                                                                                   aaggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatg                    >  1:504491‑M1/1‑2 (MQ=255)
                                                                                   aaggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatg                    <  2:504491‑M1/73‑72 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  2:483420‑M1/1‑1 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  1:111457‑M1/1‑1 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  2:125394‑M1/1‑1 (MQ=255)
                                                                                    |                                                                                         
CGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972667‑1972840

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.