New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2333629	15 (0.490)	4 (0.150)	4/112	2.1	18.7%	pseudogene (774/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
	?	NC_000913	= 2334464	22 (0.850)					coding (496/624 nt)	yfaT	DUF1175 family protein, putative host defense protein
Rejected: Frequency below/above cutoff threshold.

GCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2333556‑2333629 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aagcgccacgAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCAACAACTTATGACATGGAAGGACACCCGA  <  NC_000913/2334464‑2334402                                                                                                                                           GCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGC                                                                       <  2:356760/68‑1      GCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCAC                                                                  >  1:201495/1‑68       CTTCACCCTTTTTTTTACCCGCTGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCTCG                                                                 <  2:552555/68‑1                                   CACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGAAAACTGCCAACGGCATGCGCGCAGTCA                                     >  1:549548/1‑68                                     CCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGAAAACTGACAACGGAATGCGCGCA                                         >  1:637356/1‑62                                     CCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGAAAACTGACAACGGAATGCGCGCA                                         <  2:637355/62‑1                                          AGCTCTTTGCCTGACACTTTGCTAAGCGCCACGAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCg                              >  2:392541/1‑67                                                                 AAGCGCCACGAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCAACAACTTATGACATGGAAGGACA       <  1:658663/68‑1                                                                  AGCGCCACGAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCAACAACTTATGACATGGAAGGACAC      >  1:439592/1‑68                                                                      CCACGAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCAACAACTTATGACATGGAAGGACACCCGA  <  1:1152030/68‑1                                                                      CCACGAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCAACAACTTATGACATGGAAGGACACCCGA  <  1:1153544/68‑1                                                                                                                                           GCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2333556‑2333629 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aagcgccacgAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCAACAACTTATGACATGGAAGGACACCCGA  <  NC_000913/2334464‑2334402

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