New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4214492	63 (1.740)	5 (0.150)	5/224	2.6	7.4%	coding (213/930 nt)	metA	homoserine O‑transsuccinylase
	?	NC_000913	= 4214509	67 (1.980)					coding (230/930 nt)	metA	homoserine O‑transsuccinylase
Rejected: Frequency below/above cutoff threshold.

TTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4214340‑4214492 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgttgcgcGATTCACGGGAATCGATGCGCAACAGCTGAATATCGACCTGCAAAGGT  <  NC_000913/4214509‑4214462                                                                                                                                                                                                           TTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTT                                                      >  1:66214/1‑149     TGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCG                                                   >  5:129758/1‑148      GATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGC                                                  >  6:85531/1‑148        TGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCGAT                                               >  6:123892/1‑149                              TCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCGATTCACGGGAATCGATGCGCAACA                         <  5:123892/149‑1                                   AAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCGATTCACGGGAATCGATGCGCAACAGCTGA                    <  5:85531/149‑1                                                                                                TCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCGATTCACGGGAATCGATGCGCAACAGCTGAATATCGACCTGCAAAGGT  >  5:169354/1‑106                                                                                                TCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCGATTCACGGGAATCGATGCGCAACAGCTGAATATCGACCTGCAAAGGT  <  6:169354/106‑1                                                                                                                                                                                                           TTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4214340‑4214492 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgttgcgcGATTCACGGGAATCGATGCGCAACAGCTGAATATCGACCTGCAAAGGT  <  NC_000913/4214509‑4214462

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