New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2185162	41 (1.190)	3 (0.090)	3/262	4.1	7.2%	coding (140/789 nt)	thiM	hydoxyethylthiazole kinase
	?	NC_000913	= 2185172	38 (1.170)					coding (130/789 nt)	thiM	hydoxyethylthiazole kinase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2185008‑2185162 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gctggcgCTCGGTGCATCGCCAGCGATGGTTATCGAAACCGAAGAGGCCAGTCAGTTTGCGGCTATCGCCAGTGCCTTGTTGATTAACGTTGGCACACTGACGCAGCCACGCGCTCAGGCGATGCGTGCTGCCGTTGAGCAAGCAAAAAGCTCT  <  NC_000913/2185172‑2185026                                                                                                                                                                                                                                                                                                                AAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCG                                                                                                                                                           <  6:498/149‑1     GTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAATCTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGG                                                                                                                                                       >  7:72064/1‑149                           TTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGCTCGGTGCATCGCCAGCGA                                                                                                                                  <  3:23341/148‑1                                                                                           CAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGCTCGGTGCATCGCCAGCGATGGTTATCGAAACCGAAGAGGCCAGTCAGTTT                                                                                                  >  5:16351/1‑116                                                                                           CAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGCTCGGTGCATCGCCAGCGATGGTTATCGAAACCGAAGAGGCCAGTCAGTTT                                                                                                  <  6:16351/116‑1                                                                                                                                                      CTGGCGCTCGGTGCATCGCCAGCGATGGTTATCGAAACCGAAGAGGCCAGTCAGTTTGCGGCTATCGCCAGTGCCTTGTTGATTAACGTTGGCACACTGACGCAGCCACGCGCTCAGGCGATGCGTGCTGCCGTTGAGCAAGCAAAAA       <  3:40985/148‑1                                                                                                                                                          CGCTCGATGCATCGCCAGCGATGGTTATCGAAACCGAAGAGGCCAGTCAGTTTGCGGCTATCGCCAGTGCCTTGTTGATTAACGTTGGCACACTGACGCAGCCACTCGCTCAGGCGATGCGTGCTGCAGTTGAGCAAGCAAAAAGCTCT  >  7:106115/1‑149                                                                                                                                                                                                                                                                                                                AAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2185008‑2185162 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gctggcgCTCGGTGCATCGCCAGCGATGGTTATCGAAACCGAAGAGGCCAGTCAGTTTGCGGCTATCGCCAGTGCCTTGTTGATTAACGTTGGCACACTGACGCAGCCACGCGCTCAGGCGATGCGTGCTGCCGTTGAGCAAGCAAAAAGCTCT  <  NC_000913/2185172‑2185026

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