Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	AE015451	4,808,177	Δ103 bp	64.0%	intergenic (‑286/‑167)	PP_4234 ← / → dsbD‑II	Aldehyde oxidase and xanthine dehydrogenase/thiol/disulfide interchange protein 2

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 4808176	0 (0.000)	8 (0.190)	8/326	NT	64.0%	intergenic (‑285/‑270)	PP_4234/dsbD‑II	Aldehyde oxidase and xanthine dehydrogenase/thiol/disulfide interchange protein 2
	?	AE015451	4808280 =	9 (0.220)					intergenic (‑389/‑166)	PP_4234/dsbD‑II	Aldehyde oxidase and xanthine dehydrogenase/thiol/disulfide interchange protein 2

GCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGCAAGCCAACCCGTTCGATGTGA  >  AE015451/4808053‑4808515                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            gaaaagctgaccttggcgtcaactttacagcaattgggtgggtagaggcaaagtcaagctgacaatgtgtcgcggtgtgtatagggcaatgggcggaattcgcgctttggcgggtaggtgtttggcgggaggttttcagcgcctgggaggtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccCCTCAAAAt                                                                                                                                                                                                                                     <  1:144228‑M2/9‑1 (MQ=255)                                                                ttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcCCGCGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCtt                                                                                                                                                                      <  1:145669‑M2/163‑1 (MQ=255)                                                                           gtgtgtatagggcaatgggcggaattcgcgctttggcgggtaggtgtttggcgggaggttttcagcgcctgggaggtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTg                                                                                                                                                            <  1:553680‑M2/82‑1 (MQ=255)                                                                              ttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAAt                                                                                                                                                                                                            <  1:252732‑M2/125‑1 (MQ=255)                                                                                             gcggaattcgcgctttggcgggtaggtgtttggcgggaggttttcagcgcctgggaggtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccCCTCAAAATTCCAAcgc                                                                                                                                                                                                                             <  2:561832‑M2/17‑1 (MQ=255)                                                                                             gcggaattcgcgctttggcgggtaggtgtttggcgggaggttttcagcgcctgggaggtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccCCTCAAAATTCCAAcgc                                                                                                                                                                                                                             >  1:561834‑M2/136‑152 (MQ=255)                                                                                                   ttcgcgctttggcgggtaggtgtttggcgggaggttttcagcgcctgggaggtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttccgcgcctgggagatcgagcgccgcccCCTCAAAATTc                                                                                                                                                                                                                                   <  2:570070‑M2/11‑1 (MQ=255)                                                                                                         ctttggcgggtaggtgtttggcgggaggttttcagcgcctgggaggtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGa                                                                                                                                                                                <  2:528969‑M2/62‑1 (MQ=255)                                                                                                             ggcgggtaggtggttggcgggaggtttgcagcgcctgggaggtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAAt                                                                                                                                                                                                            <  2:105747‑M2/34‑1 (MQ=255)                                                                                                                             gcgAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTa                                                                                                                                     >  2:673026/1‑208 (MQ=255)                                                                                                                                        cGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGacc                                                                                                                                                         >  1:122386‑M2/2‑177 (MQ=255)                                                                                                                                           gggcggCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCaaa                                                                                                                                                                                      <  1:740209/145‑1 (MQ=255)                                                                                                                                           gggcggCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCaaa                                                                                                                                                                                      >  2:740207/1‑145 (MQ=255)                                                                                                                                                        tctcACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAAcccc                                                                                                                 <  1:340417/201‑1 (MQ=255)                                                                                                                                                         ctcACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAAt                                                                                                                                                                                                            <  1:685592/109‑1 (MQ=255)                                                                                                                                                         ctcACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAAt                                                                                                                                                                                                            >  2:685590/1‑109 (MQ=255)                                                                                                                                                           cgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGa                                                                                                                                                           <  2:733503‑M2/83‑1 (MQ=255)                                                                                                                                                                                     cccccTCAAAATTCCAACGCCGCCCACCCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTa                                                                                                                                     <  2:137360/149‑1 (MQ=255)                                                                                                                                                                                               aTTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCtt                                                    <  2:741530/223‑1 (MQ=255)                                                                                                                                                                                                ttCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATgc                                                                    <  1:243694/206‑1 (MQ=255)                                                                                                                                                                                                                         cAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCGACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATgc                                                                    <  1:673028/181‑1 (MQ=255)                                                                                                                                                                                                                             gccgcccCCTCAAAATTCCAAcgc                                                                                                                                                                                                                             >  2:505097‑M2/8‑24 (MQ=255)                                                                                                                                                                                                                                    ccTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGCCAGCCAACCCGTTCGATGTGa  >  1:85045/1‑236 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            GCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGCAAGCCAACCCGTTCGATGTGA  >  AE015451/4808053‑4808515

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