New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3928162	80 (0.680)	3 (0.030)	3/250	NT	3.8%	coding (1442/1452 nt)	viaA	stimulator of RavA ATPase activity; von Willebrand factor domain protein
	?	NC_000913	= 3928174	75 (0.670)					coding (1430/1452 nt)	viaA	stimulator of RavA ATPase activity; von Willebrand factor domain protein
Rejected: Frequency below/above cutoff threshold.

CTTCTCTGCTGTAATAATTTATCGCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3928134‑3928162 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgccTGC‑TCAGACGCTGGCGGCGATAAATTATTACAGCAGAGAAGGGACGCTCTCGCGAACAGCAGCTGGCCATACTCCACACTGAACCTGGCCGATATGCGGCAGTTGCAGCAGCAGCATAGTCAAACGAGACTGGCCGATACCGCCGC  <  NC_000913/3928174‑3928029                                                                                                                                                                                  CTTCTCTGCTGTAATAATTTATCGCCGCCAGCGTCTGA‑GCAGGCGGCGATAAATTATTACAGCAGAGAAGGGACGCTCTCGCGAACAGCAGCTGG                                                                                  <  3:777634/95‑1 CTTCTCTGCTGTAATAATTTATCGCCGCCAGCGTCTGA‑GCAGGCGGCGATAAATTATTACAGCAGAGAAGGGACGCTCTCGCGAACAGCAGCTGG                                                                                  >  4:777634/1‑95        GCTGTAATAATTTATCGCCGCCAGCGTCTGA‑GCAGGCGGCGATAAATTATTACAGCAGAGAAGGGACGCTCTCGCGAACAGCAGCTGGCCATACTCCACACTGAACCTGGCCGATATGCGGCAGTTGCAGCAGCAGCATAGTCAAACGA                     >  1:981952/1‑149                           GCCTGC‑TCAGACGCTGGCGGCGATAAATTATTACAGCAGAGAAGGGAAGCTCTCGCGAACAGCAGCTGGCCATACTCCACACTGAACCTGGCCGATATGCGGCAGTTGCAGCAGCAGCATAGTCAAACGAGAATGGCCGATACCGCCGC  >  4:853478/1‑149                                                                                                                                                                                  CTTCTCTGCTGTAATAATTTATCGCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3928134‑3928162 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgccTGC‑TCAGACGCTGGCGGCGATAAATTATTACAGCAGAGAAGGGACGCTCTCGCGAACAGCAGCTGGCCATACTCCACACTGAACCTGGCCGATATGCGGCAGTTGCAGCAGCAGCATAGTCAAACGAGACTGGCCGATACCGCCGC  <  NC_000913/3928174‑3928029

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