New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3968570 =	112 (1.320)	5 (0.060)	5/522	NT	4.3%	coding (655/1104 nt)	wecA	UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
	?	NC_000913	= 4499511	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region
Rejected: Frequency below/above cutoff threshold.

GCCAGCACGCCAATGGAAGCGAGCAGTGCTGCGGCAAGGGTAATCAGCACAAACGCCTGACGGGAAGTAAACCCGGCACGCATGATCAAATGGTGAATATGCTGACGGTCAGGAGAGAATGGGCTCATGCCTTTACGCAGGCGACGGTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTATCCACAAAGCGGTAACCGGGCTGATGGGATGGGTTTTGCCCTGGGTCGTTTCGAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3968850‑3968570 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACCGCCGTGTTCTTCGGTTCAAGTCCCAACATCCGGGCGAA  <  NC_000913/4499511‑4499225                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          GCCAGCACGCCAATGGAAGCGAGCAGTGCTGCGGCAAGGGTAATCAGCACAAACGCCTGACGGGAAGTAAACCCGGCACGCATGATCAAATGGTGAATATGCTGACGGTCAGGAGAGAATGGGCTCATGCCTTTACGCAGGCGACGGTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTATCCACAAAGCGGTAACCGGGCTGATGGGATGGGTTTTGCCCTGGGTCGTTTCGAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCCTGGATTTG                                                                                                                                                                                                                                                                                         >  2:350971/1‑289                                     AGGGTAATCAGCACAAACGCCTGACGGGAAGTAAACCCGGCACGCATGATCAAATGGTGAATATGCTGACGGTCAGGAGAGAATGGGCTCATGCCTTTACGCAGGCGACGGTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTATCCACAAAGCGGTAACCGGGCTGATGGGATGGGTTTTGCCCTGGGTCGTTTCGAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACCTAACCC                                                                                                                                                                                                                                                     >  1:594328/1‑289                                                                                         AATGGTGAATATGCTGACGGTCAGGAGAGAATGGGCTCATGCCTTTACGCAGGCGACGGTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTATCCACAAAGCGGTAACCGGGCTGATGGGATGGGTTTTGCCCTGGGTCGTTTCGAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCA                                                                                                                                                                                                 <  1:62545/289‑1                                                                                                                                                                        CATTAGCGGAATGGCGATTATCCACAAAGCGGTAACCGGGCTGATGGGATGGGTTTTGCCCTGGGTCGTTTCGAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCA                                                                                                                  <  1:350971/289‑1                                                                                                                                                                                                                         TGGGTTTTGCCCTGGGTCGTTTCTAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTCCAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTA                                                                  <  2:594328/288‑1                                                                                                                                                                                                                                                                                        CCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACCGCCGTGTTCTTCGGTTCAAGTCCCAACATCCGGGCGAA  <  2:273880/289‑1                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          GCCAGCACGCCAATGGAAGCGAGCAGTGCTGCGGCAAGGGTAATCAGCACAAACGCCTGACGGGAAGTAAACCCGGCACGCATGATCAAATGGTGAATATGCTGACGGTCAGGAGAGAATGGGCTCATGCCTTTACGCAGGCGACGGTACATAATCGCCACCATATCCATTAGCGGAATGGCGATTATCCACAAAGCGGTAACCGGGCTGATGGGATGGGTTTTGCCCTGGGTCGTTTCGAGCAGGATCCAGATAACGGTAAAACCAATCAGCGTACTGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3968850‑3968570 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACCGCCGTGTTCTTCGGTTCAAGTCCCAACATCCGGGCGAA  <  NC_000913/4499511‑4499225

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