New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	1504392 =	13 (0.480)	5 (0.220)	4/248	2.6	38.5%	intergenic (+13/‑153)	trpS/zapE	tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
	?	NC_002947	1504475 =	5 (0.220)					intergenic (+96/‑70)	trpS/zapE	tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

GCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGC  >  NC_002947/1504356‑1504623                                                                                                                        |                                                                                                                                                     gCCTGGGCACCACCGCCGGCCTGTGGCTGCCAACGGGTGGGTGCGGGGGTGCCCGCGAAACAGGCGCCGGGGTGCCGGGCCCGGGCTTTGCGGGGGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCt                                                                                                                        <  2:104248/150‑1 (MQ=25)               gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCCCGAAACAGGCGCCCCCGGGCATGGCACCCGCGTTGCCCGTGTTCGCGGGCACCCCCCCGCCCACCCACAGATTGCCATTAAGCGGGCCGATCGCGAAAGTGACGGcccc                                                                                                          >  2:86358/1‑149 (MQ=255)                                            CGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                                                             >  2:553233‑M1/1‑150 (MQ=255)                                             ggggggccccGCCAACCGGGCGCCGGGGGGCTGGGCAGGGGTTTGGCGGGGGGCCGCGGCCCCCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAc                                                                            <  1:54460/143‑1 (MQ=255)                                                   GCCCGCGAACCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt                                                                      >  2:450560‑M1/1‑150 (MQ=255)                                                   GCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt                                                                      >  2:450533‑M1/1‑150 (MQ=255)                                                     ccgcgaacaccggcaaagccggtgccatgcaccgcggcgcctggttcgcgggcacgcccgctcccacAGACTGATTGCCAATAAGCGGGCCTATCGCTACAGTGACGGCCCGTGTTGTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                                                                    >  2:215811‑M2/68‑150 (MQ=255)                                                      cgcgaacaccggcaaagccggtgccatgcaccgcggcgcctgtttcgcgggcacgcccgctcccacAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                                                   >  1:127814‑M2/67‑150 (MQ=255)                                                      cgcgaacaccggcaaagccggtgccatgcaccgcggcgcctgtttcgcgggcacgcccgctcccacAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCGTAGa                                                                   >  1:171102‑M2/67‑150 (MQ=255)                                                               ggCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCa                                                          >  1:161454/1‑150 (MQ=255)                                                                ggcaaagccggtgccatgcaccgcggcgcctgtttcgcgggcacgcccgctcccacAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtt                                                                                                        >  1:151824‑M2/57‑103 (MQ=255)                                                                ggcaaagccggtgccatgcaccgcggcgcctgtttcgcgggcacgcccgctcccacAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtt                                                                                                        <  2:151824‑M2/47‑1 (MQ=255)                                                                      cgGGGCCTGGCCCCGGCTGGGCGGGTGTGCGCGGCCCCCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                                                   <  1:29643/150‑1 (MQ=255)                                                                                                        gcacgcccgctcccacAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGgcgc                 <  1:292904‑M2/134‑1 (MQ=255)                                                                                                                     cacAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGt    >  1:158853‑M2/4‑150 (MQ=255)                                                                                                                     cacAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGt    >  1:158854‑M2/4‑150 (MQ=255)                                                                                                                     cacAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGt    >  2:354133‑M2/4‑150 (MQ=255)                                                                                                                       cAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGc  >  2:185998‑M2/2‑150 (MQ=255)                                                                                                                        |                                                                                                                                                     GCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGC  >  NC_002947/1504356‑1504623

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