New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727510 (0.000)31 (0.840) 22/168 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

CTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACA  >  NZ_CP009273/1972665‑1972838
                                                                                      |                                                                                       
cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                                                                                      <  1:611129‑M1/90‑4 (MQ=255)
 gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                                                     <  1:193502‑M2/9‑1 (MQ=255)
 gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                                                     <  1:350669‑M2/9‑1 (MQ=255)
  ggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaa                                                                                    <  2:748830‑M2/10‑1 (MQ=255)
  cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGGTCACACAAAACACTCAATTACTTAAggaag                                                                                    >  1:661475‑M1/1‑85 (MQ=255)
  cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                                                                                    >  2:987234‑M1/1‑85 (MQ=255)
   gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                                                   <  2:662711‑M2/11‑1 (MQ=255)
    tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                                                                                  >  1:886355‑M1/1‑83 (MQ=255)
    tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                                                                                  >  2:648435‑M1/1‑83 (MQ=255)
        tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcga                                                                              <  1:987234‑M1/90‑12 (MQ=255)
            cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                                                                          <  1:512789‑M1/90‑16 (MQ=255)
            cGCACATACAACGGAGGGGGGCTGCAATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                                                                          >  1:171389‑M1/1‑75 (MQ=255)
              cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                        >  2:248507‑M1/1‑73 (MQ=255)
               gctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCAt                                                                       <  2:790050‑M2/23‑1 (MQ=255)
                 tcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATAc                                                                     <  2:148529‑M2/25‑1 (MQ=255)
                   aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctg                                                                   <  1:430655‑M1/90‑23 (MQ=255)
                    aacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttt                                                                  <  2:532444‑M2/28‑1 (MQ=255)
                      aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                >  1:525125‑M1/1‑65 (MQ=255)
                      aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                <  2:872636‑M1/90‑26 (MQ=255)
                              gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatca                                                        >  2:66146‑M1/1‑57 (MQ=255)
                              gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatca                                                        >  1:683692‑M1/1‑57 (MQ=255)
                               atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgtg                                                       <  1:96858‑M2/39‑1 (MQ=255)
                                    cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttg                                                  <  1:66146‑M1/90‑40 (MQ=255)
                                     gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgt                                                 <  2:373444‑M1/90‑41 (MQ=255)
                                      gactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                                                   <  1:1060111‑M2/11‑1 (MQ=255)
                                      gactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                                                   >  2:1060111‑M2/45‑55 (MQ=255)
                                          tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatct                                            <  2:900625‑M1/90‑46 (MQ=255)
                                          tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatct                                            <  2:445824‑M1/90‑46 (MQ=255)
                                               aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgggatcatgtttgtcatctggagc                                       <  2:820167‑M1/90‑51 (MQ=255)
                                                gccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                      >  1:160450‑M2/35‑90 (MQ=255)
                                                  gCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                                    <  1:165345‑M1/90‑54 (MQ=255)
                                                  gCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                                    <  1:3226‑M1/90‑54 (MQ=255)
                                                         gCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacag                             <  2:630165‑M1/90‑61 (MQ=255)
                                                          cAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagg                            <  2:79397‑M1/90‑62 (MQ=255)
                                                          atttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                            <  1:191200‑M2/66‑1 (MQ=255)
                                                          atttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                            >  2:283097‑M2/25‑90 (MQ=255)
                                                                 cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                     >  1:234388‑M1/1‑22 (MQ=255)
                                                                 cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                     >  2:1034042‑M1/1‑22 (MQ=255)
                                                                  gcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                    >  1:770943‑M2/17‑90 (MQ=255)
                                                                         acTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtca             <  2:443130‑M1/90‑77 (MQ=255)
                                                                           tCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                          >  1:917087‑M1/1‑12 (MQ=255)
                                                                           tCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                          <  2:917087‑M1/75‑64 (MQ=255)
                                                                           caccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa           <  1:1011222‑M2/83‑1 (MQ=255)
                                                                                    tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct  >  1:37817‑M1/1‑3 (MQ=255)
                                                                                    tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct  >  2:352321‑M1/1‑3 (MQ=255)
                                                                                    tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct  <  2:1000579‑M1/90‑88 (MQ=255)
                                                                                      |                                                                                       
CTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACA  >  NZ_CP009273/1972665‑1972838

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