Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,580,649	IS5 (–) +4 bp	coding (811‑814/1395 nt)	hsdS ←	type I restriction enzyme EcoKI specificity subunit

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2288919 =	NA (NA)	25 (0.560)	22/482	0.7	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	= 4580652	0 (0.000)					coding (811/1395 nt)	hsdS	type I restriction enzyme EcoKI specificity subunit
*	?	NC_000913	= 3653230	NA (NA)	35 (0.780)	33/482	0.3	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	4580649 =	0 (0.000)					coding (814/1395 nt)	hsdS	type I restriction enzyme EcoKI specificity subunit

AATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCT  >  NC_000913/2288676‑2289165                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       aaTTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTg                                                                                                                                                                                                                                                 <  2:378630/251‑1 (MQ=255)         aaCCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGATGGCTATCTTTTTTCAGGACAAAAACATCAACTATTTCTCCCTACTCAACAGTTTATGCCAATGCTTTGGGTGAATATAAGGAAGGGGCGCATa                                                                                                                                                                                                                                           >  1:78645/1‑249 (MQ=255)                    aaTACAGAGTATTCAGCAGCGCCATTCATTTAACTGTCTTGTTACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTAATATGCTGTTTTACGCCATGAGTGATTTTCGTATTATTGGCGGTCCAATAAACCAGCTGGATATATTTTTTCAGGACAACGACATCAACTTGATAACCAGAATCAAAAGTTGATGCCAATGTTATGGGTAAAAATACGGAAGGTTCGAATAATCGGGGAAAtt                                                                                                                                                                                                                               <  2:261563/250‑1 (MQ=255)                                          cattcattTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCCATGCTATGGGGGAATATAAGGAAGGTGCGAATAAGc                                                                                                                                                                                                                                        >  1:385/1‑219 (MQ=255)                                          cattcattTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGc                                                                                                                                                                                                                                        <  2:385/219‑1 (MQ=255)                                                                              cAGTATTATTACTTACACCAAATTCCGTTATCAGCTGTTCGCTGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAATCAGCTGTATATCTTTTTCCAGGACAACACCATCAACTGTCTCTCCATACTGAACATTTGATGCCAATGCTATGGGTTAATATAAGGAAGGTGCGAATAAGCGGGGAAAGTCTTCTCTGCTGACTCAGTCATTTCATTTCTCCATGTGGGAGCCGATTTTTGCTCCCGTa                                                                                                                                                                    <  2:77879/251‑1 (MQ=255)                                                                                                                cTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTTCAGGACAAAACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGATATGGGGGAATATAAGGAAGGTGCGAATAAGAGGGGGAATTCTTCTCGGCTGGCTCAGTCATTTCATTTTTTCATGTTTTAGCCGATTTTTTCTCCCGTAAATGCCTTTGATCAGCCTTTTTAGACCGtttttt                                                                                                                                  >  2:407774/1‑248 (MQ=255)                                                                                                                         gCCAGGATTGATTTTCGTATTATTGGCGGGCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAGCTGTCTCACCATACTCAACAGTTGATGCCAATGGTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGTAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATTTTTGAGCGGATTTTTTTTCCCGTGAATGCCTTGAATCAGCCGATTTAGACCGtgtcggctccattta                                                                                                                         <  2:64388/251‑15 (MQ=25)                                                                                                                                         gctcttccgatcttCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGCGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGt                                                                                                         >  1:272252/14‑251 (MQ=255)                                                                                                                                                             aaCCAGCTGGATATATTTTTCCCGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTccc                                                                                      >  2:193484/1‑250 (MQ=255)                                                                                                                                                                             ttttCCAGGACAACACCATCAACTGTCGCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCGTCATGTTTGAGCCGATTTTTGCTCCCGTAAAGGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGTCATTGAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCAtt                                                                        <  1:213856/248‑1 (MQ=255)                                                                                                                                                                                                cAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTg                                                   <  1:89217/250‑1 (MQ=255)                                                                                                                                                                                                                 tCAACAGTTGATGCAACTGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATc                                 >  1:374662/1‑251 (MQ=255)                                                                                                                                                                                                                                                ataAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCt  >  1:276302/1‑251 (MQ=11)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       AATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCCT  >  NC_000913/2288676‑2289165

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