Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	293,207	IS2 (–) +5 bp	intergenic (‑259/+9)	yagK ← / ← yagL	uncharacterized protein YagK/resolvase‑like catalytic domain‑containing protein YagL

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	293207 =	0 (0.000)	18 (0.680)	17/488	0.5	100%	intergenic (‑259/+13)	yagK/yagL	uncharacterized protein YagK/resolvase‑like catalytic domain‑containing protein YagL
	?	NC_000913	381260 =	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region
*	?	NC_000913	= 293211	0 (0.000)	13 (0.490)	13/488	0.7	100%	intergenic (‑263/+9)	yagK/yagL	uncharacterized protein YagK/resolvase‑like catalytic domain‑containing protein YagL
	?	NC_000913	= 382590	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

GCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGG  >  NC_000913/382345‑382791                                                                                                                                                                                                                                                      |                                                                                                                                                                                                          gCGGAGTCCGGTGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGGTCGAGCATTATAACGAATGGCATCCGCATATTGCGCTGGGTTATCGGTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAAGGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGGAAATCCAATTAc                                                                                                                                                                                                      <  1:227336/251‑1 (MQ=255)     aGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTAt                                                                                                                                                                                                  >  1:79509/1‑251 (MQ=32)             gagTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGACTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGt                                                                                                                                                                                          >  1:13721/1‑251 (MQ=25)                         aGCAGAGAGCTTCGTTAAAACGATAAAGCGTGACTACATCAGTAGCATGTGCAAAACAGACGGGTTAACGGCAGCAAAGAACCTTTGAGATGCGTTCGACCATTATAACTAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGTGAATATCTGCTGCAGCGTGCCTGTAATGGTTTAAGTGATAACAGATGTCTGGTAATATAGGGGCAAATCCAATTCCCTATCAGGCAGTTTgcgtgaatc                                                                                                                                                                               <  2:1851/250‑7 (MQ=11)                             gagagCTTCGTGAACACGATAAAGCGTGCCTACATCCGTATCATGCCCAAACAAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAc                                                                                                                                                                          >  2:31413/1‑251 (MQ=255)                                agCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCgg                                                                                                                                                                       <  1:205180/251‑1 (MQ=255)                                 gCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGa                                                                                                                                                                      <  1:219821/251‑1 (MQ=255)                                           aaCGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGTTCCCTGAATaa                                                                                                                                                            <  2:40718/251‑1 (MQ=255)                                                                            cAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGa                                                                                                                           <  2:117016/251‑1 (MQ=255)                                                                                 cAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGGGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGAGGTCTGGAAATATAGGGGCAAATCCAATTACCTAGCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaaca                                                                                                                      <  1:258944/251‑1 (MQ=255)                                                                                           aaCGGCAGCAAAGAACCTTGCAGAGGCGTTCGTGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATCACGGCACCATGACAAATAAGGGTATTCGCGCCGATAACAACTTCACCTg                                                                                                            >  2:22450/1‑251 (MQ=255)                                                                                                  gcaAAGAACCTTGCAGAGGAGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAGTCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGAATTCGCGCCGATAACAACTTCACCTGCAGTTtc                                                                                                     >  2:195654/1‑251 (MQ=255)                                                                                                                       ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGc                                                                                >  2:226882/1‑251 (MQ=255)                                                                                                                         cGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGGCACGGGAATATCTGCGGAAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCTCCAATCACTACCGGACCCTGAATAACGGCTCCATGACAAATACGGGTATTCGCGCCGATAACAATTTCACCTGCAGGTTCATCAATATTTACCTGGTCTGCa                                                                               <  2:256894/250‑1 (MQ=25)                                                                                                                           aGCATGATAACGACTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGACACGGGAATATCTGCGGCAGCGGGCTTGAAATGGGTTAAGAGATAAGAGATGTCTGGAAATATAGGGGCAAATCAAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAAAGGCAACATGACAAATACGGGTATGAGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATc                                                                            >  1:214367/1‑251 (MQ=25)                                                                                                                                  taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAatat                                                                     >  1:126881/1‑251 (MQ=255)                                                                                                                                      gAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATAttttt                                                                 <  1:189727/251‑1 (MQ=255)                                                                                                                                          ggCAGCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGTCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTctct                                                             <  1:217339/251‑1 (MQ=255)                                                                                                                                           gCATCCGCATAGTGCGCTGGGATATCGATCGCCAAGGGAATATCTGAGGCAGCGGGCTTGTACTGGGTTAAGTGATAACAGATGTCTGGAACTATAGGGGCCAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACATGGTCTGCAATCTGAATATTTTTCTCtt                                                            >  2:242699/1‑251 (MQ=255)                                                                                                                                                             gggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGAt                                          <  1:22450/251‑1 (MQ=255)                                                                                                                                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTAGCAGGCAGTTTGCGCCAATCACGACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGGTTTTCAGTTCCTGTGACAGAAGGTTGGGAAt        <  1:155078/251‑1 (MQ=255)                                                                                                                                                                                                     tGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGc    <  1:31413/249‑1 (MQ=255)                                                                                                                                                                                                      gTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCgg  <  1:29471/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                          GCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGG  >  NC_000913/382345‑382791

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.