Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	CP009273	3,178,254	IS5 (+) +4 bp	55.6%	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)	5 (0.280)	5/268	NT	44.7%	noncoding (1195/1195 nt)	IS5	repeat region
	?	CP009273	3178254 =	6 (0.350)					coding (56/291 nt)	yqiC	BMFP family putative fusogenic protein
*	?	CP009273	= 2283592	NA (NA)	10 (0.560)	10/270	NT	61.8%	noncoding (1/1195 nt)	IS5	repeat region
	?	CP009273	= 3178257	6 (0.340)					coding (59/291 nt)	yqiC	BMFP family putative fusogenic protein

CTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  CP009273/2283487‑2283592 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAACCAGAATTATGTGAAGCTAT  <  CP009273/3178257‑3178116                                                                                                                                                                                                                                                          CTGAATTCACTGTCGGCGAAGGTAAGATGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGAAAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGCTATGCTCAATTT                                                                                                     >  1:313338/1‑149         ACTTTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGG                                                                                             <  1:331259/149‑1                        AAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATT                                                                                    >  2:206725/1‑143                                   TCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGAGTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGgtctg                                                                   >  1:109081/1‑144                                            ACCCTGTCCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTT                                                          <  1:36913/149‑1                                                  TTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGC                                                    >  1:288055/1‑149                                                     TATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGTAGGCATTGATTGGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTCCTGATAGGTGTGGAG                                                       >  2:26550/1‑143                                                                  GACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGCCAATCATTGTCGTTTACCTGTACCGATAGTTGTTGAGTGCCATTGTAGTG                                          <  2:141481/143‑1                                                                   ACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAG                                   <  1:329029/149‑1                                                                                                    ACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAACCAGAATTATGTGAAGCTAT  <  1:280299/149‑1                                                                                                                                                                                                                                                          CTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  CP009273/2283487‑2283592 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAACCAGAATTATGTGAAGCTAT  <  CP009273/3178257‑3178116

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.