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

TGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGT  >  CP009273/2282264‑2282531                                                                                                                                       |                                                                                                                                      tGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGCTATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATAGAAGGAAGGTGCg                                                                                                                              >  1:55157/1‑144 (MQ=255)             ccAGGATTGATTTTCTTATTAGTTGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCTTCAACTTTCTCACCATACTCAACAGTTGATTCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGa                                                                                                                  <  1:111372/144‑1 (MQ=255)                ggATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGaaac                                                                                                               >  1:138292/1‑143 (MQ=255)                                                    aGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGa                                                                      >  2:271410/1‑149 (MQ=255)                                                                     cAGGACAACACCATCCACTGTCTCACCATACTCAACAATAAATGCCAATGCTATGGGTGAAAATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTGCATGTTTGAGCCGAGTTTGTCTCCCg                                                     >  2:2356/1‑149 (MQ=21)                                                                       ggACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtttc                                                        >  1:174642/1‑144 (MQ=255)                                                                                  tCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGGTTGAGCCGACTTTTTCTCCCGTCACTGCCTTGaa                                        >  2:130354/1‑149 (MQ=255)                                                                                                      aaCAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGTTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCAATTTAGACCGttt                    >  2:41003/1‑149 (MQ=255)                                                                                                        cAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCg                       >  1:24620/1‑144 (MQ=255)                                                                                                             gATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGCCCGTttct                  >  1:9035/1‑144 (MQ=255)                                                                                                                gCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGGGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCg               <  1:116861/144‑1 (MQ=255)                                                                                                                      gCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGTTTCCTTCCACGGTTAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGc    <  2:259216/149‑1 (MQ=14)                                                                                                                        tATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTATTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGt  >  2:32078/1‑149 (MQ=255)                                                                                                                                       |                                                                                                                                      TGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGT  >  CP009273/2282264‑2282531

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.