Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	AE015451	1,531,732	Δ107 bp	32.7%	coding (38‑144/630 nt)	PP_1344 ←	conserved protein of unknown function

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 1531731	5 (0.120)	7 (0.180)	7/322	NT	32.7%	coding (145/630 nt)	PP_1344	conserved protein of unknown function
	?	AE015451	1531839 =	24 (0.610)					coding (37/630 nt)	PP_1344	conserved protein of unknown function

CACATGGCAGTGCTCGCGGGCAGCCGGTTGCAACTGGCTTTCCAGCAGGCGCTGCAAGTGCTCCAGGCGCTCGGCCTGGTTCAGCAGCAGGCGCAGCGGGCGCACCTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAGCCAAAGACAAGACTGGTTCTCGATAGCTTTCATGTTGCCAAGCCTTGAACTCAGTGAAAATGCCCCTATCTAAAACAGGCCCCTGCCAAAACGCTGTATCAGCATCGTGGAAATCCACCAC  >  AE015451/1531546‑1531963                                                                                                                                                                                          |                                                                                                                                                                                                                                         cacaTGGCGGTGCTCGCGGGCAGCCGGTTGCAACTGGCTTTCCAGCAGGCGCTGCAAGTGCTCCAGGCGCTCGGCCTGGTTCAGCAGCAGGCGCAGCGGGCGCACCTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagcctt                                                                                                                                                                                        <  1:161293‑M1/236‑51 (MQ=255)                                                                                    gcagcagGCGCAGCGGGCGCACCTGGCTCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatg                                                                                                                                                                        <  2:168842‑M1/169‑67 (MQ=255)                                                                                                          cTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgtt                                                                                                                                                                                                  <  2:484455‑M1/121‑41 (MQ=255)                                                                                                                              ggggCTGGGCGGGGGGGGGTTTGTAGGCCGTGGGGGTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAGCCAAAGACAAGACTGGTTCTCGATAGCTTTCATGTTGCCAAGCCTTGAACTCAGTGAAAATg                                                             <  1:657195/232‑1 (MQ=255)                                                                                                                                                         ccATGGGGGTACGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatgcccctatctaaaacaggcccctgccaaaacgctgtatcagcatcgtggaaatccaccactttcctcaccatcgtttccgggtagaatgctcgttcgcatgcggcctatgggctgcacgggcgactcatggggccgc                                <  1:104342‑M1/236‑203 (MQ=255)                                                                                                                                                             gggggTACGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatgcccctatctaaaacaggcccctgccaaaacgctgtatcagcatcgtggaaatccaccactttcctcaccatcgtttccgggtagaatgctcgttcgcatgaggcct                                                              >  2:327185‑M1/1‑30 (MQ=255)                                                                                                                                                              ggggTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGc                                                                                                                                                           <  1:427803/108‑1 (MQ=255)                                                                                                                                                              ggggTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGc                                                                                                                                                           >  2:427802/1‑108 (MQ=255)                                                                                                                                                               gggTACGCGTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatgcccctatctaaaacaggcccctgccaaaacgctgtatcagcatcgtggaaatccaccactttcctcaccatcgtttccgggtagaatgctcgttcgcatgcggcctatgggctgcacgggcgactcatggggccgccctcc                           >  1:755986‑M1/1‑28 (MQ=255)                                                                                                                                                                      cgTGAAGTGACGTGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatgcccctatctaaaacaggcccctgccaaaacgctgtatcagcatcgtggaaatccaccactttcctcaccatcgtttccgggtagaatgctcgttcgcatgcggcctatgggctgcacgggcgactcatggggccgccctccatccctaa                   >  1:159303‑M1/1‑21 (MQ=255)                                                                                                                                                                             tGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTCCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAGCCAAAGACAAGACTGGTTCTCGATAGCTTTCATGtt                                                                                       >  1:39893/1‑161 (MQ=255)                                                                                                                                                                             tGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTCCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAGCCAAAGACAAGACTGGTTCTCGATAGCTTTCATGtt                                                                                       <  2:39893/161‑1 (MQ=255)                                                                                                                                                                                  tGGGCACAGgcagccaaagacaagactggttctcgatagctttcatgttgccaagccttgaactcagtgaaaatgcccctatctaaaacaggcccctgccaaaacgctgtatcagcatcgtggaaatccaccactttcctccccatcgtttccgggtagaatgctcgttcgcatgcggcctatgggctgcacgggcgactcatggggccgccctccctccctaagtgtgga            >  1:38897‑M1/1‑9 (MQ=255)                                                                                                                                                                                        aggTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAGCCAAAGACAAGACTGGTTCTCGATAGCTTTCATGTTGCCAAGCCTTGAACTCAGTGAAAATGCCCCTATCTAAAACAGGCCCCTGCCAAAACGCTGTATCAGCATCGTGGAAATCcaccac  <  1:442036/233‑1 (MQ=255)                                                                                                                                                                                          |                                                                                                                                                                                                                                         CACATGGCAGTGCTCGCGGGCAGCCGGTTGCAACTGGCTTTCCAGCAGGCGCTGCAAGTGCTCCAGGCGCTCGGCCTGGTTCAGCAGCAGGCGCAGCGGGCGCACCTGGCGCAGCAGGGCGGAAGGTGGCTGGGCGGGGGAGGGTTTGTAGGCCATGGGGGTACGCGTGAAGTGACGTGGGCACAGTCTAACATCTCCGTTGCCTGTGCCGGCCTCTTCGCGGGGGAACCCGCTCCTGCCAGGGTTATCGGTATCTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAGCCAAAGACAAGACTGGTTCTCGATAGCTTTCATGTTGCCAAGCCTTGAACTCAGTGAAAATGCCCCTATCTAAAACAGGCCCCTGCCAAAACGCTGTATCAGCATCGTGGAAATCCACCAC  >  AE015451/1531546‑1531963

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