Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,720,136	IS3 (–) +3 bp	coding (88‑90/213 nt)	cspA →	cold shock protein CspA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2170171 =	NA (NA)	71 (0.910)	67/492	0.1	98.6%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	3720136 =	1 (0.010)					coding (88/213 nt)	cspA	cold shock protein CspA
*	?	NC_000913	= 2171428	NA (NA)	90 (1.150)	83/492	0.0	100%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	= 3720138	0 (0.000)					coding (90/213 nt)	cspA	cold shock protein CspA

AGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCC  >  NC_000913/2169933‑2170344                                                                                                                                                                                                                                               |                                                                                                                                                                              aGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGGGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTAccc                                                                                                                                                                     >  1:15197/1‑249 (MQ=255)                                                                                                                                                                   cACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATGCGCCTGAATTTCGCAGTGAAGCCCTGAAGCGTGCTGAACGCATCGGTGTTACTGCCGCAggcc  <  2:15197/250‑4 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                              AGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCC  >  NC_000913/2169933‑2170344

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