New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 2256677	0 (0.000)	7 (0.120)	5/250	NT	100%	coding (34/1302 nt)	PP_1990	transposase
	?	NC_002947	= 2413151	NA (NA)					intergenic (+415/‑6)	rlmM/PP_2114	23S rRNA (cytidine(2498)‑2'‑O)‑methyltransferase RlmM/transposase

GAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTGGAGCAGGCCATTGCGCCGGAGTGGGTCGATCAGGTTTTCGAGGAGCACCGCCAACGGCAGTATTCTCGCGAGCTGCTGTTCTCGACCATTATCAAGTTGATGTCGCTTGTTTCATTGGGCTTGAAGCCATCGCTGCATGCTGCGGCAAGACAACTGGACGAC  >  NC_002947/2256620‑2256823                                                          |                                                                                                                                                   gagaAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTGGAGCAGGCCATTGCGCcggagtgggtgcttttttctcgaatcgcgctagtacatcagcccacatcgctcgggtcttcctgagtggatttcagcgaattctaccg                                                            >  5:153829‑M1/1‑58 (MQ=255)              ccGCCAGCGTAATGGCCAAACTTGCCTTGGAGCAGGCCATTGCGCtggcgggtgcttttttctcgaatcgcgctagtacatcagcccacatcgctcgggtcttcctgactggatttcagcgaattctaccgaagaccttgaaagggatg                                            >  6:280026‑M1/1‑45 (MQ=255)              ccGCCAGCGTAATGGCCAAACTTGCCTTGAAGCAGGCCATTACGCtggcgggtgcttttttctcgaatcgcgctagtacatcagcccacatcgctcgggtcttcctgactggatttcagcgaattctaccgaagaccttgaaagggatg                                            >  1:220288‑M1/1‑45 (MQ=255)                gCCAGCGTAATGGCCAAACTTGCCTTGGAGCAGGCCATTGCGCtggcgggtgcttttttctcgaatcgcgctagtacatcagcccacatcgctcgggtcttcctgactggatttcagcgaattctaccgaagaccttgaaagggatggc                                          >  5:359111‑M1/1‑43 (MQ=255)                  cAGCGTAATGGCCAAACTTGCCTTGGAGCAGGCCATTGCGCtggcgggtgcttttttctcgaatcgcgctagtacatcagcccacatcgctcgggtcttcctgactggatttcagcgaattctaccgaagaccttgaaagggatggccc                                        >  2:293145‑M1/1‑41 (MQ=255)                            gCCAAACTTGCCTTGGAGCAGGCCATTGCGCtggcgggtgcttttttctcgaatcgcgctagtacatcagcccacatcgctcgggtcttcctgactggatttcagcgaattctaccgaagaccttgaaagggatggccgtaagggcga                               >  1:91631‑M1/1‑31 (MQ=255)                                      ccTTGGAGCAGGCCATTGCGCtggcgggtgcttttttctcgaatcgcgctagtacatcagcccacatcgctcgggtcttcctgactggatttcagcgaattctaccgaagaccttgaaagggatggcccttacggcgggtcacttttt                     <  6:359111‑M1/148‑128 (MQ=255)                                                         gctggcgggtgcttttttctcgaatcgcgctagtacatcagcccacatcgctcgggtcttcctgactggatttcagcgaattctaccgaagaccttgaaagggatggcccttacggcgggtcactttttgtcaaacgcgacaaaaagt  >  2:125590‑M1/1‑2 (MQ=255)                                                         gctggcgggtgcttttttctcgaatcgcgctagtacatcagcccacatcgctcgggtcttcctgactggatttcagcgaattctaccgaagaacttgaaagggatggcccttacggcgggtcactttttgtcaaacgcgacaaaaatt  >  8:128607‑M1/1‑2 (MQ=255)                                                          |                                                                                                                                                   GAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTGGAGCAGGCCATTGCGCCGGAGTGGGTCGATCAGGTTTTCGAGGAGCACCGCCAACGGCAGTATTCTCGCGAGCTGCTGTTCTCGACCATTATCAAGTTGATGTCGCTTGTTTCATTGGGCTTGAAGCCATCGCTGCATGCTGCGGCAAGACAACTGGACGAC  >  NC_002947/2256620‑2256823

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