New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4617850 =	109 (1.200)	5 (0.060)	4/268	NT	4.4%	coding (528/780 nt)	deoC	2‑deoxyribose‑5‑phosphate aldolase, NAD(P)‑linked
	?	NC_000913	4617871 =	113 (1.290)					coding (549/780 nt)	deoC	2‑deoxyribose‑5‑phosphate aldolase, NAD(P)‑linked
Rejected: Frequency below/above cutoff threshold.

AGTCAGCACCGAACAGTTCATCTGCAATGGCGAGATATTTCTGCGCATCTTCCGCAGTACGCACGCCGCCCGCCGGTTTGAAACCAACGGTTTTTTCTACGCCCATATCACGGATCACTTCCATCATGATGCGCGCGCTTTCCGGCGTCGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4618001‑4617850 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgcgCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCAC  >  NC_000913/4617871‑4618018                                                                                                                                                                                                                                                                                                              AGTCAGCACCGAACAGTTCATCTGCAATGGCGAGATATTTCTGCGCATCTTCCGCAGTACGCACGCCGCCCGCCGGTTTGAAACCAACGGTTTTTTCTACGCCCATATCACGGATCACTTCCATCATGATGCGCGCGCTTTCCGGCGTC                                                                                                                                                         <  2:936842/149‑1   TCAGCACCGAACAGTTCATCTGCAATGGCGAGATATTTCTGCGCATCTTCCGCAGTACGCACGCCGCCCGCCGGTTTGAAACCAACGGTTTTTTCTACGCCCATCTCACGGATCACTTCGATCATGATGCGCGCGCTTTCCGGCGTCGC                                                                                                                                                       >  1:252765/1‑149   TCAGCACCGAACAGTTCATCTGCAATGGCGAGATATTTCTGCGCATCTTCCGCAGTACGCACGCCGCCCGCCGGTTTGAAACCAACGGTTTTTTCTACGCCCATATCACGGATCACTTCCATCATGATGCGCGCGCTTTCCGGCGTCGC                                                                                                                                                       <  1:283080/149‑1    CAGCACCGAACAGTTCATCTGCAATGGCGAGATATTTCTGCGCATCTTCCGCAGTACGCACGCCGCCCGCCGGTTTGAAACCAACGGTTTTTTCTACGCCCATATCACGGATCACTTCCATCATGATGCGCGCGCTTTCCGGCGTCGCG                                                                                                                                                      <  1:269992/149‑1                                                                                          GTTTTTTCTACGCCCATATCACGGATCACTTCCATCATGATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACAGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATG                                                                >  2:728054/1‑149                                                                                                                      CTTCCATCATGATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGCTGA                                    >  1:1244961/1‑149                                                                                                                           ATCATGATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGT                               >  1:646784/1‑149                                                                                                                                GATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGCACTGTTCGGT                          >  1:834581/1‑149                                                                                                                                GATGCGCGCGCTTTCCGGCGTCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGT                          >  1:1489625/1‑149                                                                                                                                                       CGCATCATGATGGAAGTGATCCGTGATATTGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGTGCAGATGCGCGTCA   <  2:963042/149‑1                                                                                                                                                        GCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCAC  <  1:1144111/149‑1                                                                                                                                                                                                                                                                                                              AGTCAGCACCGAACAGTTCATCTGCAATGGCGAGATATTTCTGCGCATCTTCCGCAGTACGCACGCCGCCCGCCGGTTTGAAACCAACGGTTTTTTCTACGCCCATATCACGGATCACTTCCATCATGATGCGCGCGCTTTCCGGCGTCGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4618001‑4617850 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgcgCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCAC  >  NC_000913/4617871‑4618018

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.