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

TACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972663‑1972840
                                                                                        |                                                                                         
tacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAAc                                                                                          <  1:93433‑M2/6‑1 (MQ=255)
tacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAAc                                                                                          <  1:142274‑M2/6‑1 (MQ=255)
tacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAAc                                                                                          <  1:23663‑M2/6‑1 (MQ=255)
   gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                                                       <  1:130095‑M2/9‑1 (MQ=255)
   gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                                                       <  2:244762‑M2/9‑1 (MQ=255)
    ggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaa                                                                                      <  1:281805‑M2/10‑1 (MQ=255)
    cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                                                                                      >  2:297522‑M1/1‑85 (MQ=255)
     gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                                                     <  1:242777‑M2/11‑1 (MQ=255)
         cTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                                                                                 >  1:237803‑M1/1‑80 (MQ=255)
                cacaTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                          >  2:267101‑M1/1‑73 (MQ=255)
                cacaTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                          <  1:339486‑M1/90‑18 (MQ=255)
                cacaTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                          >  1:251503‑M1/1‑73 (MQ=255)
                   aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                       <  2:252477‑M1/90‑21 (MQ=255)
                        aCGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                  <  2:251503‑M1/90‑26 (MQ=255)
                          gggggggggctgCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatg                                                                >  2:358381‑M1/4‑63 (MQ=255)
                            gggggggcTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgag                                                              >  1:30617‑M1/2‑61 (MQ=255)
                            gggggggcTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaaggccctgatatgag                                                              >  1:630027‑M1/2‑61 (MQ=255)
                             agaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttttgttttg                                                             <  1:425350‑M2/35‑1 (MQ=255)
                                 atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgtg                                                         <  2:555938‑M2/39‑1 (MQ=255)
                                 atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgtg                                                         <  1:120414‑M2/39‑1 (MQ=255)
                                    tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtt                                                      <  1:267101‑M1/90‑38 (MQ=255)
                                          ctgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                >  1:269920‑M2/43‑90 (MQ=255)
                                            tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatct                                              <  2:15441‑M1/90‑46 (MQ=255)
                                            tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatct                                              <  2:99975‑M1/90‑46 (MQ=255)
                                               ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                           <  1:161909‑M1/90‑49 (MQ=255)
                                                   tGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagcca                                       <  1:11410‑M1/90‑53 (MQ=255)
                                                    gCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                                      <  2:237803‑M1/90‑54 (MQ=255)
                                                          tGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaaca                                <  2:587842‑M1/90‑60 (MQ=255)
                                                           gCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacag                               <  2:630027‑M1/90‑61 (MQ=255)
                                                            cAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagg                              <  1:297522‑M1/90‑62 (MQ=255)
                                                            atttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                              >  2:147341‑M2/25‑90 (MQ=255)
                                                              gATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                            <  2:168794‑M1/90‑64 (MQ=255)
                                                                   cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                       >  1:556640‑M1/1‑22 (MQ=255)
                                                                   cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                       <  1:358381‑M1/90‑69 (MQ=255)
                                                                      ttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                    >  1:39810‑M2/15‑90 (MQ=255)
                                                                       aaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                  >  2:443270‑M1/1‑18 (MQ=255)
                                                                       aaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                  <  1:443270‑M1/43‑26 (MQ=255)
                                                                            cTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcat              >  1:405778‑M1/1‑13 (MQ=255)
                                                                            cTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcat              >  1:451653‑M1/1‑13 (MQ=255)
                                                                                aTTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaac          <  2:404433‑M1/90‑82 (MQ=255)
                                                                                 ttACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                                      >  1:339465‑M1/1‑8 (MQ=255)
                                                                                 ttACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                                      <  2:339465‑M1/61‑54 (MQ=255)
                                                                                   aCTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaactta       >  1:340217‑M1/1‑6 (MQ=255)
                                                                                    cTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttac      >  2:387870‑M1/1‑5 (MQ=255)
                                                                                      tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct    >  1:472599‑M1/1‑3 (MQ=255)
                                                                                      tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct    >  1:48689‑M1/1‑3 (MQ=255)
                                                                                        aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  2:634954‑M1/1‑1 (MQ=255)
                                                                                        aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  <  2:86236‑M1/90‑90 (MQ=255)
                                                                                        |                                                                                         
TACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972663‑1972840

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