Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	CP009273	3,178,254	IS5 (+) +4 bp	33.3%	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)	6 (0.290)	6/276	NT	26.7%	noncoding (1195/1195 nt)	IS5	repeat region
	?	CP009273	3178254 =	16 (0.810)					coding (56/291 nt)	yqiC	BMFP family putative fusogenic protein
*	?	CP009273	= 2283592	NA (NA)	10 (0.490)	9/278	NT	37.8%	noncoding (1/1195 nt)	IS5	repeat region
	?	CP009273	= 3178257	16 (0.800)					coding (59/291 nt)	yqiC	BMFP family putative fusogenic protein

GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  CP009273/2283448‑2283592 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAACCAGAATTATGTG  <  CP009273/3178257‑3178123                                                                                                                                                                                                                                                                                          GACAAGAAAGTCTCTTTTCTGGTCTGCCGGCGCTTACTGCTGAATTCACTCTCGGCGAAGGTCAGTTGATGACTCATTATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGG                                                                                                                                    <  1:35412/150‑1              CTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGA                                                                                                                       >  2:123749/1‑150                                             TTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCG                                                                                        <  2:271756/150‑1                                                     CGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGACCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGGGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGTTCAATC                                                                                >  2:10696/1‑150                                                                 GTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCC                                                                    >  1:155775/1‑150                                                                          TCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGA                                                           <  1:24148/150‑1                                                                                     CCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAG                                                <  1:357432/150‑1                                                                                                               CATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAG                      <  1:326709/150‑1                                                                                                                      TCATATCAGGGGCTTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAA               <  1:225717/150‑1                                                                                                                           TCAGGGACTTGTTCGCACCTTCCTTAGGCATAGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTAGGGACACTCATTGTAGTTTACCTGTACTGATAGGTGTTGAGAGCCCTCGGAGTGCGATAAGGGCAAGcaatagaatagg          >  1:56536/1‑138                                                                                                                                   TTGTTCGCACCTTCCTTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAACCAGAATTATGTG  <  2:75949/150‑1                                                                                                                                                                                                                                                                                          GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  CP009273/2283448‑2283592 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAGGCATTGATTCGTGAACCTGGCGAGCGATTTGCTCAATTTTTTTCGGGTCAATCATTGTAGTTTCCCTGTACTGATAGGTGTTGAGTGCCATTGTAGTGCGATAAGGGTAAGTCATAAACCAGAATTATGTG  <  CP009273/3178257‑3178123

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.