New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727511 (0.080)16 (1.320) 12/170 0.2 94.1% 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

GGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAG  >  NZ_CP009273/1972692‑1972839
                                                           |                                                                                        
ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga                                                            <  1:324268‑M1/90‑31 (MQ=255)
  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatc                                                          <  2:313731‑M1/90‑33 (MQ=255)
     gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatg                                                       >  1:197192‑M1/1‑55 (MQ=255)
                     atGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagcc                                       >  2:246725‑M1/1‑39 (MQ=255)
                            aTGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaac                                <  1:104549‑M1/90‑59 (MQ=255)
                               cAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagg                             >  2:122971‑M1/1‑29 (MQ=255)
                                 gATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                           >  2:287229‑M1/1‑27 (MQ=255)
                                   tCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttc                         <  2:328117‑M1/90‑66 (MQ=255)
                                      cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                      >  1:127885‑M1/1‑22 (MQ=255)
                                      cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                      >  2:162578‑M1/1‑22 (MQ=255)
                                                tCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatc            <  2:207703‑M1/90‑79 (MQ=255)
                                                tCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatc            <  2:347090‑M1/90‑79 (MQ=255)
                                                    ttACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaact        >  1:91532‑M1/1‑8 (MQ=255)
                                                      aCTTAAggaaggtgcgtacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaactta      >  2:262190‑M1/1‑6 (MQ=255)
                                                        ttAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                          >  1:59371‑M1/1‑4 (MQ=255)
                                                        ttAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                          <  2:59371‑M1/52‑49 (MQ=255)
                                                          aaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg  >  1:228878/1‑90 (MQ=255)
                                                           |                                                                                        
GGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAG  >  NZ_CP009273/1972692‑1972839

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