New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	3645087 =	17 (0.930)	4 (0.230)	4/274	1.9	19.9%	coding (96/2043 nt)	prlC	oligopeptidase A
	?	NC_000913_3_vch_pgi	3645120 =	16 (0.920)					coding (63/2043 nt)	prlC	oligopeptidase A
Rejected: Frequency below/above cutoff threshold.

CATTATGTAAAGCGCCTGTTGAGCGCTTCCTTAACCTCTTTAACCAGGACTGCGCGAATGACGAATCCGTTACTGACTCCCTTTGAATTGCCTCCGTTTTCTAAAATTCTCCCGGAACATGTCGTTCCAGCCGTGACTAAGGCATTGAACGAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_vch_pgi/3645239‑3645087 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaacgacATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAG  >  NC_000913_3_vch_pgi/3645120‑3645204                                                                                                                                                                                                                                                CATTATGTAAAGCGCCTGTTGAGCGCTTCCTTAACCTCTTTAACCAGGACTGCGCGAATGACGAATCCGTGACTGACTCCCTTTGAATTGCCTCCGTTTTCTAAAATTCTCCCGGAACATGTCGTTCCAGCCGTGACTAAGCAATTGAA                                                                                           >  1:108518/1‑149  ATTATGTAAAGCGCCTGTTGAGCGCTTCCTTAACCTCTTTAACCAGGACTGCGCGAATGACGAATCCGTTACTGACTCCCTTTGAATTGCCTCCGTTTTCTAAAATTCTCCCGGAACATGTCGTTCCAGCCGTGACTAAGGCATTGAAC                                                                                          <  1:4493/149‑1                                                                                                              TCCCGGAACATGTCGTTCCAGCCGTGACTAAGGCATTGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAG  >  1:21726/1‑129                                                                                                              TCCCGGAACATGTCGTTCCAGCCGTGACTAAGGCATTGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAG  <  2:21726/129‑1                                                                                                                        TGTCGTTCCAGCCGTGACTAAGGCATTGAACGACA                                                                                      >  1:48128/1‑35                                                                                                                        TGTCGTTCCAGCCGTGACTAAGGCATTGAACGACA                                                                                      <  2:48128/35‑1                                                                                                                                                                                                                                                CATTATGTAAAGCGCCTGTTGAGCGCTTCCTTAACCTCTTTAACCAGGACTGCGCGAATGACGAATCCGTTACTGACTCCCTTTGAATTGCCTCCGTTTTCTAAAATTCTCCCGGAACATGTCGTTCCAGCCGTGACTAAGGCATTGAACGAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_vch_pgi/3645239‑3645087 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaacgacATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAG  >  NC_000913_3_vch_pgi/3645120‑3645204

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