New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4055336	167 (1.200)	4 (0.030)	4/428	NT	2.3%	coding (1004/1050 nt)	glnL	sensory histidine kinase in two‑component regulatory system with GlnG
	?	NC_000913	4055340 =	174 (1.260)					coding (1000/1050 nt)	glnL	sensory histidine kinase in two‑component regulatory system with GlnG
Rejected: Frequency below/above cutoff threshold.

GCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4055208‑4055336 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cTGGTAAATTCAATTTTGCCTGAATGCTGATCAATCAAATTACGAGCGATGGATAAGCCAAGCCCGGTGCCACCTTCGCGGCCGCTGACCATCGGGTAAAACAGCGTATCCTGCAAATGAGGCGGAATGCCCGGCCCGTTATCTTCCACATCAATCCGCGCCGCCAGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCG  >  NC_000913/4055340‑4055576                                                                                                                                                                                                                                                                                                                                                                                GCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCTGGTAAATTCAATTTTGCCTGAATGCTGATCAATCAAATTACGAGCGATGGATAAGCCAAGCCCGGTGCCACCTTCGCGGCCGCTGACCATCGGGTAAAACAGCGTAT                                                                                                                                   >  2:242587/1‑237                AGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCTGGTAAATTCAATTTTGCCTGAATGCTGATCAATCAAATTACGAGCGATGGATAAGCCAAGCCCGGTGCCACCTTCGCGGCCGCTGACCATCGGGTAAAACAGCGTATCCTGCAAATGAGGCGGA                                                                                                                  >  1:1114816/1‑239                    CCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCTGGTAAATTCAATTTTGCCTGAATGCTGATCAATCAAATTACGAGCGATGGATAAGCCAAGCCCGGTGCCACCTTCGCGGCCGCTGACCATCGGGTAAAACAGCGTATCCTGCAAATGAGGCGGAATGC                                                                                                              <  1:242587/239‑1                                                                                           TGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCTGGTAAATTCAATTTTGCCTGAATGCTGATCAATCAAATTACGAGCGATGGATAAGCCAAGCCCGGTGCCACCTTCGCGGCCGCTGACCATCGGGTAAAACAGCGTATCCTGCAAATGAGGCGGAATGCCCGGCCCGTTATCTTCCACATCAATCCGCGCCGCCAGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAAC                                       >  2:1008349/1‑239                                                                                                                                 CTGGTAAATTCAATTTTGCCTGAATGCTGATCAATCAAATTACGAGCGATGGATAAGCCAAGCCCGGTGCCACCTTCGCGGCCGCTGACCATCGGGTAAAACAGCGTATCCTGCAAATGAGGCGGAATGCCCGGCCCGTTATCTTCCACATCAATCCGCGCCGCCAGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCG  >  1:1015381/1‑238                                                                                                                                 CTGGTAAATTCAATTTTGCCTGAATGCTGATCAATCAAATTACGAGCGATGGATAAGCCAAGCCCGGTGCCACCTTCGCGGCCGCTGACCATCGGGTAAAACAGCGTATCCTGCAAATGAGGCGGAATGCCCGGCCCGTTATCTTCCACATCAATCCGCGCCGCCAGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCG  >  1:400223/1‑238                                                                                                                                 CTGGTAAATTCAATTTTGCCTGAATGCTGATCAATCAAATTACGAGCGATGGATAAGCCAAGCCCGGTGCCACCTTCGCGGCCGCTGACCATCGGGTAAAACAGCGTATCCTGCAAATGAGGCGGAATGCCCGGCCCGTTATCTTCCACATCAATCCGCGCCGCCAGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCG  >  1:511465/1‑238                                                                                                                                 CTGGTAAATTCAATTTTGCCTGAATGCTGATCAATCAAATTACGAGCGATGGATAAGCCAAGCCCGGTGCCACCTTCGCGGCCGCTGACCATCGGGTAAAACAGCGTATCCTGCAAATGAGGCGGAATGCCCGGCCCGTTATCTTCCACATCAATCCGCGCCGCCAGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCG  >  2:365762/1‑238                                                                                                                                                                                                                                                                                                                                                                                GCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4055208‑4055336 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cTGGTAAATTCAATTTTGCCTGAATGCTGATCAATCAAATTACGAGCGATGGATAAGCCAAGCCCGGTGCCACCTTCGCGGCCGCTGACCATCGGGTAAAACAGCGTATCCTGCAAATGAGGCGGAATGCCCGGCCCGTTATCTTCCACATCAATCCGCGCCGCCAGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCG  >  NC_000913/4055340‑4055576

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