New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 1391494NA (NA)8 (0.460) 7/166 0.8 100% intergenic (+145/‑128) BW25113_RS06980/ynaJ IS5‑like element IS5 family transposase/DUF2534 family protein
?NZ_CP009273 1972747 = 0 (0.000)intergenic (‑293/‑487) flhD/uspC flagellar transcriptional regulator FlhD/universal stress protein UspC

AACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGA  >  NZ_CP009273/1972658‑1972834
                                                                                         |                                                                                       
gcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcc                                                                                         <  1:190981‑M2/1‑1 (MQ=255)
  atttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcctt                                                                                       <  2:146610‑M2/3‑1 (MQ=255)
         cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                                                                                  <  1:189647‑M1/88‑4 (MQ=255)
         cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                                                                                  >  2:189647‑M1/1‑85 (MQ=255)
                 cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaac                                                                        <  1:232672‑M1/90‑14 (MQ=255)
                     cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                    >  1:114852‑M1/1‑73 (MQ=255)
                        aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                 <  1:10079‑M1/90‑21 (MQ=255)
                                gaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttttgttt                                                         <  2:454848‑M2/33‑1 (MQ=255)
                                 aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaac                                                                        <  1:251372‑M1/74‑14 (MQ=255)
                                 aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaac                                                                        >  2:251372‑M1/1‑61 (MQ=255)
                                  ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga                                                       >  2:319051‑M1/1‑60 (MQ=255)
                                  agaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttttgttttg                                                       <  1:447473‑M2/35‑1 (MQ=255)
                                    aaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTggg                                                     <  2:153316‑M2/37‑1 (MQ=255)
                                      atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgtg                                                   <  1:72468‑M2/39‑1 (MQ=255)
                                                        ccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaa                                 <  1:320275‑M2/57‑1 (MQ=255)
                                                                        cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                 >  1:291264‑M1/1‑22 (MQ=255)
                                                                        cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                 >  2:171090‑M1/1‑22 (MQ=255)
                                                                                     cttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt    <  1:218998‑M2/86‑1 (MQ=255)
                                                                                       tcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa  >  2:336921‑M2/3‑90 (MQ=255)
                                                                                         |                                                                                       
AACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGA  >  NZ_CP009273/1972658‑1972834

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.