New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727511 (0.070)21 (1.530) 16/168 0.1 95.5% 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

TCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAG  >  NZ_CP009273/1972666‑1972839
                                                                                     |                                                                                        
tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaa                                                                                      <  2:99685‑M1/90‑5 (MQ=255)
 cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                                                                                     >  1:128425‑M1/1‑85 (MQ=255)
           cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                                                                           >  2:138149‑M1/1‑75 (MQ=255)
           cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                                                                                       >  2:42192‑M1/1‑75 (MQ=255)
           cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                                                                                       <  1:42192‑M1/78‑4 (MQ=255)
             cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                         <  2:184223‑M1/90‑18 (MQ=255)
             cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATACAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                         >  2:134677‑M1/1‑73 (MQ=255)
                     aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                 >  1:13891‑M1/1‑65 (MQ=255)
                     aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                 <  2:74177‑M1/90‑26 (MQ=255)
                          ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga                                                            >  2:241451‑M1/1‑60 (MQ=255)
                          ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga                                                            >  1:152964‑M1/1‑60 (MQ=255)
                                   cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttg                                                   <  2:287445‑M1/90‑40 (MQ=255)
                                                tGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagcca                                      <  2:334913‑M1/90‑53 (MQ=255)
                                                           gATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                           >  1:341448‑M1/1‑27 (MQ=255)
                                                                  caAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcat                    <  2:152964‑M1/90‑71 (MQ=255)
                                                                   aaaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatcaggagccatagaacagggttcatcatg                   <  1:241451‑M1/90‑72 (MQ=255)
                                                                                aCTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                          >  1:320677‑M1/1‑6 (MQ=255)
                                                                                aCTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                          <  2:320677‑M1/54‑49 (MQ=255)
                                                                                  ttAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc    <  2:347092/90‑1 (MQ=255)
                                                                                   tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                          >  1:56284‑M1/1‑3 (MQ=255)
                                                                                   tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                          <  2:56284‑M1/51‑49 (MQ=255)
                                                                                    aaggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacctt  >  1:234555‑M1/1‑2 (MQ=255)
                                                                                     |                                                                                        
TCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAG  >  NZ_CP009273/1972666‑1972839

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