Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	AE015451	1,044,759	Δ21 bp	7.5%	coding (467‑487/702 nt)	PP_0904 →	InaA protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 1044758	42 (1.050)	3 (0.080)	3/324	NT	7.5%	coding (466/702 nt)	PP_0904	InaA protein
	?	AE015451	1044780 =	33 (0.860)					coding (488/702 nt)	PP_0904	InaA protein

ACAGTTTCGAGCAGTTGGGCGTGCGTGTACCGCGCATCGTCTTCTGTGGCGCGCAGCGCGATGCCGAGCACCAGTGGCGCGCGCTGCTGGTCAGCGAAGCGCTGGACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/1044522‑1044758 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcacgTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGCTACGCCGCCATTCGTCGCTCAATGACGCGGAATGGCGGACGCTGCTCTACTTTTACAAGATGGCGTTTGGCAGCGCTGTCAAA  >  AE015451/1044780‑1044982                                                                                                                                                                                                                                                                                                                                                                                                                                                          ACAGTTTCGAGCAGTTGGGCGTGCGTGTACCGCGCATCGTCTTCTGTGGCGCGCAGCGCGATGCCGAGCACCAGTGGCGCGCGCTGCTGGTCAGCGAAGCGCTGGACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCAC                                                                                                                                                                                                              <  1:333266/236‑1 ACAGTTTCGAGCAGTTGGGCGTGCGTGTACCGCGCATCGTCTTCTGTGGCGCGCAGCGCGATGCCGAGCACCAGTGGCGCGCGCTGCTGGTCAGCGAAGCGCTGGACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCAC                                                                                                                                                                                                              <  1:637487/236‑1 ACAGTTTCGAGCAGTTGGGCGTGCGTGTACCGCGCATCGTCTTCTGTGGCGCGCAGCGCGATGCCGAGCACCAGTGGCGCGCGCTGCTGGTCAGCGAAGCGCTGGACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCAC                                                                                                                                                                                                              <  1:663485/236‑1                                                                                                                                                                     GTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGCTACGCCGCCATTCGTCGCTCAATGACGCGGAATGGCGGACGC                                           >  1:688191/1‑235                                                                                                                                                                           CATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGCTACGCCGCCATTCGTCGCTCAAT                                                              >  2:426312/1‑210                                                                                                                                                                                                          ACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGCTACGCCGCCATTCGTCGCTCAATGACGCGGAATGGCGGACGCTGCTCTACTTTTACAAGATGGCGTTTGGCAGCGCTGT      <  1:426313/235‑1                                                                                                                                                                                                                                          CACGTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGCTACGCCGCCATTCGTCGCTCAATGACGCGGAATGGCGGACGCTGCTCTACTTTTACAAGATGGCGTTTGGCAGCGCTGTCAAA  >  1:599561/1‑207                                                                                                                                                                                                                                                                                                                                                                                                                                                          ACAGTTTCGAGCAGTTGGGCGTGCGTGTACCGCGCATCGTCTTCTGTGGCGCGCAGCGCGATGCCGAGCACCAGTGGCGCGCGCTGCTGGTCAGCGAAGCGCTGGACGGTTTTGTCGAACTCGACACCTGGCATGCCGAAGGTGCCCGCGAGCGTTACCCGCAAGTGGTGCATGAGCGCATGCTCAAGGACCTGGCGGACAACCTGGCGCGGATGCACCTGGGGCACTGGCAGCACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/1044522‑1044758 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcacgTATTCATCAAGGTAATTGGCGAAGGCGAGCAGGCCCGCGTCGAAGTGGCGCTGCTGGACCTGGAAAAGTGCAGGCGGCGTATCAGCTGCCAGCGGGCGGCCGGCAATGACTTGCGTCAGCTACGCCGCCATTCGTCGCTCAATGACGCGGAATGGCGGACGCTGCTCTACTTTTACAAGATGGCGTTTGGCAGCGCTGTCAAA  >  AE015451/1044780‑1044982

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