New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_007779	2519722 =	NA (NA)	3 (0.010)	3/176	NT	1.1%	intergenic (+32/‑55)	nupC/insL	nucleoside (except guanosine) transporter/transposase
	?	NC_007779	= 3175166	282 (1.010)					coding (324/828 nt)	cpdA	cyclic 3',5'‑adenosine monophosphate phosphodiesterase

GTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATG‑gg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_007779/2519791‑2519722 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGGTATCTCCCCGGCGAAGCGCGTGTTTATTGGTGAGCAATGGCAAATCCTGTTGCTGGATAGCCAGGTGTTTGGCGTGCCGCACGGTGAGCTGAGCGAGTT  <  NC_007779/3175166‑3175065                                                                                                                                                                               GTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATG‑GGGTATCTCCCCGGCGAAGCGCGTGTTTATT                                                                         <  1:147950/101‑1                              ATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATG‑GGGTATCTCCCCGGCGAAGCGCGTGTTTATTGGTGAGCAATGGCAAATCCTGTTGCTGGA                                            >  1:8523897/1‑101                                                 AACCCTTAAGTTAGCGCTTATGCGGGTATCTCCCCGGCGAAGCGCGTGTTTATTGGTGAGCAATGGCAAATCCTGTTGCTGGATAGCCAGGTGTTTGGCGT                          >  1:5741343/1‑101                                                                        GGGTATCTCCCCGGCGAAGCGCGTGTTTATTGGTGAGCAATGGCAAATCCTGTTGCTGGATAGCCAGGTGTTTGGCGTGCCGCACGGTGAGCTGAGCGAGT   <  1:1737674/101‑1                                                                        GGGTATCTCCCCGGCGAAGCGCGTGTTTATTGGTGAGCAATGGCAAATCCTGTTGCTGGATAGCCAGGTGTTTGGCGTGCCGCACGGTGAGCTGAGCGAGT   <  1:5224385/101‑1                                                                         GGTATCTCCCCGGCGAAGCGCGTGTTTATTGGTGAGCAATGGCAAATCCTGTTGCTGGATAGCCAGGTGTTTGGCGTGCCGCACGGTGAGCTGAGCGAGTT  <  1:8907955/101‑1                                                                                                                                                                               GTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATG‑gg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_007779/2519791‑2519722 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGGTATCTCCCCGGCGAAGCGCGTGTTTATTGGTGAGCAATGGCAAATCCTGTTGCTGGATAGCCAGGTGTTTGGCGTGCCGCACGGTGAGCTGAGCGAGTT  <  NC_007779/3175166‑3175065

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