Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	NC_000913	3,802,542	IS5 (–) +4 bp	8.7%	coding (514‑517/1020 nt)	waaO ←	UDP‑D‑glucose:(glucosyl)LPS alpha‑1,3‑glucosyltransferase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3652036 =	NA (NA)	11 (0.120)	10/176	NT	10.4%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	= 3802545	90 (1.030)					coding (514/1020 nt)	waaO	UDP‑D‑glucose:(glucosyl)LPS alpha‑1,3‑glucosyltransferase
*	?	NC_000913	= 3653230	NA (NA)	6 (0.070)	6/174	NT	6.0%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	3802542 =	89 (1.030)					coding (517/1020 nt)	waaO	UDP‑D‑glucose:(glucosyl)LPS alpha‑1,3‑glucosyltransferase

GAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3653142‑3653230 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cCTAACGAATGGGCGCGTTTTTCCCACCAGTCAGCTTGCCCTTCTGTAACAACCATTGCGACTTTATCGTCAGGAAATGAGAAATTAATTAGTGGCTCA  >  NC_000913/3802542‑3802639                                                                                                                                                                                             GAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAACGAATG                                                                                          <  1:1546484/99‑1             ATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAACGAATGGGCGCGTTTTTC                                                                              >  2:1519600/1‑99                                  TCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAACGAATGGGCGCGTTTTTCCCACCAGTCAGCTTGCCCTTC                                                         <  2:1704263/99‑1                                                TGACAAACATGATCTCAGATCAGGGACTTGTTCGCACCTTCCCTAACGAATGGGCGCGTTTTTCCCACCAGTCAGCTTGCCCTTCTGTAACAACCATTG                                           >  2:2208162/1‑99                                                           ATCTCATATCAGGGACTTGTTCGCACCTTCCCTAACGAATGGGCGCGTTTTTCCCACCAGTCAGCTTGCCCTTCTGTAACAACCATCGCGACTTTATCG                                >  2:878915/1‑99                                                                            TGTTCGCACCTTCCCTAACGAATGGGCGCGTTTTTCCCACCAGTCAGCTTGCCCTTCTGTAACAACCATTGCGACTTTATCGTCAGGAAATGAGAAATT               >  1:1470099/1‑99                                                                                         CCTAACGAATGGGCGCGTTTTTCCCACCAGTCAGCTTGCCCTTCTGTAACAACCATTGCGACTTTATCGTCAGGAAATGAGAAATTAATTAGTGGCTCA  >  1:1574431/1‑99                                                                                                                                                                                             GAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3653142‑3653230 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cCTAACGAATGGGCGCGTTTTTCCCACCAGTCAGCTTGCCCTTCTGTAACAACCATTGCGACTTTATCGTCAGGAAATGAGAAATTAATTAGTGGCTCA  >  NC_000913/3802542‑3802639

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