New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	3375666 =	19 (1.040)	4 (0.230)	4/272	1.9	18.6%	coding (798/1368 nt)	dcuD	putative transporter
	?	NC_000913_3_vch_pgi	3375686 =	17 (0.980)					coding (818/1368 nt)	dcuD	putative transporter
Rejected: Frequency below/above cutoff threshold.

GTTACCACGTTGGCAAACTGCGTACCCATGCCGTCAAAAAACGCCTGCACATCGTCCATTGTTTCGCGCAAGTTATGCTTGCGGAAGAACTCAACAAACATCGTCACAGTCAAACTCAGTAACATCACCACCACCAGATGCAGTTCTGACTGCATC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_vch_pgi/3375821‑3375666 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑actgcatcTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAGTTTGCCAACGTGGTAACGCTGGTGGTC  >  NC_000913_3_vch_pgi/3375686‑3375831                                                                                                                                                                                                                                                                                                                GTTACCACGTTGGCAAACTGCGTACCCATGCCGTCAAAAAACGCCTGCACATCGTCCATTGTTTCGCGCAAGTTATGCTTGCGGAAGAACTCAACAAACATCGTCACAGTCAAACTCAGTAACATCACCACCACCAGATTCAGTTCTGA                                                                                                                                                           >  1:13159/1‑149  TTACCACGTTGGCAAACTGCGTACCCATGCCGTCAAAAAACGCCTGCACATCGTCCATTGTTTCGCGCAAGTTATGCTTGCGGAAGAACTCAACAAACATCGTCACAGCCAAACTCAGTAACATCACCACCACCAAATGCAGTTCTGAC                                                                                                                                                          >  2:72644/1‑149                                                                                                                                    CACCAGATGCAGTTCTGACTGCATCTGGTGGTGGTGATGTTACTGAGT                                                                                                                             <  1:25605/48‑1                                                                                                                                    CACCAGATGCAGTTCTGACTGCATCTGGTGGTGGTGATGTTACTGAGT                                                                                                                             >  2:25605/1‑48                                                                                                                                      CCAGATGCAGTTCTGACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTT                                               >  1:243026/1‑124                                                                                                                                      CCAGATGCAGTTCTGACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTT                                               <  2:243026/124‑1                                                                                                                                                     ACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTCAGTTCTTCCGCAAGCATTACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAGTTTGCCAACGTGGTAACGCTGG       <  1:299042/149‑1                                                                                                                                                          ATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTCAGCAGTTTGCCCACGTGTTAACGCTGtgtgtc  >  2:48005/1‑143                                                                                                                                                                                                                                                                                                                GTTACCACGTTGGCAAACTGCGTACCCATGCCGTCAAAAAACGCCTGCACATCGTCCATTGTTTCGCGCAAGTTATGCTTGCGGAAGAACTCAACAAACATCGTCACAGTCAAACTCAGTAACATCACCACCACCAGATGCAGTTCTGACTGCATC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_vch_pgi/3375821‑3375666 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑actgcatcTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAGTTTGCCAACGTGGTAACGCTGGTGGTC  >  NC_000913_3_vch_pgi/3375686‑3375831

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