New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 707978	NA (NA)	3 (0.030) +GGATTGTTTGCTGTTTCT	4/142	NT	NA	intergenic (+221/‑356)	glnS/chiP	glutamyl‑tRNA synthetase/chitoporin, uptake of chitosugars
	?	NC_000913	4137788 =	NA (NA)					intergenic (‑131/+144)	yijF/gldA	DUF1287 family protein/glycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase

TATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/707881‑707978 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATAT  >  NC_000913/4137788‑4137845                                                                                                   ||||||||||||||||||                                                           TATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCTGAACGCGGCGTTTTGG                                                                            <  2:356908/100‑1 TATTGAGTTTCCGTGCTTTTGTAGGCCGGATAAGTCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCTGAACGCGGCGTTTTG                                                                             <  1:1459505/99‑1    TGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCTGAACGCGGCGTTTTGGAT                                                                          <  2:2953419/99‑1                    TGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCTGAACGCGGCGTTTTGGATTGTTTGCTGTTTCTTGC                                                         <  1:1943803/100‑1                                                                        CCAAAACGCCGCGTTCAGCGGCGTTTTGGATTGTTTGCTGTTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAACTCAATA     >  1:2207664/1‑100                                                                         CAAAACGCCGCGTTCAGCGGCGTTTTGGATTGTTTGCTGTTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAACTCAATAT    >  2:882669/1‑100                                                                           AAACGCCGCGTTCAGCGGCGTTTTGGATTGTTTGCTGTTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAACTCAATATAT  >  1:1338949/1‑100                                                                                                   ||||||||||||||||||                                                           TATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAACAGCAAACAATCCAAAACGCCGCGTTCAGCGGCGTTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/707881‑707978 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATAT  >  NC_000913/4137788‑4137845

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