New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	508683 =	99 (0.830)	5 (0.070) +AAAAGCAATAATTTCAACAGGTTTTGTTGCCCAGCGCATCAGGCAATCTTGTTCTCAGCACATTG	4/168	8.6	6.0%	noncoding (102/236 nt)	REP48 (repetitive extragenic palindromic) element; contains 5 REP sequences	REP48 (repetitive extragenic palindromic) element; contains 5 REP sequences
	?	NC_000913	4053426 =	179 (1.500)					intergenic (+8/+221)	hemN/yshB	coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent/uncharacterized protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAATGGTGCTCGGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/508731‑508683 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAAT  >  NC_000913/4053426‑4053468                                                  |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||                                            gaatcaTGCTCAGATGCCTGTTGCTGCGTAAACGCCTTATCCGGCCTACAAAAGCAATAATTTCAACAGGTTTTGTTGCCCAGCGCATCAGGCAATCTTGTTCTCAGCACATTGATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGC         <  2:1479890/144‑1   atcaTGCTCAGATGCCTGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCAATAATTTCAACAGGTTTTGTTGCCCAGCGCATCAGGCAATCTTGTTCTCAGCACATTGATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCA       <  1:541101/146‑1     caTGCTCAGATGCCTGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCAATAATTTCAACAGGTTTTGTTGCCCAGCGCATCAGGCAATCTTGTTCTCAGCACATTGATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATC     <  1:1158709/148‑1     caTGCTCAGATGCCTGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCAATAATTTCAACAGGTTTTGTTGCCCAGCGCATCAGGCAATCTTGTTCTCAGCACATTGATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATC     <  1:659457/148‑1        GCTCAGATGCCTGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCAATAATTTCAACAGGTTTTGTTGCCCAGCGCATCAGGCAATCTTGTTCTCAGCACATTGATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAAT  <  1:1867149/150‑1                                                  |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||                                            CAATGGTGCTCGGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/508731‑508683 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAAT  >  NC_000913/4053426‑4053468

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