New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727510 (0.000)10 (0.600) 8/160 0.6 100% intergenic (‑297/‑483) flhD/uspC flagellar transcriptional regulator FlhD/universal stress protein UspC
?NZ_CP009273 = 2096424 NA (NA)intergenic (‑68/+3) BW25113_RS10605/wbbL IS5‑like element IS5 family transposase/N‑acetylglucosaminyl‑diphospho‑decaprenol L‑rhamnosyltransferase WbbL

CGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972667‑1972840
                                                                                    |                                                                                         
cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                                                                                        >  2:189647‑M1/1‑85 (MQ=255)
cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                                                                                        <  1:189647‑M1/88‑4 (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                       >  2:171090‑M1/1‑22 (MQ=255)
                                                               cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                       >  1:291264‑M1/1‑22 (MQ=255)
                                                                            cttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt          <  1:218998‑M2/86‑1 (MQ=255)
                                                                              tcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa        >  2:336921‑M2/3‑90 (MQ=255)
                                                                                 ttAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacc     >  1:371744‑M1/1‑4 (MQ=255)
                                                                                  tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct    >  1:157193‑M1/1‑3 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  1:359661‑M1/1‑1 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  1:447676‑M1/1‑1 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  1:170515‑M1/1‑1 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  2:427187‑M1/1‑1 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaaca                                >  1:305642‑M1/1‑1 (MQ=255)
                                                                                    aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaaca                                <  2:305642‑M1/60‑60 (MQ=255)
                                                                                    |                                                                                         
CGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972667‑1972840

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.