Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	CP009273	3,178,254	IS5 (+) +4 bp	40.0%	coding (56‑59/291 nt)	yqiC →	BMFP family putative fusogenic protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	2282398 =	NA (NA)	9 (0.340)	9/274	NT	42.1%	noncoding (1195/1195 nt)	IS5	repeat region
	?	CP009273	3178254 =	12 (0.470)					coding (56/291 nt)	yqiC	BMFP family putative fusogenic protein
*	?	CP009273	= 2283592	NA (NA)	7 (0.260)	7/276	NT	36.2%	noncoding (1/1195 nt)	IS5	repeat region
	?	CP009273	= 3178257	12 (0.460)					coding (59/291 nt)	yqiC	BMFP family putative fusogenic protein

ACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  CP009273/2283449‑2283592 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAACCAG  <  CP009273/3178257‑3178132                                                                                                                                                                                                                                                                                ACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACCCATGATGGACCCTGTGCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGG                                                                                                                           <  1:205119/149‑1         ATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGACTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGggacgcta                                                                                                                   >  2:406388/1‑141                TTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAAC                                                                                                            <  1:466689/149‑1                                              CACTGCGGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAACCATGATCTGAGATCAGGGAATTATTTGCACGTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGGCTTGCTCAATATTTTTCGG                                                                              <  1:343355/149‑1                                                               AGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTGCCTTATGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTGTTTCGGGTCAATCATTGTAGTTT                                                             <  1:245800/149‑1                                                                     TGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCCCCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGT                                                       >  1:68787/1‑149                                                                                   CCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTG                                         >  1:419175/1‑149                                                                                                                        TATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAACC    <  1:158771/149‑1                                                                                                                          TCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGATGAGTGCCATTGTAGTGCGATAAGGGTAAGTGATAgcacag  >  1:97900/1‑143                                                                                                                                                                                                                                                                                ACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  CP009273/2283449‑2283592 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAACCAG  <  CP009273/3178257‑3178132

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