Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,652,623	IS2 (–) +5 bp	intergenic (+2638/+101)	ydfD → / ← ynfP	lysis protein/protein YnfP

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1651548	NA (NA)	20 (0.780)	20/486	0.3	79.7%	noncoding (706/706 nt)	IS2	repeat region
	?	NC_000913	= 1652627	5 (0.200)					intergenic (+2642/+101)	ydfD/ynfP	lysis protein/protein YnfP
*	?	NC_000913	1652623 =	3 (0.120)	15 (0.590)	12/484	0.9	83.0%	intergenic (+2638/+105)	ydfD/ynfP	lysis protein/protein YnfP
	?	NC_000913	= 2070271	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

CGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGAAAAA  >  NC_000913/1651304‑1651781                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          cGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAAAGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAgggg                                                                                                                                                                                                                                     >  1:237686/1‑249 (MQ=14)       ccGGCGAGGAACGGAATAGCATAGAGCTTCGGGAAAACGGTAAGGCGTGATTACTTCAGTATCATGCCGAAACCAGCCGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCTAGCATTATAACTAATGGCCTCCGCATAGTGCGCTGGGTTATCGCGCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGTTTAAGTGATAACAGAGTTCTGGAAATATAGGTACAAATCCACAGGTCAACTAc                                                                                                                                                                                                                               <  2:119545/251‑1 (MQ=11)                                            cGATAACGCGTGACTACATCAGTATCATGCCCACACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGCATATCTGCTGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAACTCCACAGGTGAACTACGCTACTCTGTTTTTTACGCAATAAAAACTGGTGtttt                                                                                                                                                                                          >  2:190777/1‑251 (MQ=17)                                                            cATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCAGCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGTAAATATATGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAg                                                                                                                                                                          <  1:49921/251‑1 (MQ=255)                                                                aGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTc                                                                                                                                                                      >  1:156882/1‑251 (MQ=255)                                                                 gtagcATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCAGCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTcc                                                                                                                                                                     <  2:156882/247‑1 (MQ=255)                                                                        gCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCa                                                                                                                                                              <  1:30335/251‑1 (MQ=255)                                                                                                                                              ccGCATAGTGCGCTGGGTTAGCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAAGTTCTg                                                                                         <  2:210040/250‑1 (MQ=255)                                                                                                                                                         gcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTgaaagaa                                                                             <  1:223451/251‑1 (MQ=255)                                                                                                                                                                                aaTATCTGCGGCAGAGGGCTTGAAATGGGTTAAGTGATAACAGATGTCTGGCAATATAGGGGCAAATCCACAGGTGAACTACGATCCTCTGTTGGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGTGACCCGACAAGATTTAACAATAAAAGAACTGATTGAAAAAGTTTTGGCACTGAAAGAAAAGTAAGTCGCCAAAACATCACt                                                      >  2:263719/1‑251 (MQ=255)                                                                                                                                                                                 atatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTc                                                     <  1:64488/251‑1 (MQ=255)                                                                                                                                                                                                                   gATAACAGATGTCTGGAAAAATAGGGGCAAATCCACAGGTGAACTACGATCCTATGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCGCCCAATCTTGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGa                   <  2:63765/251‑1 (MQ=255)                                                                                                                                                                                                                        cAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCAt               >  1:96725/1‑250 (MQ=255)                                                                                                                                                                                                                         aGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATaa             <  1:250904/251‑1 (MQ=255)                                                                                                                                                                                                                             gTCTGGAAATAGAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTTGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCATTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCg          <  1:14419/250‑1 (MQ=255)                                                                                                                                                                                                                                   aaaTATAGGGGCCAATCCACATGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATCCCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGaaaa   <  1:125046/251‑1 (MQ=255)                                                                                                                                                                                                                                     atatAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGaaaaa  >  1:264902/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          CGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGAAAAA  >  NC_000913/1651304‑1651781

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