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

GTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972662‑1972840
                                                                                         |                                                                                         
gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTaa                                                                                           >  2:853711/1‑90 (MQ=255)
     cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                                                                                      >  2:693877‑M1/1‑85 (MQ=255)
     cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                                                                                      >  2:594508‑M1/1‑85 (MQ=255)
     cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                                                                                      >  2:533727‑M1/1‑85 (MQ=255)
        cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtg                                                                                   >  2:391942‑M1/1‑82 (MQ=255)
        cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtg                                                                                   >  2:186952‑M1/1‑82 (MQ=255)
          cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                                                                                 <  2:571148‑M1/90‑11 (MQ=255)
          aaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGta                                                                                 <  2:77941‑M2/15‑1 (MQ=255)
               cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                                                                            >  2:22599‑M1/1‑75 (MQ=255)
                 cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                          <  2:96979‑M1/90‑18 (MQ=255)
                 cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                          >  1:591463‑M1/1‑73 (MQ=255)
                   caTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtcc                                                                        <  1:170305‑M1/90‑20 (MQ=255)
                    aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                       <  1:693877‑M1/90‑21 (MQ=255)
                    aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                       >  1:446955‑M1/1‑70 (MQ=255)
                     tACAACGGAGTGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccct                                                                      <  2:577522‑M1/90‑22 (MQ=255)
                         aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                  >  1:901673‑M1/1‑65 (MQ=255)
                             aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgag                                                              >  2:12484‑M1/1‑61 (MQ=255)
                              ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga                                                             >  2:40803‑M1/1‑60 (MQ=255)
                              agaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttttgttttg                                                             <  1:659686‑M2/35‑1 (MQ=255)
                               gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagat                                                            <  1:22599‑M1/90‑32 (MQ=255)
                                  atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgtg                                                         <  2:207422‑M2/39‑1 (MQ=255)
                                   gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatg                                                        >  2:692307‑M1/1‑55 (MQ=255)
                                      gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgttt                                                     >  2:276302‑M1/1‑52 (MQ=255)
                                          actgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttttgttttg                                                             <  1:69293‑M2/35‑1 (MQ=255)
                                          actgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttttgttttg                                                             >  2:69293‑M2/44‑78 (MQ=255)
                                             tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatct                                              <  2:591463‑M1/90‑46 (MQ=255)
                                             gagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                              <  2:473480‑M2/50‑1 (MQ=255)
                                                ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                           <  2:683444‑M1/90‑49 (MQ=255)
                                                ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                           <  2:446955‑M1/90‑49 (MQ=255)
                                                             cAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagg                              <  1:168063‑M1/90‑62 (MQ=255)
                                                                    cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                       >  2:245696‑M1/1‑22 (MQ=255)
                                                                    cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                       >  1:22730‑M1/1‑22 (MQ=255)
                                                                       aaaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatg                    <  1:391942‑M1/90‑72 (MQ=255)
                                                                       aaaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatg                    <  1:829015‑M1/90‑72 (MQ=255)
                                                                              tCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaac                                 <  2:641315‑M1/70‑59 (MQ=255)
                                                                              tCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaac                                 >  1:641315‑M1/1‑12 (MQ=255)
                                                                                  ttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg         >  2:480761‑M2/4‑90 (MQ=255)
                                                                                       tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct    >  1:410514‑M1/1‑3 (MQ=255)
                                                                                       tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct    >  2:158278‑M1/1‑3 (MQ=255)
                                                                                        aaggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccct   >  2:622612‑M1/1‑2 (MQ=255)
                                                                                        aaggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggtt                           >  2:326949‑M1/1‑2 (MQ=255)
                                                                                        aaggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggtt                           <  1:326949‑M1/66‑65 (MQ=255)
                                                                                         aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  2:798618‑M1/1‑1 (MQ=255)
                                                                                         aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  2:915182‑M1/1‑1 (MQ=255)
                                                                                         |                                                                                         
GTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972662‑1972840

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