New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727510 (0.000)17 (0.940) 11/170 0.4 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:91520‑M2/6‑1 (MQ=255)
   tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaa                                                                                       <  2:423844‑M1/90‑5 (MQ=255)
   gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                                                       <  2:183924‑M2/9‑1 (MQ=255)
   gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                                                       <  2:438834‑M2/9‑1 (MQ=255)
     gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                                                     <  1:474732‑M2/11‑1 (MQ=255)
       gaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATg                                                                                   >  2:205632‑M2/78‑90 (MQ=255)
         cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                                                                                 >  2:469162‑M1/1‑80 (MQ=255)
                   aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                       >  1:165924‑M1/1‑70 (MQ=255)
                   aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                       >  1:346928‑M1/1‑70 (MQ=255)
                    tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccct                                                                      <  1:466279‑M1/90‑22 (MQ=255)
                         cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatat                                                                 <  2:73743‑M1/90‑27 (MQ=255)
                               ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatc                                                           <  1:210377‑M1/90‑33 (MQ=255)
                                aatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgt                                                          >  2:278740‑M2/53‑90 (MQ=255)
                                 atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgtg                                                         <  2:424367‑M2/39‑1 (MQ=255)
                                 atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgtg                                                         <  1:89879‑M2/39‑1 (MQ=255)
                                 atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgtg                                                         <  1:205632‑M2/39‑1 (MQ=255)
                                  gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatg                                                        >  1:185225‑M1/1‑55 (MQ=255)
                                      cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttg                                                    <  2:243526‑M1/90‑40 (MQ=255)
                                        gactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                  >  2:34949‑M2/45‑90 (MQ=255)
                                         actgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                                 <  2:170814‑M2/47‑1 (MQ=255)
                                                  atGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagcc                                        >  2:480420‑M1/1‑39 (MQ=255)
                                                      agaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                    >  1:453515‑M2/31‑90 (MQ=255)
                                                        aagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATTCTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                  <  1:270129‑M2/62‑1 (MQ=255)
                                                         aTGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaac                                 <  1:195643‑M1/90‑59 (MQ=255)
                                                                       tattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                   <  2:297923‑M2/77‑1 (MQ=255)
                                                                          cacTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtc                >  1:90131‑M1/1‑15 (MQ=255)
                                                                                      tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                         <  1:478346‑M1/53‑51 (MQ=255)
                                                                                      tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                         >  2:478346‑M1/1‑3 (MQ=255)
                                                                                        aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  2:348785‑M1/1‑1 (MQ=255)
                                                                                        aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  <  2:57481‑M1/90‑90 (MQ=255)
                                                                                        aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  1:173502‑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 < 31 ≤ 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.