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

TACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972663‑1972840
                                                                                        |                                                                                         
tacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAAc                                                                                          <  2:9398‑M2/6‑1 (MQ=255)
    cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                                                                                        >  1:173636‑M1/1‑85 (MQ=255)
    cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                                                                                        <  2:173636‑M1/88‑4 (MQ=255)
     gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                                                     <  2:238160‑M2/11‑1 (MQ=255)
         cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                                                                                 >  1:375003‑M1/1‑80 (MQ=255)
           aaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGtata                                                                               <  2:357751‑M2/17‑1 (MQ=255)
                cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                          >  1:509306‑M1/1‑73 (MQ=255)
                            aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgag                                                              >  1:84278‑M1/1‑61 (MQ=255)
                             agaaacgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttttgttttg                                                             <  1:92008‑M2/35‑1 (MQ=255)
                               aaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTggg                                                           <  2:493511‑M2/37‑1 (MQ=255)
                                gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaac                                                                              >  1:346044‑M1/1‑57 (MQ=255)
                                gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaac                                                                              <  2:346044‑M1/70‑14 (MQ=255)
                                               ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                           <  2:391640‑M1/90‑49 (MQ=255)
                                                 aatGCGTGATGCAGATCACACATAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                         <  2:489595‑M1/90‑51 (MQ=255)
                                                       tGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccataga                                   >  2:281470‑M1/1‑34 (MQ=255)
                                                           gCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacag                               <  2:475842‑M1/90‑61 (MQ=255)
                                                            atttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCAt                                                                         >  1:325602‑M2/25‑47 (MQ=255)
                                                            atttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCAt                                                                         <  2:325602‑M2/23‑1 (MQ=255)
                                                            atttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                              >  2:483964‑M2/25‑90 (MQ=255)
                                                                tCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttc                          >  2:2053‑M1/1‑25 (MQ=255)
                                                                   cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccct                                                                      <  2:26590‑M1/43‑22 (MQ=255)
                                                                   cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccct                                                                      >  1:26590‑M1/1‑22 (MQ=255)
                                                                       aaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatga                   >  1:101355‑M1/1‑18 (MQ=255)
                                                                                      tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct    <  1:294661‑M1/90‑88 (MQ=255)
                                                                                      tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct    <  1:35284‑M1/90‑88 (MQ=255)
                                                                                        aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  2:80879‑M1/1‑1 (MQ=255)
                                                                                        |                                                                                         
TACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972663‑1972840

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