New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5201627 =	11 (0.840)	4 (0.420)	4/200	0.8	42.0%	intergenic (‑65/‑115)	PP_4582/PP_4583	membrane protein/putative Peptidase
	?	NC_002947	5201722 =	3 (0.320)					intergenic (‑160/‑20)	PP_4582/PP_4583	membrane protein/putative Peptidase

AGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCG  >  NC_002947/5201618‑5201835                                                                                                         |                                                                                                                  aGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTAGAGCGGGTTTACCCGCGAAGACGCCGGCACAGGCAACTCCCATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGAt                                                                      >  1:295618/1‑150 (MQ=255)                GTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCa                                                       <  2:202982‑M1/150‑1 (MQ=255)                        CTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGc                                                                                   >  1:210348‑M1/1‑114 (MQ=255)                        CTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGc                                                                                   <  2:210348‑M1/114‑1 (MQ=255)                                    AAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCt                                   >  1:202982‑M1/1‑150 (MQ=255)                                      acccgctccaacagggatcgcaccgaggtaggagcgggtttacccgcgaagaggccggcacaggcaaCTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGc                                                                   <  1:111723‑M2/49‑1 (MQ=255)                                      acccgctccaacagggatcgcaccgaggtaggagcgggtttacccgcgaagaggccggcacaggcaaCTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGc                                                                   >  2:111723‑M2/68‑116 (MQ=255)                                       CCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAGGCCCCGACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGc                                <  2:272407‑M1/150‑1 (MQ=255)                                               aacagggatcgcaccgaggtaggagcgggtttacccgcgaagaggccggcacaggcaaCTCCAATCAAGAGACCTTCAATGccaa                                                                                         >  1:35347‑M2/59‑85 (MQ=255)                                               aacagggatcgcaccgaggtaggagcgggtttacccgcgaagaggccggcacaggcaaCTCCAATCAAGAGACCTTCAATGccaa                                                                                         <  2:35347‑M2/27‑1 (MQ=255)                                                      gTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGc                                                             <  1:188267/106‑1 (MQ=255)                                                      gTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGc                                                             >  2:188267/1‑106 (MQ=255)                                                                     gagcgggtttacccgcgaagaggccggcacaggcaaCTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCg  >  2:88719‑M2/37‑150 (MQ=255)                                                                                                       aaCTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAgcgc                     >  1:221160‑M2/3‑97 (MQ=255)                                                                                                       aaCTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAgcgc                     <  2:221160‑M2/95‑1 (MQ=255)                                                                                                         |                                                                                                                  AGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCG  >  NC_002947/5201618‑5201835

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