New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3950464 =	71 (0.980)	17 (0.240)	9/68	2.8	25.6%	intergenic (+44/‑43)	ilvL/ilvX	ilvXGMEDA operon leader peptide/uncharacterized protein IlvX
	?	NC_000913	= 4631789	28 (0.390)					coding (1057/1938 nt)	slt	soluble lytic murein transglycosylase

AGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTG  >  NC_000913/4631755‑4631822                                   |                                  aGCGAAAGAGAAAGATGAATGGCGTTACTGGCAgg                                   >  1:7800832/1‑35 (MQ=255) aGCGAAAGAGAAAGATGAATGGCGTTACTGGCAgg                                   >  1:1626305/1‑35 (MQ=255) aGCGAAAGAGAAAGATGAATGGCGTTACTGGCAgg                                   >  1:3660146/1‑35 (MQ=255)   cGAAAGAGAAAGATGAATGGCGTTACTGGCAGGgg                                 <  1:7557459‑M2/35‑3 (MQ=255)    gaaagagaaagaTGAATGGCGTTACTGGCAGGggg                                <  1:449134‑M2/35‑4 (MQ=255)     aaagagaaagaTGAATGGCGTTACTGGCAGGgggg                               <  1:444287‑M2/35‑5 (MQ=255)     aaagagaaagaTGAATGGCGTTACTGGCAGGgggg                               <  1:4191996‑M2/35‑5 (MQ=255)     aaagagaaagaTGAATGGCGTTACTGGCAGGgggg                               <  1:6487079‑M2/35‑5 (MQ=255)      aagagaaagaTGAATGGCGTTACTGGCAGGggggg                              <  1:7116817‑M2/35‑6 (MQ=255)      aagagaaagaTGAATGGCGTTACTGGCAGGggggg                              <  1:5644267‑M2/35‑6 (MQ=255)       agagaaagaTGGATGGCGTTACTGGCAGGgggggg                             <  1:5394456‑M2/35‑7 (MQ=255)       agagaaagaTGAATGGCGTTACTGGCAGGgggggg                             <  1:4381184‑M2/35‑7 (MQ=255)       agagaaagaTGAATGGCGTTACTGGCAGGCGGAtt                             <  1:3488254/35‑1 (MQ=255)          gaaagaTGAATGGCGTTACTGGCAGGggggggttt                          <  1:2684545‑M2/35‑10 (MQ=255)          gaaagaTGAATGGCGTTACTGGCAGGCGGATTTAt                          <  1:2906331/35‑1 (MQ=255)           aaagaTGAATGGCGTTACTGGCAGGgggggttttt                         <  1:7226863‑M2/35‑11 (MQ=255)            aagaTGAATGGCGTTACTGGCAGGgggggtttttt                        <  1:6907900‑M2/35‑12 (MQ=255)             agaTGAATGGCGTTACTGGCAGGggggtgtttttg                       <  1:7748630‑M2/35‑13 (MQ=255)              gaTGAATGGCGTTACTGGCAGGggggtttattttt                      <  1:184535‑M2/35‑14 (MQ=255)              gaTGAATGGCGTTACTGGCAGGgggggtttttttt                      <  1:5529506‑M2/35‑14 (MQ=255)               aTGAATGGCGTTACTGGCAGGgggggttttttttg                     <  1:5557622‑M2/35‑15 (MQ=255)               aTGAATGGCGTTACTGGCAGGcgggggtttttttg                     <  1:5457279‑M2/35‑15 (MQ=255)                tGAATGGCGTTACTGGCAGGCGGATTTATTGCTgg                    <  1:821111/35‑1 (MQ=255)                tGAATGGCGTTACTGGCAGGCGGATTTATTGCTgg                    <  1:5199047/35‑1 (MQ=255)                tGAATGGCGTTACTGGCAGGCGGATTTATTGCTgg                    <  1:1258159/35‑1 (MQ=255)                  aaTGGCGTTACTGGCAGGCGGATTTATTGCTGGaa                  <  1:2647507/35‑1 (MQ=255)                      gCGTTACTGGCAGGCGGATTTATTGCTGGAAcgcg              <  1:1887932/35‑1 (MQ=255)                      gCGTTACTGGCAGGCGGATTTATTGCTGGAAcgcg              <  1:4169617/35‑1 (MQ=255)                      gCGTTACTGGCAGGCGGATTTATTGCTGGAAcgcg              <  1:3085778/35‑1 (MQ=255)                        gTTACTGGCAGGCGGATTTATTGCTGGAACGCGGa            <  1:782028/35‑1 (MQ=255)                        gTTACTGGCAGGCGGATTTATTGCTGGAACGCGGa            <  1:7383191/35‑1 (MQ=255)                        gTTACTGGCAGGCGGATTTATTGCTGGAACGCGGa            <  1:3241996/35‑1 (MQ=255)                          tACTGGCAGGCGGATTTATTGCTGGAACGCGGACg          <  1:6092208/35‑1 (MQ=255)                            cTGGCAGGCGGATTTATTGCTGGAACGCGGACGtg        <  1:4729922/35‑1 (MQ=255)                            cTGGCAGGCGGATTTATTGCTGGAACGCGGACGtg        <  1:2137262/35‑1 (MQ=255)                            cTGGCAGGCGGATTTATTGCTGGAACGCGGACGtg        <  1:5793860/35‑1 (MQ=255)                            cTGGCAGGCGGATTTATTGCTGGAACGCGGACGtg        <  1:4970730/35‑1 (MQ=255)                            cTGGCAGGCGGATTTATTGCTGGAACGCGGACGtg        <  1:4704103/35‑1 (MQ=255)                             tGGCAGGCGGATTTATTGCTGGAACGCGGACGtga       >  1:1929048/1‑35 (MQ=255)                              ggCAGGCGGATTTATTGCTGGAACGCGGACGtgaa      >  1:2266038/1‑35 (MQ=255)                              ggCAGGCGGATTTATTGCTGGAACGCGGACGtgaa      <  1:3494060/35‑1 (MQ=255)                              ggCAGGCGGATTTATTGCTGGAACGCGGACGtgaa      <  1:4479472/35‑1 (MQ=255)                               gCAGGCGGATTTATTGCTGGAACGCGGACGtgaag     >  1:3340757/1‑35 (MQ=255)                               gCAGGCGGATTTATTGCTGGAACGCGGACGtgaag     >  1:6780615/1‑35 (MQ=255)                               gCAGGCGGATTTATTGCTGGAACGCGGACGtgaag     >  1:1553052/1‑35 (MQ=255)                                  ggCGGATTTATTGCTGGAACGCGGACGtgaagctg  >  1:2346668/1‑35 (MQ=255)                                  ggCGGATTTATTGCTGGAACGCGGACGtgaagctg  <  1:7424215/35‑1 (MQ=255)                                  ggCGGATTTATTGCTGGAACGCGGACGtgaagctg  >  1:1578088/1‑35 (MQ=255)                                   |                                  AGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTG  >  NC_000913/4631755‑4631822

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