New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3571080	17 (0.840)	4 (0.210)	4/278	1.7	17.8%	coding (240/1974 nt)	glgX	glycogen debranching enzyme
	?	NC_000913	= 3571279	21 (1.110)					coding (41/1974 nt)	glgX	glycogen debranching enzyme
Rejected: Frequency below cutoff threshold.

ACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3570922‑3571080 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cccctcggcgCGCATTACGACGGTCAGGGCGTCAACTTCACACTTTTCTCCGCTCATGCCGAGCGGGTAGAACTGTGTGTCTTTGACGCCAATGGCCAGGAACATCGCTATGACTTGCCAGGGCACAGTGGCGACATTTGGCACGGTTATCTGC  <  NC_000913/3571279‑3571136                                                                                                                                                                                                                                                                                                                 ACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAATCCGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGC                                                                                                                                                           >  2:10718/1‑150     CCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTC                                                                                                                                                      <  2:152335/151‑1        CTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGG                                                                                                                                                    >  2:132482/1‑150                                                                                   CTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGCGCATTACGACGGTCAGGGCGTCAACTTCACACTTTTCTCCGCTCATGCCGAGCGGGTAGAACTGTGTGTCTTT                                                                        >  2:87850/1‑151                                                                                                                   GGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGCGCATTACGACGGTCAGGGCGTCAACTTCACACTTTTCTCCGCTCATGCCGAGCGGGTAGAACTGTGTGTCTTTGACGCCAATGGCCAGGAACATCGCTATGACTT                                        >  2:309515/1‑151                                                                                                                         ATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGCGCATTACGACGGTCAGGGCGTCAACTTCACACTTTTCTCCGCTCATGCCGAGCGGGTAGAACTGTGTGTCTTTGACGCCAATGGCCAGGAACATCGCTATGACTTGCCAGG                                  >  1:277706/1‑151                                                                                                                                        GGGTTAAAGCGATGCCCCTCGGCGCGCATTACGACGGTCAGGGCGTCAACTTCACACTTTTCTCCGCTCATGCCGAGCGGGTAGAACTGTGTGTCTTTGACGCCAATGGCCAGGAACATCGCTATGACTTGCCAGGGCACAGTGGCGACAT                   >  1:72730/1‑151                                                                                                                                                         CTCGGCGCGCATTACGACGGTCAGGGCGTCAACTTCACACTTTTCTCCGCTCATGCCGAGCGGGTAGAACTGTGTGTCTTTGACGCCAATGGCCAGGAACATCGCTATGACTTGCCAGGGCACAGTGGCGACATTTGGCACGGTTATCTGC  <  2:1120/151‑1                                                                                                                                                                                                                                                                                                                 ACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3570922‑3571080 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cccctcggcgCGCATTACGACGGTCAGGGCGTCAACTTCACACTTTTCTCCGCTCATGCCGAGCGGGTAGAACTGTGTGTCTTTGACGCCAATGGCCAGGAACATCGCTATGACTTGCCAGGGCACAGTGGCGACATTTGGCACGGTTATCTGC  <  NC_000913/3571279‑3571136

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