Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	136,615	Δ5 bp	3.5%	coding (299‑303/348 nt)	yacC ←	PulS_OutS family protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 136614	111 (0.870)	4 (0.030)	4/506	NT	3.5%	coding (304/348 nt)	yacC	PulS_OutS family protein
	?	NC_000913	136620 =	111 (0.880)					coding (298/348 nt)	yacC	PulS_OutS family protein

TAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/136493‑136614 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGGCTTTCATGTCGAAGGTGTCGTAATTACTGAGGTCCCACTGGTTTTGCTGAGCGAAAAAGACCAGTGCGCGACGAATTTGCCCGTTAGGTAAGTTCTGGTAACCACAATCGTTCTTCAGAAAGACAAAAACTGCCGTTAAATCGGCCATATCTTCGGCTTCAGACTCGCTGAGCGCGTAACTGTTTGCGCAGACGGCCATCAGGCTGCCGAATAACA  >  NC_000913/136620‑136895                                                                                                                                                                                                                                                                                                                                                                                                                TAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGAAGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCATTTTTTAGCGACGGGAATGCCAATGCCGATGAGATAGCGGTAGCTGTCTTCACCGAGGGCTTTCATGTCGAAGGTGTCGTAATTACTGAGGTCCCACTGGTTTTGCTGAGCGAGAAAGAACAGTGCGAGACGAATTTGCCCGTTAGGTAAGTTCTGGTTACCACAA                                                                                                               >  1:490296/1‑289                                                           TACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGGCTTTCATGTCGAAGGTGTCGTAATTACTGAGGTCCCACTGGTTTTGCCGAGCGAAAAAGACCAGTGCGAGACGAATTGGCCCGTTAGGGAAGTTCTGATAACCACAATCGTTCTTGAGCAAGCCAAAAACTGCAGTAAAATCGGCAATATATTCGGCTTCAGGCT                                                     >  2:327743/1‑289                                                                                       GTAGCCAAGCAGGCTGAAGCAATCGCGGGCGATGGCATTGTTTAGCACCGGGAATGCCAATGCCGCTGATATCGCGGTATCTGTCTTGACAGCGGGCTTTCAGGTCCAAGTTGTCGGAATTCCTCGGGTCCCACTGGTTTTGCGGAGCGAAAAATACCAGCGCGCGACTAATCTGCCCGTTAGGTAAGTTCTGGTAACCACAATCGTTCTTCAGAAATACAAAAACTGCCGTTAAATCGGCCATATCTTCGGCTTCAGACTCGCTGAGCGCGTAACTGTGTGCGCAGA                          <  2:490296/288‑1                                                                                                              CGCGGGCCAGAGAATTTTTTAGCTACGGGACGGCCAATGCCGCTGAGATCGCGGTATCTGTCTTCACCGAGGGCTTTCATGTCGAAGTGGTCGTAATTACTGAGGTCCCACTGGTTTTGCTGAGCGAAAAAGACCAGTGCGCGACGAATTTGCCCGTTAGGTAAGTTCTGGTAACCACAATCGTTCTTCAGAAAGACAAAAACTGCCGTTAAATCGGCCATATCTTCGGCTTCAGACTCGCTGAGCGCGTAACTGTTTGCGCAGACGGCCATCAGGCTGCCGAATAACA  <  1:327743/289‑1                                                                                                                                                                                                                                                                                                                                                                                                                TAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/136493‑136614 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGGCTTTCATGTCGAAGGTGTCGTAATTACTGAGGTCCCACTGGTTTTGCTGAGCGAAAAAGACCAGTGCGCGACGAATTTGCCCGTTAGGTAAGTTCTGGTAACCACAATCGTTCTTCAGAAAGACAAAAACTGCCGTTAAATCGGCCATATCTTCGGCTTCAGACTCGCTGAGCGCGTAACTGTTTGCGCAGACGGCCATCAGGCTGCCGAATAACA  >  NC_000913/136620‑136895

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