New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727512 (0.100)20 (1.000) 12/170 0.4 90.9% 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

CGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACA  >  NZ_CP009273/1972667‑1972838
                                                                                    |                                                                                       
cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                                                                                    >  2:478727‑M1/1‑85 (MQ=255)
          cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                                                                          >  1:5637‑M1/1‑75 (MQ=255)
            cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                        >  1:175418‑M1/1‑73 (MQ=255)
                    aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                >  1:461349‑M1/1‑65 (MQ=255)
                      ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatg                                                              >  1:108208‑M1/1‑63 (MQ=255)
                              gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatg                                                      >  1:440386‑M1/1‑55 (MQ=255)
                                             aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                       <  2:461349‑M1/90‑51 (MQ=255)
                                             aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                       <  1:565420‑M1/90‑51 (MQ=255)
                                              atGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagcc                                      >  1:436563‑M1/1‑39 (MQ=255)
                                                               cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                          <  2:562522‑M1/85‑64 (MQ=255)
                                                               cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                          >  1:562522‑M1/1‑22 (MQ=255)
                                                                  aaaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatg                  <  1:564422‑M1/90‑72 (MQ=255)
                                                                      cacTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtc              >  2:133997‑M1/1‑15 (MQ=255)
                                                                                 ttAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc   <  2:166346/90‑1 (MQ=255)
                                                                                  tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct  <  1:265963‑M1/90‑88 (MQ=255)
                                                                                  tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                         >  2:438808‑M1/1‑3 (MQ=255)
                                                                                  tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                         <  1:438808‑M1/51‑49 (MQ=255)
                                                                                  tAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa  <  1:309777/90‑1 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccataggacagggttcatcatg                  >  2:35730‑M1/1‑1 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccataggacagggttcatcatg                  <  1:35730‑M1/72‑72 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacag                             >  2:365541‑M1/1‑1 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacag                             <  1:365541‑M1/61‑61 (MQ=255)
                                                                                    |                                                                                       
CGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACA  >  NZ_CP009273/1972667‑1972838

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