Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	AE015451	4,892,114	1412 bp→GCGGCCGC	coding (5‑1416/1416 nt)	pyk →	Pyruvate kinase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	AE015451	4892114	4893525	1412	48 [0]	[0] 41	pyk	Pyruvate kinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 4892113	0 (0.000)	41 (0.580) +GCGGCCGC	37/584	0.7	100%	coding (1275/1275 nt) coding (4/1416 nt)	ttuD pyk	putative hydroxypyruvate reductase Pyruvate kinase
	?	AE015451	4893526 =	0 (0.000)					intergenic (+1/‑207)	pyk/PP_4302	Pyruvate kinase/Urea transporter

CGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGACGCCTGATAAAAAAGTTAAGATCCTTGCCACCCTTGGCCCTGCGATCAACGGCATCGACGATATCCGCCAGCTGGTCGAAGCCGGGGTGAACATCTTCCGCCTCAACTTCAGCCACGGCGAACATGCCGACCACGCCCTGCGCTACCAGTGGATCCGCGAAGTCGAGCAACAGCTGAACTACCCGCTGGGCATCCTCATGGACCTGCAAGGGCCGAAGCTGCGCGTTGGCCGCTTCGCCGAAGGCAAGGTGCAGTTGCAACGCGGGCAGGCCCTGCGCCTGGATCT  >  AE015451/4891813‑4892399                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                               cgggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGa                                                                                                                                                                                                                                                                                                >  4:463034/1‑301 (MQ=255) cgggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCAACAACGTCAACGACTTCCGCCCCCTCCTGATCCTTGATCCTGCCCAATCATGa                                                                                                                                                                                                                                                                                                >  2:243489/1‑301 (MQ=255)  agttcagacgtgtgctcttccgatctcctcaccgaaagcctgaaaggcctgccgGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAg                                                                                                                                                                                                                                                                                               <  4:136225‑M1/275‑2 (MQ=255)     cGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcgg                                                                                                                                                                                                                                                                                            >  4:478966‑M1/1‑297 (MQ=255)     cGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcgg                                                                                                                                                                                                                                                                                            <  4:352438‑M1/301‑5 (MQ=255)        aaCGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccg                                                                                                                                                                                                                                                                                         >  3:396580‑M1/1‑294 (MQ=255)             cGAGTTCCTGCTCAGCCTAACCGAAAGCCTGAAAGGCCGGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgatt                                                                                                                                                                                                                                                                                    <  4:610365‑M1/301‑13 (MQ=255)                    cTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggg                                                                                                                                                                                                                                                                             >  1:215869‑M1/1‑282 (MQ=255)                            ccTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCCCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCCCCATCCTGATCCTTGAGACTGCCCAACCATGAgagatcggaagagcgtcggtgagggca                                                                                                                                                                                                                                                                     >  2:136190‑M1/1‑274 (MQ=255)                             cTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCGGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgc                                                                                                                                                                                                                                                                    <  1:296075‑M1/301‑29 (MQ=255)                                     aaGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgccccttt                                                                                                                                                                                                                                                            <  1:396566‑M1/301‑37 (MQ=255)                                      aGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttc                                                                                                                                                                                                                                                           <  1:458921‑M1/301‑38 (MQ=255)                                                gcgtgccGGGGGTGTAGGCCCGGGCAGGTGAAACCGCCGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaag                                                                                                                                                                                                                                                 <  2:618710‑M1/298‑48 (MQ=255)                                                 ccTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaagg                                                                                                                                                                                                                                                <  4:243554‑M1/301‑49 (MQ=255)                                                     ccGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCCCCATCCTGATCCTTGAACCTGCCCAATCATGAccggccgccattgccgtgggcgcgttgcgcccctttcgcgccacacggcccc                                                                                                                                                                                                                                            >  2:446835‑M1/1‑249 (MQ=255)                                                       gggCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGCCTGCCCAATCATGAgcggccgcgattgccgggggcgcgtttcgcgcctttcgcgacccaaggccgccc                                                                                                                                                                                                                                          >  1:194290‑M1/1‑247 (MQ=255)                                                                   ccgcgccggTGACACCGAGCGCTTCGATGGCTCCGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCATCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccac                                                                                                                                                                                                                              <  2:462888‑M1/297‑67 (MQ=255)                                                                   ccTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccac                                                                                                                                                                                                                              <  3:215875‑M1/301‑67 (MQ=255)                                                                          ggTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcc                                                                                                                                                                                                                       >  4:489435‑M1/1‑228 (MQ=255)                                                                                    aCGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGCCTGCCCAATCATGAacggccgggattgccgggggcgagttgcgcccctttccgcacacaagcccgctccccagtggccccgcgggcaacggggagcg                                                                                                                                                                                                             >  2:589225‑M1/1‑218 (MQ=255)                                                                                    aCGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatg                                                                                                                                                                                                             >  3:365375‑M1/1‑218 (MQ=255)                                                                                        cATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagt                                                                                                                                                                                                         <  3:444734‑M1/301‑88 (MQ=255)                                                                                              tGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCCTGAgcggccgcgatagccgggggcgcgctgcgcccctttcgcgacccaaggccgcccccaccgggccacgctggcacaagggaatggcgtcccttt                                                                                                                                                                                                   >  2:64924‑M1/1‑208 (MQ=255)                                                                                                   gggaAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggag                                                                                                                                                                                              <  1:135805‑M1/300‑99 (MQ=255)                                                                                                       ggagacagccggggcCGTCCTGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggc                                                                                                                                                                                          <  2:448775‑M1/294‑103 (MQ=255)                                                                                                            aCGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgatagccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgt                                                                                                                                                                                     <  4:64939‑M1/301‑108 (MQ=255)                                                                                                               ccGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcaggccctgtgggagcggccttgcgtc                                                                                                                                                                                  >  2:183466‑M1/1‑191 (MQ=255)                                                                                                               ccGGGGCCTCCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtc                                                                                                                                                                                  <  2:489284‑M1/301‑111 (MQ=255)                                                                                                                           tGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgcccccacagggccacgctgggccaagggaatgcagtccctgttgcagcggccttgtgtcccgattgggggg                                                                                                                                                                      >  2:603941‑M1/1‑179 (MQ=255)                                                                                                                                                         cctggggccggcgggccagcgaccagctGGACAACAACAACGGCTACGGCTGCTTCGCGCCGCTGGATGGGCTTGTCGTCACCGAACCGACCCGCACCACCGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccg                                                                                                                                        <  2:478815‑M1/288‑153 (MQ=255)                                                                                                                                                                       cccttaccctacccgccgccctgccgatcttacggctacttagccgcgctggatacgctgatcgtcaccGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccg                                                                                                                          <  1:200278‑M1/271‑167 (MQ=255)                                                                                                                                                                         cAGCGACGCGCTTGACCACAACCCCTGCGACGGCTACTTCGCCGCGCTGGATGCGCTGAGCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgag                                                                                                                        <  1:365365‑M1/301‑169 (MQ=255)                                                                                                                                                                                         aacaacaacGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccac                                                                                                        <  4:589427‑M1/301‑185 (MQ=255)                                                                                                                                                                                                  ggctacggctacTTCGCCGCGCTGGATGCGCTGATCTTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctccccc                                                                                               <  4:183516‑M1/301‑194 (MQ=255)                                                                                                                                                                                                     tacggctacTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgagatcggaagagcacacgtctgaactccag                                                                                            >  3:200284‑M1/1‑105 (MQ=255)                                                                                                                                                                                                          ctacTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctaaccgaccgcggaaaaccgaggcccaaagaggcccacccctcccccccccaccc                                                                                       >  2:373151‑M1/1‑100 (MQ=255)                                                                                                                                                                                                             cTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAAGCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacaggggcacgctggcacaagggaatgaagtccctgtgggagcggccttgtgtcgcgatgggagggcaaagcccgccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccataaaccgcca                                                                                    >  4:690063‑M1/1‑97 (MQ=255)                                                                                                                                                                                                                    gcgcTGGAGGCGCTGAGCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactg                                                                             <  3:194296‑M1/301‑212 (MQ=255)                                                                                                                                                                                                                      gcTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctccccgaccgcggaaaaccgaggcccacagcgtcccaccgctcccccatccacctcccactcccgcc                                                                           >  2:527835‑M1/1‑88 (MQ=255)                                                                                                                                                                                                                          cctcttccgatcttcacCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccat                                                                       <  3:300846‑M1/293‑218 (MQ=255)                                                                                                                                                                                                                           aTGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctgggacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatg                                                                      <  4:604143‑M1/301‑219 (MQ=255)                                                                                                                                                                                                                             gcgcTGATCGTCACCGAACCGACCCGCACCAACGTCAGCGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggccacccgctcccccatcaacctccactgactgccccatgta                                                                    >  1:303688‑M1/1‑81 (MQ=255)                                                                                                                                                                                                                                       tCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaacgaggcccaccgctcccccatcaacctcccctgactgccccatagatcggcagcgg                                                          >  1:300839‑M1/1‑71 (MQ=255)                                                                                                                                                                                                                                          ccgaaccgaCCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttc                                                       <  4:373254‑M1/301‑234 (MQ=255)                                                                                                                                                                                                                                                     cGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtg                                            <  4:446973‑M1/301‑245 (MQ=255)                                                                                                                                                                                                                                                                    aacTTGCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgagattgccgggggcgggttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccac                             <  4:528002‑M1/300‑260 (MQ=255)                                                                                                                                                                                                                                                                                           tGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggcaaaccgaggcccaaagaggcccaccgctcccccatcaccctccactcactgccccctgtacaccaaacctttcgtaaccccgctccccgtctcggcccagccgttactcaacggcttcgccc      >  2:17393‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                                                                                               aCTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactcaacggcttcagccaggt  >  1:360902‑M1/1‑15 (MQ=255)                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                               CGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGACGCCTGATAAAAAAGTTAAGATCCTTGCCACCCTTGGCCCTGCGATCAACGGCATCGACGATATCCGCCAGCTGGTCGAAGCCGGGGTGAACATCTTCCGCCTCAACTTCAGCCACGGCGAACATGCCGACCACGCCCTGCGCTACCAGTGGATCCGCGAAGTCGAGCAACAGCTGAACTACCCGCTGGGCATCCTCATGGACCTGCAAGGGCCGAAGCTGCGCGTTGGCCGCTTCGCCGAAGGCAAGGTGCAGTTGCAACGCGGGCAGGCCCTGCGCCTGGATCT  >  AE015451/4891813‑4892399

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