New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	951204 =	31 (1.150)	5 (0.190)	5/278	2.9	13.4%	intergenic (‑124/+68)	pflA/pflB	pyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase/formate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1
	?	NC_000913_3_pae_pgi	951236 =	34 (1.280)					intergenic (‑156/+36)	pflA/pflB	pyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase/formate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1
Rejected: Frequency below/above cutoff threshold.

CACTCAATCTATGTAATTAGATTTGACTGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAGTGGGGCCTTTTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_pae_pgi/951287‑951204 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGC  >  NC_000913_3_pae_pgi/951236‑951351                                                                                                                                                                                                          CACTCAATCTATGTAATTAGATTTGACTGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAATGGGGCCTTTTTTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAATTTGATTGAAGGTAAGAGTAA                                                     >  2:426458/1‑149                        TGACTGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAGTGGGGCCTTTTTTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTA                              <  1:426458/149‑1                                                    CAATAAAGGCTCCACGAAAGTGGGGCCTTTTTTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTTCTGCTGTTCTTTCGTGAGCGAGTTGAAGCGTACTGCCTAGC  >  1:134435/1‑149                                                              TCCACGAAAGTGGGGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGA                                                                     <  1:78044/72‑1                                                              TCCACGAAAGTGGGGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGA                                                                     >  2:78044/1‑72                                                                                                                                                                                                          CACTCAATCTATGTAATTAGATTTGACTGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAGTGGGGCCTTTTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_pae_pgi/951287‑951204 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGC  >  NC_000913_3_pae_pgi/951236‑951351

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