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

GTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATG  >  NZ_CP009273/1972668‑1972833
                                                                                   |                                                                                  
gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaattcccccgcttattcgcaccttcCTTAACATaaa                                                                              <  2:538766‑M2/11‑1 (MQ=255)
      acacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaa                                                                        >  1:389294‑M1/1‑78 (MQ=255)
           cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                   >  2:75712‑M1/1‑73 (MQ=255)
              tcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATAc                                                                <  1:5007‑M2/25‑1 (MQ=255)
              aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                >  2:586970‑M1/1‑70 (MQ=255)
              aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                <  2:266349‑M1/90‑21 (MQ=255)
              aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaa                                                                                >  1:128687‑M1/1‑70 (MQ=255)
              aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaa                                                                                <  2:128687‑M1/74‑5 (MQ=255)
                 cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctga                                                             >  1:441447‑M1/1‑67 (MQ=255)
                       aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgag                                                       >  2:105123‑M1/1‑61 (MQ=255)
                       aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgag                                                       >  2:320799‑M1/1‑61 (MQ=255)
                       aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                >  1:104662‑M1/1‑61 (MQ=255)
                       aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                <  2:104662‑M1/81‑21 (MQ=255)
                            atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgtg                                                  <  2:521035‑M2/39‑1 (MQ=255)
                               tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtt                                               >  1:608104‑M1/1‑53 (MQ=255)
                                   gactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgtg                                                  >  2:624478‑M2/45‑83 (MQ=255)
                                   gactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgtg                                                  <  1:624478‑M2/39‑1 (MQ=255)
                                      tgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                        <  1:133433‑M2/49‑1 (MQ=255)
                                          ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                    <  2:608104‑M1/90‑49 (MQ=255)
                                           cagccgagaagaatttccccgcttattcgcaccttcCTTAACATaa                                                                               >  2:256570‑M2/37‑46 (MQ=255)
                                           cagccgagaagaatttccccgcttattcgcaccttcCTTAACATaa                                                                               <  1:256570‑M2/10‑1 (MQ=255)
                                            agccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaa                                  >  2:225431‑M2/36‑90 (MQ=255)
                                             gccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                 >  1:102451‑M2/35‑90 (MQ=255)
                                             atGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagcc                                 >  2:160950‑M1/1‑39 (MQ=255)
                                               gCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                               <  1:128022‑M1/90‑54 (MQ=255)
                                                       cAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagg                       >  1:225433‑M1/1‑29 (MQ=255)
                                                        tttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                      <  2:503749‑M2/67‑1 (MQ=255)
                                                        tttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                      <  1:523490‑M2/67‑1 (MQ=255)
                                                              cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                >  2:355996‑M1/1‑22 (MQ=255)
                                                                 aaaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatg             <  2:225433‑M1/90‑72 (MQ=255)
                                                                  aaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcgtcatga            >  1:398289‑M1/1‑18 (MQ=255)
                                                                     cacTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgt                                            <  2:564038‑M1/55‑41 (MQ=255)
                                                                     cacTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgt                                            >  1:564038‑M1/1‑15 (MQ=255)
                                                                            ttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg  >  1:364043‑M2/4‑90 (MQ=255)
                                                                                   |                                                                                  
GTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATG  >  NZ_CP009273/1972668‑1972833

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.