New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19030724 (0.210)18 (0.960) 15/168 0.2 81.8% intergenic (‑49/‑110) mgrB/yobH PhoP/PhoQ regulator MgrB/YobH family protein
?NZ_CP009273 = 2096424 NA (NA)intergenic (‑68/+3) BW25113_RS10605/wbbL IS5‑like element IS5 family transposase/N‑acetylglucosaminyl‑diphospho‑decaprenol L‑rhamnosyltransferase WbbL

AGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATC  >  NZ_CP009273/2096343‑2096510
                                                                                 |                                                                                      
aGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATccc                                                                                <  2:259936/90‑1 (MQ=33)
                  aaCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAAt                                                              >  1:290917/1‑90 (MQ=255)
                        gTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATc                                                        >  2:91198/1‑90 (MQ=255)
                            tATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGaa                                                    <  2:149214/90‑1 (MQ=255)
                               ggCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGc                                                 <  2:290917/90‑1 (MQ=255)
                                              aCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCtt                                  <  2:96141/90‑1 (MQ=255)
                                                           tCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCAt                     >  1:500738/1‑90 (MQ=255)
                                                              gggACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATg                    >  1:340803/1‑88 (MQ=255)
                                                              gggACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATg                    <  2:340803/88‑1 (MQ=255)
                                                              gggACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTa                  >  2:311776/1‑90 (MQ=255)
                                                               ggACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTaa                 <  1:84650/90‑1 (MQ=255)
                                                               ggACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTaa                 <  2:500738/90‑1 (MQ=255)
                                                                   ttGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCAAGTAAAGAt             >  2:89645/1‑90 (MQ=255)
                                                                          cACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCtt      <  1:266099/90‑1 (MQ=255)
                                                                              ttccttAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCAtc  >  2:114022/1‑90 (MQ=255)
                                                                              ttccttAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCAtc  >  2:154748/1‑90 (MQ=255)
                                                                              ttccttAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCAtc  >  1:304705/1‑90 (MQ=255)
                                                                                 |                                                                                      
AGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATC  >  NZ_CP009273/2096343‑2096510

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