Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	CP009273	3,178,254	IS5 (+) +4 bp	46.7%	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.230)	5/270	NT	37.7%	noncoding (1195/1195 nt)	IS5	repeat region
	?	CP009273	3178254 =	8 (0.380)					coding (56/291 nt)	yqiC	BMFP family putative fusogenic protein
*	?	CP009273	= 2283592	NA (NA)	9 (0.420)	9/272	NT	52.2%	noncoding (1/1195 nt)	IS5	repeat region
	?	CP009273	= 3178257	8 (0.380)					coding (59/291 nt)	yqiC	BMFP family putative fusogenic protein

AAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  CP009273/2283455‑2283592 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAACCAGAATTAT  <  CP009273/3178257‑3178126                                                                                                                                                                                                                                                                                AAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGA                                                                                                                           <  1:328738/149‑1                       gtGCTTACGGCTGAATGCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGATTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATT                                                                                                     <  1:18228/147‑1                                                                      TGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTAC                                                           >  2:339072/1‑144                                                                                  TTCTATGGCTCCAGCTCACAAACATGACCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCACTGATTCGTGAACCTTGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGC                                          <  1:5659/149‑1                                                                                   gCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGA                                              <  2:189464/143‑1                                                                                           TCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCCTTGTAGT                                 >  1:171520/1‑149                                                                                                               TCATATCAGGGACGTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGT                  <  2:156600/144‑1                                                                                                                              GTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAACCAGAATTA   <  2:102664/144‑1                                                                                                                               TTCGCACCTTCCTTAGACATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCAGTGTAGTGCGATAAGGGTAAGTCATAAAACAGACTTAT  >  2:41861/1‑144                                                                                                                                                                                                                                                                                AAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  CP009273/2283455‑2283592 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAACCAGAATTAT  <  CP009273/3178257‑3178126

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