New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727512 (0.070)33 (1.230) 22/172 0.2 94.3% 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:544695‑M2/6‑1 (MQ=255)
 acgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACa                                                                                         <  2:547694‑M2/7‑1 (MQ=255)
   gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                                                       <  2:205560‑M2/9‑1 (MQ=255)
   gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                                                       <  2:747268‑M2/9‑1 (MQ=255)
    cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                                                                                      >  2:407092‑M1/1‑85 (MQ=255)
     gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                                                     <  1:594341‑M2/11‑1 (MQ=255)
     gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                                                     <  1:360271‑M2/11‑1 (MQ=255)
      tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                                                                                    >  1:176145‑M1/1‑83 (MQ=255)
      tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                                                                                    >  2:400841‑M1/1‑83 (MQ=255)
       gaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATg                                                                                   <  1:588926‑M2/13‑1 (MQ=255)
         cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                                                                                 >  2:599211‑M1/1‑80 (MQ=255)
           aaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGtata                                                                               <  2:577212‑M2/17‑1 (MQ=255)
              cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCTATTACTTAAggaaggtgcgaacaa                                                                            >  2:586686‑M1/1‑75 (MQ=255)
                cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                          >  2:602805‑M1/1‑73 (MQ=255)
                 gcttaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCAt                                                                         <  2:275373‑M2/23‑1 (MQ=255)
                  caTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtcc                                                                        >  1:266972‑M1/1‑71 (MQ=255)
                   aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                       >  2:517382‑M1/1‑70 (MQ=255)
                   aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                       <  1:517382‑M1/90‑21 (MQ=255)
                        aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                  <  2:40265‑M1/90‑26 (MQ=255)
                          ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatg                                                                >  2:275668‑M1/1‑63 (MQ=255)
                             agaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttttgttttg                                                             <  1:33286‑M2/35‑1 (MQ=255)
                             agaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttttgttttg                                                             <  2:322743‑M2/35‑1 (MQ=255)
                               aaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTggg                                                           <  2:148078‑M2/37‑1 (MQ=255)
                               aaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTggg                                                           <  1:550652‑M2/37‑1 (MQ=255)
                                gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                                                                            <  2:713182‑M1/72‑16 (MQ=255)
                                gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                                                                            >  1:713182‑M1/1‑57 (MQ=255)
                                       gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgt                                                   <  1:263494‑M1/90‑41 (MQ=255)
                                               ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                           <  1:441564‑M1/90‑49 (MQ=255)
                                               ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                           <  1:702625‑M1/90‑49 (MQ=255)
                                               ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                           <  2:243708‑M1/90‑49 (MQ=255)
                                                 aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                         <  1:602805‑M1/90‑51 (MQ=255)
                                                       gaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                   >  1:189388‑M2/30‑90 (MQ=255)
                                                         aTGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaac                                 <  1:692224‑M1/90‑59 (MQ=255)
                                                               aTCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggtt                           <  2:196699‑M1/90‑65 (MQ=255)
                                                                   cacaAAACACTCAGTTACTTAAggaagatgcgaacaagtccctgatatgagatcatgtttgtcatctagagccatagaacagggttcatc                       >  1:379687‑M1/1‑22 (MQ=255)
                                                                   cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                       >  1:222256‑M1/1‑22 (MQ=255)
                                                                    gcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                      >  1:459571‑M2/17‑90 (MQ=255)
                                                                      aaaaCACTCAATTGCTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatg                    <  1:361510‑M1/90‑72 (MQ=255)
                                                                      aaaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatg                    <  2:748078‑M1/90‑72 (MQ=255)
                                                                      aaaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatg                    <  2:73909‑M1/90‑72 (MQ=255)
                                                                          cacTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                            <  2:540279‑M1/78‑64 (MQ=255)
                                                                          cacTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                            >  1:540279‑M1/1‑15 (MQ=255)
                                                                             tCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatc             <  1:275668‑M1/90‑79 (MQ=255)
                                                                                     ttAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                         >  1:429171‑M1/1‑4 (MQ=255)
                                                                                     ttAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                         <  2:429171‑M1/54‑51 (MQ=255)
                                                                                     ttAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc     <  2:50422/90‑1 (MQ=255)
                                                                                      tAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa    <  1:141058/90‑1 (MQ=255)
                                                                                        aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  <  2:626336‑M1/90‑90 (MQ=255)
                                                                                        aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  2:57363‑M1/1‑1 (MQ=255)
                                                                                        aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  2:138569‑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 < 26 ≤ ATCG/ATCG < 33 ≤ 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.