New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1169369 =	25 (0.760)	3 (0.100)	3/248	3.8	12.3%	intergenic (+39/+43)	bhsA/ldtC	biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
	?	NC_000913	1169396 =	20 (0.670)					intergenic (+66/+16)	bhsA/ldtC	biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AAAGATGCAGCACAAACTGACGCTGAAGTGATGCAACATGTGATGGATGTCCGTTCCGGGATGCCGGTGGATGTCCGCCGTCATCAAGTGAGCCCACAAACGCTGTAAAGCCCGGATAAAAAAATCCCCGCGACACGGCATGTCGCGGGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1169519‑1169369 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgtcgcggggaTTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACG  >  NC_000913/1169396‑1169441                                                                                                                                                                                                       AAAGATGCAGCACAAACTGACGCTGAAGTGATGCAACATGTGATGGATGTCCGTTCCGGGATGCCGGTGGATGTCCGCCGTCATCAAGTGAGCCCACAAACGCTGTAAAGCCCGGATAAAAAAATCCCCGCGACACGGCATGTCGCGGG                                                  >  6:86466/1‑149    GATGCAGCACAAACTGACGCTGAAGTGATGCAACATGTGATGGATGTCCGTTCCGGGATGCCGGTGGATGTCCGCCGTCATCAAGTGAGCCCACAAACGCTGTAAAGCCCGGATAAAAAAATCCCCGCGCCACGGCATGTCGCGGGGA                                                >  5:36244/1‑148                                                     GTTCCGGGATGCCGGTGGATGTCCGCCGTCATCAAGTGAGCCCACAAACGCTGTAAAGCCCGGATAAAAAAATCCCCGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACG  >  3:24061/1‑145                                                     GTTCCGGGATGCCGGTGGATGTCCGCCGTCATCAAGTGAGCCCACAAACGCTGTAAAGCCCGGATAAAAAAATCCCCGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACG  <  4:24061/145‑1                                                                                             CCCACAAACGCTGTAAAGCCCGGATAAAAAAATCCCCGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGAT      <  6:62933/101‑1                                                                                                                                                                                                       AAAGATGCAGCACAAACTGACGCTGAAGTGATGCAACATGTGATGGATGTCCGTTCCGGGATGCCGGTGGATGTCCGCCGTCATCAAGTGAGCCCACAAACGCTGTAAAGCCCGGATAAAAAAATCCCCGCGACACGGCATGTCGCGGGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1169519‑1169369 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgtcgcggggaTTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACG  >  NC_000913/1169396‑1169441

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