Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	4,620,828	43 bp→AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT	intergenic (+2/‑13)	deoB → / → deoD	phosphopentomutase/purine nucleoside phosphorylase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4620827	2 (0.070)	15 (0.700) +AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT	15/356	0.4	93.2%	intergenic (+1/‑56)	deoB/deoD	phosphopentomutase/purine nucleoside phosphorylase
	?	NC_000913	4620871 =	1 (0.030)					intergenic (+45/‑12)	deoB/deoD	phosphopentomutase/purine nucleoside phosphorylase

ACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTAC  >  NC_000913/4620621‑4621101                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                       cacggttgtgatccgacctggaccggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgata                                                                                                                                                                                                                                        <  1:372486‑M2/1‑1 (MQ=255)                    ggaccggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTa                                                                                                                                                                                                                      <  1:332559‑M2/19‑1 (MQ=255)                         ggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCcac                                                                                                                                                                                                                <  1:4978‑M2/25‑1 (MQ=255)                                    cacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAAc                                                                                                                                                                                                                              <  1:36646‑M2/11‑1 (MQ=255)                                                    ccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTc                                                                                                                                                                                      <  1:5594‑M2/51‑1 (MQ=255)                                                                tacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGtt                                                                                                                                                                          <  1:159288‑M2/63‑1 (MQ=255)                                                                 acggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGtttt                                                                                                                                                                        <  1:168401‑M2/65‑1 (MQ=255)                                                                    gccagaaagtaaaaccgggctcactgggtcaccgtgaaaccgtcgcggatagtggtcagactctggtataatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcgccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTGCCCCACTCATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGAt                                                                                                                                                                     <  2:365885‑M2/68‑1 (MQ=255)                                                                       cgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGcc                                                                                                                                                                  <  1:273144‑M2/71‑1 (MQ=255)                                                                                   cgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCt                                                                                                                                                      <  1:210379‑M2/83‑1 (MQ=255)                                                                                        tcactgggtccccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTAcc                                                                                                                                                                                                                    <  2:69239‑M2/21‑1 (MQ=255)                                                                                          actgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCtgcgtgcg                                                                                                                                               <  1:109560‑M2/90‑1 (MQ=255)                                                                                                                                    ggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCgtgaagt                                                                                                      <  1:200006‑M2/131‑1 (MQ=255)                                                                                                                                    ggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCgtgaa                                                                                                        <  1:340567‑M2/129‑1 (MQ=255)                                                                                                                                                                           atccatgttctgacaatgggtgcggcttgacgctctctccccggcccacttccacagggagagggtgaaaacatgttgatAAGGACAAAACAATGGCTACCCCACACATCAATGCAGAAATGGGCGATTTCGCTGACGTAGTTGTGATTCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTAc                                                              <  1:357053‑M2/171‑1 (MQ=255)                                                                                                                                                                                                                 ccccggccctctcccgcagggagagtgtgaaatcatgttgatAAGGATAAAACAATTGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTTAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATg                          <  2:265939‑M2/207‑1 (MQ=255)                                                                                                                                                                                                                               gTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCa                                                                         <  1:97971/188‑1 (MQ=255)                                                                                                                                                                                                                                        agggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTac  >  1:263968‑M2/20‑250 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                       ACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTAC  >  NC_000913/4620621‑4621101

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