New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3757352	209 (3.160)	3 (0.050)	3/264	NT	1.5%	coding (476/1152 nt)	yiaY	L‑threonine dehydrogenase
	?	NC_000913	= 3757393	194 (3.000)					coding (435/1152 nt)	yiaY	L‑threonine dehydrogenase
Rejected: Frequency below/above cutoff threshold.

GTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTCGGCATACCAATCATCAGAGAGGAGTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3757202‑3757352 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgaAATGACCCGTTTCTGCATCATCACTGACGAAGCGCGTCATATCAAAATGGCGATTGTTGATAAACATGTCACTCCGCTGCTTTCTGTCAATGACTCCTCTCTGATGATTGGTATGCCGAAGTCACTGACCGCCGCAACGGGTATGGAT  <  NC_000913/3757393‑3757246                                                                                                                                                                                                                                                                                                             GTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTCGGCATACCAATCATCAGAGAGGAGTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTT                                                                                                                                                        <  2:237854/149‑1   GGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTCGGCATACCAATCATCAGAGAGGAGTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGA                                                                                                                                                      >  1:517552/1‑149   GGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTCGGCATACCAATCATCAGAGAGGAGTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGA                                                                                                                                                      <  2:865/149‑1                                        AAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTCGGCATACCAATCATCAGAGAGGAGTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGATATGACCCGTTTCTGCATCATCACTGACGAAGCGCGT                                                                                                                 <  2:528986/149‑1                                                                                            GAGAGGAGTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGAAATGACCCGTTTCTGCATCATCACTGACGAAGCGCGTCATATCAAAATG                                                                                                     <  3:502250/109‑1                                                                                            GAGAGGAGTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGAAATGACCCGTTTCTGCATCATCACTGACGAAGCGCGTCATATCAAAATG                                                                                                     >  4:502250/1‑109                                                                                                                                                     TGAAATGACCCGTTTCTGCATCATCACTGACGAAGCGCGTCATATCAAAATGGCGATTGTTGATAAACATGTCACTCCGCTGCTTTCTGTCAATGACTCCTCTCTGATGATTGGTATGCCGAAGTCACTGACCGCCGCAACGGGTATGG    <  4:180573/149‑1                                                                                                                                                       AAATGACCCGTTTCTGCATCATCACTGACGAAGCGCGTCATATCAAAATGGCGATTGTTGATAAACATGTCACTCCGCTGCTTTCTGTCAATGACTCCTCTCTGATGATTGGTATGCCGAAGTCACTGACCGCCGCAACGGGTATGGAT  >  2:235844/1‑149                                                                                                                                                       AAATGACCCGTTTCTGCATCATCACTGACGAAGCGCGTCATATCAAAATGGCGATTGTTGATAAACATGTCACTCCGCTGCTTTCTGTCAATGACTCCTCTCTGATGATTGGTATGCCGAAGTCACTGACCGCCGCAACGGGTATGGAT  <  4:105846/149‑1                                                                                                                                                                                                                                                                                                             GTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTCGGCATACCAATCATCAGAGAGGAGTCATTGACAGAAAGCAGCGGAGTGACATGTTTATCAACAATCGCCATTTTGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3757202‑3757352 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgaAATGACCCGTTTCTGCATCATCACTGACGAAGCGCGTCATATCAAAATGGCGATTGTTGATAAACATGTCACTCCGCTGCTTTCTGTCAATGACTCCTCTCTGATGATTGGTATGCCGAAGTCACTGACCGCCGCAACGGGTATGGAT  <  NC_000913/3757393‑3757246

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