Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	NC_000913	1,292,151	IS4 (+) +12 bp	7.9%	coding (695‑706/909 nt)	galU →	UTP‑‑glucose‑1‑phosphate uridylyltransferase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1292151 =	76 (0.960)	8 (0.090)	8/174	NT	8.4%	coding (695/909 nt)	galU	UTP‑‑glucose‑1‑phosphate uridylyltransferase
	?	NC_000913	4502090 =	NA (NA)					noncoding (1426/1426 nt)	IS4	repeat region
*	?	NC_000913	= 1292162	76 (0.960)	5 (0.050)	5/176	NT	5.4%	coding (706/909 nt)	galU	UTP‑‑glucose‑1‑phosphate uridylyltransferase
	?	NC_000913	= 4503515	NA (NA)					noncoding (1/1426 nt)	IS4	repeat region

TGTTGCGCCGTCTAATCTCGCTATTGTGGGTCGTTACGTACTTAGCGCGGATATTTGGCCGTTGCTGGCAAAAACCCCTCCGGGAGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1292074‑1292162 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TAATGCCGATCAGTTAAGGATCAGTTGACCGATCCAGTGGCTGTGTAAGAATCCGGAAACGCTCACTTGTTTCCGGATTTT  <  NC_000913/4503515‑4503435                                                                                                                                                                            TGTTGCGCCGTCTAATCTCGCTATTGTGGGTCGTTACGTACTTAGCGCGGATATTTGGCCGTTGCTGGCAAAAACCCCTCCGGGAGCTGTAATGCCGAT                                                                         <  2:1026774/99‑1 TGTTGCGCCGTCTAATCTCGCTATTGTGGGTCGTTACGTACTTAGCGCGGATATTTGGCCGTTGCTGGCAAAAACCCCTCCGGGAGCTGTAATGCCGAT                                                                         >  2:1079246/1‑99                                                          GCCGTTGCTGGCAAAAACCCCTCCGGGAGCTGTAATGCCGATCAGTTAAGGATCAGTTGACCGATCCAGTGGCTGTGTAAGAATCCGGAAACGCTCACT                <  2:1898349/99‑1                                                           CCGTTGCTGGCAAAAACCCCTCCGGGAGCTGTAATGCCGATCAGTTAAGGATCAGTTGACCGATCCAGTGGCTGTGTAAGAATCCGGAAACGCTCACTT               <  2:876871/99‑1                                                                        AAACCCCTCCGGGAGCTGTAATGCCGATCAGTTAAGGATCAGTTGACCGATCCAGTGGCTGTGTAAGAATCCGGAAACGCTCACTTGTTTCCGGATTTT  >  1:550201/1‑99                                                                                                                                                                            TGTTGCGCCGTCTAATCTCGCTATTGTGGGTCGTTACGTACTTAGCGCGGATATTTGGCCGTTGCTGGCAAAAACCCCTCCGGGAGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1292074‑1292162 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TAATGCCGATCAGTTAAGGATCAGTTGACCGATCCAGTGGCTGTGTAAGAATCCGGAAACGCTCACTTGTTTCCGGATTTT  <  NC_000913/4503515‑4503435

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