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	1 (0.020)	22 (0.500) +AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT	17/350	1.2	96.8%	intergenic (+1/‑56)	deoB/deoD	phosphopentomutase/purine nucleoside phosphorylase
	?	NC_000913	4620871 =	1 (0.020)					intergenic (+45/‑12)	deoB/deoD	phosphopentomutase/purine nucleoside phosphorylase

CACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACC  >  NC_000913/4620634‑4621119                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         cgacatggaccggtaatgaccacacgagtgaaaacattccggtaattgtttacggcccgaaagtaaaaccgggcaaactgggtaaacgtgaaaccttcgcggatattggtcagactcgggcaaaatatttaggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcacaccggcactctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAAt                                                                                                                                                                                                                                             <  2:513447‑M2/14‑1 (MQ=255)         accggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccccagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTAccc                                                                                                                                                                                                                                     >  1:378312‑M2/230‑251 (MQ=255)                                    attccggtactggtttacggcccgaaagtaaaaccggtctcactgggtctccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatattgaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAGACAATGGCTACCCCACACATTAATGCAGAAATGGGCGAtt                                                                                                                                                                                                          <  2:573123‑M2/49‑1 (MQ=255)                                    attccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGAtt                                                                                                                                                                                                          <  1:666678‑M2/49‑1 (MQ=255)                                                    acggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggacagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggggcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGttt                                                                                                                                                                                           <  2:408635‑M2/64‑1 (MQ=255)                                                        cccgaaagtaaaaccgggctcactgggtcaccgtgaaactgtcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATg                                                                                                                                                                                      <  1:330344‑M2/69‑1 (MQ=255)                                                                  aaaccgggctcactggggcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATACAAAAATGAGCGAtttcgctgaagtagtttttaagaccgacgaac                                                                                                                                                                            >  2:441669‑M2/173‑237 (MQ=255)                                                                                    caccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctccccccagccctctcccacagggggggggtgaaaacatgttgatAAGGATAAAACAATTGATACCCCACACATTAATGCAGAAATGGGCGATTTCCCTGCCGTCGTTTTTATGCCAggcggcccgctgcggggggaatat                                                                                                                                                           >  1:39396‑M2/155‑240 (MQ=255)                                                                                    caccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGtat                                                                                                                                                           <  1:505988‑M2/96‑1 (MQ=255)                                                                                    caccgtgaaaccttcgcggatattggtcagactctggcaaaataatttggtacttctgatattgaatatggcaaagccatgttctgataataggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGtat                                                                                                                                                           <  2:546951‑M2/96‑1 (MQ=255)                                                                                      ccgtgaaaccttcgcggatattggtcagactctggcagaatattttggtacttctgatatgtaatatggcaaagccatgttctgataatgggttcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTTCCCCACACATTAATGCAGAAATGGGAGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATAtt                                                                                                                                                        <  2:511092‑M2/99‑1 (MQ=255)                                                                                               cttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagggggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTa                                                                                                                                                                                                                            >  1:174499‑M2/144‑174 (MQ=255)                                                                                               cttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTa                                                                                                                                                                                                                            <  2:174499‑M2/31‑1 (MQ=255)                                                                                                  cgtggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACttt                                                                                                                                             <  2:291151‑M2/110‑1 (MQ=255)                                                                                                  cgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACttt                                                                                                                                             <  1:276264‑M2/110‑1 (MQ=255)                                                                                                           tggtcagactctggtaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGa                                                                                                                                    <  1:644309‑M2/119‑1 (MQ=255)                                                                                                             ggcagactctgtcaaaatattttggtgcttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaacccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACCATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGc                                                                                                                                 <  2:378312‑M2/122‑1 (MQ=255)                                                                                                                                     gtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCgg                                                                                                         <  1:601797‑M2/146‑1 (MQ=255)                                                                                                                                                  ggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCtt                                                                                             <  1:619381‑M2/158‑1 (MQ=255)                                                                                                                                                       atggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCg                                                                                        <  1:535314‑M2/163‑1 (MQ=255)                                                                                                                                                          tGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAAAAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGTGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGGTGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTAc                                                                                    >  1:659433/1‑251 (MQ=255)                                                                                                                                                                 catgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACaaa                                                                             <  1:575972‑M2/174‑1 (MQ=255)                                                                                                                                                                 catgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACaaa                                                                             <  1:356465‑M2/174‑1 (MQ=255)                                                                                                                                                                                    cggcgtgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTaa                                                          <  2:144726‑M2/193‑1 (MQ=255)                                                                                                                                                                                                      ccggcactctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGTCCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTAcaaaggccgcaaaatttccgttgctggcggcagtttgttcg                                        <  1:29564‑M2/211‑21 (MQ=255)                                                                                                                                                                                                                            gggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATAGGTATCCCGTCCTGCTCCATCTacac                  >  2:650805‑M2/19‑251 (MQ=255)                                                                                                                                                                                                                                  aaccatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGCCCCGCTGCGTGCGACGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGGGAACAACGTTCGAGGTATGCTGGGCTTCGCCGGTACTTACAAAGACCGCAAAATTTACGTACTGGGTCACGGTATGGGTATCCCGTCCTGCTCCATATACACCaa               >  2:49427‑M2/13‑248 (MQ=255)                                                                                                                                                                                                                                   aacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGGATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAACATTTCCGTAATGGGTCACGGACTGGGTATCCCGTCCTGCTCCATCTACACCAAAGa            >  1:285943‑M2/12‑250 (MQ=255)                                                                                                                                                                                                                                      atgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTg        >  1:35702‑M2/9‑251 (MQ=255)                                                                                                                                                                                                                                        gttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGtata                                                                                                                                                          <  1:536872‑M2/97‑1 (MQ=255)                                                                                                                                                                                                                                        gttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGtata                                                                                                                                                          >  2:536872‑M2/7‑103 (MQ=255)                                                                                                                                                                                                                                            atAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGACGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATACCGTCCTGCTCCATCTACACCAAAGAACTGATCAcc  >  2:269599‑M2/3‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         CACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACC  >  NC_000913/4620634‑4621119

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