Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	NC_000913	4,542,308	IS5 (–) +4 bp	7.8%	coding (272‑275/597 nt)	fimE →	regulator for fimA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3652036 =	NA (NA)	6 (0.070)	5/176	NT	6.0%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	= 4542311	89 (1.050)					coding (275/597 nt)	fimE	regulator for fimA
*	?	NC_000913	= 3653230	NA (NA)	9 (0.100)	9/176	NT	8.8%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	4542308 =	89 (1.050)					coding (272/597 nt)	fimE	regulator for fimA

CTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3653134‑3653230 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGC  >  NC_000913/4542308‑4542399                                                                                                                                                                                               CTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCT                                                                                            >  2:978263/1‑99        CGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAACTGGA                                                                                     >  1:1107057/1‑99              TAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAACTGGAAAGGCG                                                                               >  2:1699668/1‑99                      TGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAACTGGAAAGGCGCTGACCGG                                                                       >  1:382637/1‑99                       GACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAACTGGAAAGGCGCTGACCGGA                                                                      >  1:1065627/1‑99                                  AACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAACTGGAAAGGCGCTGACCGGACTGACGCTATA                                                           <  2:312006/99‑1                                             ATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCG                                                <  1:564193/99‑1                                                                         TATCAGGGACTTGTTCGCACCTTCCCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAG                    <  2:571633/99‑1                                                                                  CTTATTCGCACCTTCCCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGC           >  2:129874/1‑99                                                                                           ACCTTCCCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGC  >  1:1208000/1‑99                                                                                                                                                                                               CTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3653134‑3653230 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGC  >  NC_000913/4542308‑4542399

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