New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	223597 =	NA (NA)	13 (0.090)	6/168	NT	NA	intergenic (+189/‑174)	gmhB/rrsH	D,D‑heptose 1,7‑bisphosphate phosphatase/16S ribosomal RNA of rrnH operon
	?	NC_000913	223624 =	NA (NA)					intergenic (+216/‑147)	gmhB/rrsH	D,D‑heptose 1,7‑bisphosphate phosphatase/16S ribosomal RNA of rrnH operon

GTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAACAATTTATCAGA  >  NC_000913/223525‑223648                                                                                                    |                         gttcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATACACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCa                          <  1:3643548/100‑1 (MQ=255)  ttcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCaa                         >  2:207088/1‑100 (MQ=255)   tcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAAc                        <  2:2321959/100‑1 (MQ=255)   tcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAAc                        <  2:2320921/100‑1 (MQ=255)   tcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAAc                        <  1:2530162/100‑1 (MQ=255)     ggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGAGTCGCAACTg                      >  1:3416420/1‑100 (MQ=255)      gttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGc                     <  2:619528/100‑1 (MQ=40)      gttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGc                     >  2:2963638/1‑100 (MQ=40)      gttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGc                     >  1:198054/1‑100 (MQ=40)      gttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGc                     <  1:2082861/100‑1 (MQ=40)       ttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGct                    >  2:2230558/1‑100 (MQ=38)       ttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGct                    >  2:2196465/1‑100 (MQ=38)        tGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGctc                   <  1:3022487/100‑1 (MQ=35)         gAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGctct                  <  1:1860558/100‑1 (MQ=35)           agagAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                >  2:758747/1‑100 (MQ=32)           agagAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                >  1:1068388/1‑100 (MQ=32)           agagAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                >  1:2430504/1‑100 (MQ=32)           agagAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                <  1:2213993/100‑1 (MQ=33)            gagaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                >  2:3837951/1‑99 (MQ=32)            gagaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTa               <  2:1831743/100‑1 (MQ=32)            gagaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTa               <  2:2587412/100‑1 (MQ=32)            gagaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTa               <  1:191253/100‑1 (MQ=32)             agaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTaa              >  1:1334319/1‑100 (MQ=31)              gaAAAATCCTGTAATGCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAAc             <  1:1881203/100‑1 (MQ=2)              gaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAAc             <  1:3973607/100‑1 (MQ=11)              gaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAAc             <  2:2220350/100‑1 (MQ=11)                          aagtCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAACAATTTATCAGa  <  1:3329750/96‑1 (MQ=1)                                                                                                    |                         GTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAACAATTTATCAGA  >  NC_000913/223525‑223648

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