Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	3,484,549	(CCA)4→3	6.3%	coding (62‑64/870 nt)	prkB →	putative phosphoribulokinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3484542 =	108 (1.010)	7 (0.070)	7/410	NT	6.3%	coding (55/870 nt)	prkB	putative phosphoribulokinase
	?	NC_000913	= 3484538	104 (1.020)					coding (51/870 nt)	prkB	putative phosphoribulokinase

GGTATGCAGATATTTGCGAGATTTCCCTTTGCCGCTCTGACCGTATTCAATGAAGGTTTGTTCCAGCAGGCCGAAGTCGTTAGCCTCGGGGCCGAAGTAGCTGATATGCCGCCCGGCGTCGCGCGCTTTGCGGATCGCCATGTCCATTTCCGGGCGGGTGTAACGGTGAAAACTGTCGCCTTCCACCTCAGCTGCATGCAGATTTAACTGCGCGAATATTTTACGAAACGCGAGGCTGGTGGTGGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3484787‑3484542 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggtggtggtCCCCGCGCCGCTGGATCCTGTTACCGCAATGACCGGATGTTTGGCAGACATAGTAACTCCATGACTTAGCTATCATTATTATCTCGTCCGCCGACGGTTGCAGGTTATTCGCGAAACTGACTGCGCGGCATAATATTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTAGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGT  <  NC_000913/3484538‑3484305                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  GGTATGCAGATATTTGCGAGATTTCCCTTTGCCGCTCTGACCGTATTCAATGAAGGTTTGTTCCAGCAGGCCGAAGTCGTTAGCCTCGGGGCCGAAGTAGCTGATATGCCGCCCGGCGTCGCGCGCTTTGCGGATCGCCATGTCCATTTCCGGGCGGGTGTAACGGTGAAAACTGTCGCCTTCCACCTCAGCTGCATGCAGATTTAACTGCGCGAATATTTTACGAAACGCGAGGCTGG                                                                                                                                                                                                                                                   >  2:185477/1‑239     TGCAGATATTTGCGAGATTTCCCTTTGCCGCTCTGACCGTATTCAATGAAGGTTTGTTCCAGCAGGCCGAAGTCGTTAGCCTCGGGGCCGAAGTAGCTGATATGCCGCCCGGCGTCGCGCGCTTTGCGGATCGCCATGTCCATTTCCGGGCGGGTGTAACGGTGAAAACTGTCGCCTTCCACCTCAGCTGCATGCAGATTTAACTGCGCGAATATTTTACGAAACGCGAGGCTGGTGGT                                                                                                                                                                                                                                               >  2:792997/1‑239                                                                               GTTAGCCTCGGGGCCGAAGTAGCTGATATGCCGCCCGGCGTCGCGCGCTTTGCGGATCGCCATGTCCATTTCCGGGCGGGTGTAACGGTGAAAACTGTCGCCTTCCACCTCAGCTGCATGCAGATTTAACTGCGCGAATATTTTACGAAACGCGAGGCTGGTGGTGGTCCCCGCGCCGCTGGATCCTGTTACCGCAATGACCGGATGTTTGGCAGACATAGTAACTCCATGACTTAGCT                                                                                                                                                                     <  2:1069186/239‑1                                                                                                    GCTGATATGCCGCCCGGCGTCGCGCGCTTTGCGGATCGCCATGTCCATTTCCGGGCGGGTGTAACGGTGAAAACTGTCGCCTTCCACCTCAGCTGCATGCAGATTTAACTGCGCGAATATTTTACGAAACGCGAGGCTGGTGGTGGTCCCCGCGCCGCTGGATCCTGTTACCGCAATGACCGGATGTTTGGCAGACATAGTAACTCCATGACTTAGCTAT                                                                                                                                                                   <  1:543050/220‑1                                                                                                    GCTGATATGCCGCCCGGCGTCGCGCGCTTTGCGGATCGCCATGTCCATTTCCGGGCGGGTGTAACGGTGAAAACTGTCGCCTTCCACCTCAGCTGCATGCAGATTTAACTGCGCGAATATTTTACGAAACGCGAGGCTGGTGGTGGTCCCCGCGCCGCTGGATCCTGTTACCGCAATGACCGGATGTTTGGCAGACATAGTAACTCCATGACTTAGCTAT                                                                                                                                                                   >  2:543050/1‑220                                                                                                                                                                                                     TGCAGATTTAACTGCGCGAATATTTTACGAAACGCGAGGCTGGTGGTGGTCCCCGCGCCGCTGGATCCTGTTACCGCAATGACCGGATGTTTGGCA                                                                                                                                                                                              >  1:421868/1‑96                                                                                                                                                                                                     TGCAGATTTAACTGCGCGAATATTTTACGAAACGCGAGGCTGGTGGTGGTCCCCGCGCCGCTGGATCCTGTTACCGCAATGACCGGATGTTTGGCA                                                                                                                                                                                              <  2:421868/96‑1                                                                                                                                                                                                                     CGAATATTTTACGAAACGCGAGGCTGGTGGTGGTCCCCGCGCCGCTGGATCCTGTTACCGCAATGACCGGATGTTTGGCAGACATAGT                                                                                                                                                                                      <  1:212801/88‑1                                                                                                                                                                                                                     CGAATATTTTACGAAACGCGAGGCTGGTGGTGGTCCCCGCGCCGCTGGATCCTGTTACCGCAATGACCGGATGTTTGGCAGACATAGT                                                                                                                                                                                      >  2:212801/1‑88                                                                                                                                                                                                                                             TGGTGGTGGTCCCCGCGCCGCTGGATCCTGTTACCGCAATGACCGGATGTTTGGCAGACATAGTAACTCCATGACTTAGCTATCATTATTATCTCGTCCGCCGACGGTTGCAGGTTATTCGCGAAACTGACTGCGCGGCATAATATTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTAGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGC        <  2:621209/238‑1                                                                                                                                                                                                                                                   TGGTCCCCGCGCCGCTGGATCCTGTTACCGCAATGACCGGATGTTTGGCAGACATAGTAACTCCATGACTTAGCTATCATTATTATCTCGTCCGCCGACGGTTGCAGGTTATTCGCGAAACTGACTGCGCGGCATAATATTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTAGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGT  >  2:974852/1‑238                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  GGTATGCAGATATTTGCGAGATTTCCCTTTGCCGCTCTGACCGTATTCAATGAAGGTTTGTTCCAGCAGGCCGAAGTCGTTAGCCTCGGGGCCGAAGTAGCTGATATGCCGCCCGGCGTCGCGCGCTTTGCGGATCGCCATGTCCATTTCCGGGCGGGTGTAACGGTGAAAACTGTCGCCTTCCACCTCAGCTGCATGCAGATTTAACTGCGCGAATATTTTACGAAACGCGAGGCTGGTGGTGGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3484787‑3484542 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggtggtggtCCCCGCGCCGCTGGATCCTGTTACCGCAATGACCGGATGTTTGGCAGACATAGTAACTCCATGACTTAGCTATCATTATTATCTCGTCCGCCGACGGTTGCAGGTTATTCGCGAAACTGACTGCGCGGCATAATATTGACCGTTTCGTGCAGCTCCGACCATACCAGCACCGCTTCTCCGCACTGTAGCTGGCGTTTGACGTCGGCGACTTTCTGTTCAAGTGTACGCTCATGT  <  NC_000913/3484538‑3484305

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