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	63 [0]	[0] 61	pyk	Pyruvate kinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 4892113	0 (0.000)	60 (0.560) +GCGGCCGC	54/584	0.9	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

CGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGACGCCTGATAAAAAAGTTAAGATCCTTGCCACCCTTGGCCCTGCGATCAACGGCATCGACGATATCCGCCAGCTGGTCGAAGCCGGGGTGAACATCTTCCGCCTCAACTTCAGCCACGGCGAACATGCCGACCACGCCCTGCGCTACCAGTGGATCCGCGAAGTCGAGCAACAGCTGAACTACCCGCTGGGCATCCTCATGGACCTGCAAGGGCCGAAGCTGCGCGTTGGCCGCTTCGCCGAAGGCAAGGTGCAGTTGCAACGCGGGCAGGCCCTGCGCCTGGATCTGGACAAG  >  AE015451/4891813‑4892406                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                      ggggcgTAACGACGAGTTCGTGCTCAGCCTCACCAGAGGGCTGAAAGGCCTGGCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACTAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGa                                                                                                                                                                                                                                                                                                       <  2:539585/300‑1 (MQ=255)     cGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcgg                                                                                                                                                                                                                                                                                                   <  3:143836‑M1/301‑5 (MQ=255)       aaacGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggcc                                                                                                                                                                                                                                                                                                 <  1:539472‑M1/300‑7 (MQ=255)                    ctgcccagccgaaccgaaagcctgaAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggg                                                                                                                                                                                                                                                                                    <  1:365067‑M1/289‑20 (MQ=255)                     tGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccggggg                                                                                                                                                                                                                                                                                   >  3:193335‑M1/1‑281 (MQ=255)                            ccccaccGAAAGCCGGAACGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttg                                                                                                                                                                                                                                                                            <  2:440973‑M1/298‑28 (MQ=255)                                     aaGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGCgcggccgcgattgccgggggcgccttgcgccccttt                                                                                                                                                                                                                                                                   >  2:559343‑M1/1‑265 (MQ=255)                                     aaGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcccgttgcgccccttt                                                                                                                                                                                                                                                                   >  1:133354‑M1/1‑265 (MQ=255)                                             aaGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacac                                                                                                                                                                                                                                                           >  1:568338‑M1/1‑257 (MQ=255)                                             aaGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcccgttgcgcccctttcgcgacac                                                                                                                                                                                                                                                           >  1:368407‑M1/1‑257 (MQ=255)                                               ggCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaa                                                                                                                                                                                                                                                         <  3:568353‑M1/301‑47 (MQ=255)                                                ccctGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaag                                                                                                                                                                                                                                                        <  3:424528‑M1/300‑48 (MQ=255)                                                   ggcCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggcc                                                                                                                                                                                                                                                     <  1:193333‑M1/300‑51 (MQ=255)                                                   agacgtgtgctcttccgatctccggtgacaccgacggcatcgatGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggcc                                                                                                                                                                                                                                                     <  3:250302‑M1/280‑51 (MQ=255)                                                                   ccTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccac                                                                                                                                                                                                                                     >  1:316756‑M1/1‑235 (MQ=255)                                                                    cTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacg                                                                                                                                                                                                                                    <  4:434891‑M1/301‑68 (MQ=255)                                                                       gCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctg                                                                                                                                                                                                                                 <  4:716990‑M1/301‑71 (MQ=255)                                                                        ccGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgccccgttcgcgacacaaggccagatcggaagagcgtcgtgtg                                                                                                                                                                                                                                >  1:250296‑M1/1‑230 (MQ=255)                                                                                 ccGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggg                                                                                                                                                                                                                       >  3:109383‑M1/1‑221 (MQ=255)                                                                                     cGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggcacaccgggcacaagggaaagg                                                                                                                                                                                                                   >  2:419108‑M1/1‑217 (MQ=255)                                                                                      ggCATCGATGGCTCGGAAGAAAACGGCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgca                                                                                                                                                                                                                  <  4:392103‑M1/301‑86 (MQ=255)                                                                                       ccttccgatcttcggaagaaaacGCCGGTGCCTTCATGCCCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcag                                                                                                                                                                                                                 <  1:416231‑M1/290‑87 (MQ=255)                                                                                         aTCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagac                                                                                                                                                                                                               >  1:563215‑M1/1‑213 (MQ=255)                                                                                          cacacgacgctcttccgatctccggtgccttcatgaccccggccAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtcc                                                                                                                                                                                                              <  1:118520‑M1/280‑90 (MQ=255)                                                                                              gggcTCGGAAGAAAACGCCGGGGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgt                                                                                                                                                                                                          <  1:38263‑M1/300‑94 (MQ=255)                                                                                                  tCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagagatcggaaga                                                                                                                                                                                                      >  3:416241‑M1/1‑204 (MQ=255)                                                                                                          aaaCGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcccaagggaatgcagtccctgtgggagcggcctt                                                                                                                                                                                              >  1:146300‑M1/1‑196 (MQ=255)                                                                                                          aaaCGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggcctt                                                                                                                                                                                              <  3:631241‑M1/301‑106 (MQ=255)                                                                                                           aaCGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttg                                                                                                                                                                                             <  4:419210‑M1/301‑107 (MQ=255)                                                                                                               ccGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccagatcggaagagcacacgtct                                                                                                                                                                                         >  3:118524‑M1/1‑191 (MQ=255)                                                                                                                       ttCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattgg                                                                                                                                                                                 <  3:316764‑M1/301‑119 (MQ=255)                                                                                                                          aTGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagg                                                                                                                                                                              <  4:559481‑M1/301‑122 (MQ=255)                                                                                                                          aTGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagg                                                                                                                                                                              >  3:264586‑M1/1‑180 (MQ=255)                                                                                                                           gtgtgctcttccgatctgccagcgccgaagccttgGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggaggg                                                                                                                                                                             <  4:549417‑M1/284‑123 (MQ=255)                                                                                                                                  ggCCAGCTAAGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagcc                                                                                                                                                                      <  3:563230‑M1/301‑130 (MQ=255)                                                                                                                                     cAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtggcgcgattggagggcaaagccccc                                                                                                                                                                   >  1:300562‑M1/1‑169 (MQ=255)                                                                                                                                      aGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtcccttggggagcggcctttgttcgcgcttgaggcgcaagcacctcc                                                                                                                                                                  >  2:286582‑M1/1‑168 (MQ=255)                                                                                                                                      aGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgccggcacaagggaatgcatgccacgtgggagcggccgtgtgtccggctttgagggcaaaagccccc                                                                                                                                                                  >  2:143716‑M1/1‑168 (MQ=255)                                                                                                                                       cctACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccc                                                                                                                                                                 <  3:300569‑M1/300‑135 (MQ=255)                                                                                                                                        cTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccct                                                                                                                                                                >  3:7407‑M1/1‑166 (MQ=255)                                                                                                                                            gCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctcgccccagggaatgcagtccctgtgggagcggccttgtgtcgcgagcggagggatggtgggagagcgtgg                                                                                                                                                            >  2:549279‑M1/1‑162 (MQ=255)                                                                                                                                              cAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaatac                                                                                                                                                          >  1:351037‑M1/1‑160 (MQ=255)                                                                                                                                                    ggaAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGCTCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctct                                                                                                                                                    <  1:109381‑M1/300‑148 (MQ=255)                                                                                                                                                    cGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctct                                                                                                                                                    <  3:146306‑M1/301‑148 (MQ=255)                                                                                                                                                     gAAGCCGTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctcta                                                                                                                                                   <  1:264581‑M1/301‑149 (MQ=255)                                                                                                                                                       aGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacc                                                                                                                                                 <  3:133360‑M1/301‑151 (MQ=255)                                                                                                                                                        gCCTGGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctaccc                                                                                                                                                <  3:368418‑M1/301‑152 (MQ=255)                                                                                                                                                                        ccAGCGACGAGCTGGACAACAACAGCGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgggtcgcgattggaggccaaaccccccgcggaaaaacgctcgaccctaccgcgggaacccag                                                                                                                                >  2:10776‑M1/1‑134 (MQ=255)                                                                                                                                                                          aGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcggttggagggcaaggacctccgtgaataccaccctaccacacgccggaaaaacgaag                                                                                                                              >  2:22890‑M1/1‑132 (MQ=255)                                                                                                                                                                                                          ctacTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaacgccttccctgactaccgctctacccgactgccgcaacccgagccccaagtgggccacccgctcccccatcacccc                                                                                              >  2:503005‑M1/1‑100 (MQ=255)                                                                                                                                                                                                            acTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctcc                                                                                            >  1:328294‑M1/1‑98 (MQ=255)                                                                                                                                                                                                               ccgccgcgCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccact                                                                                         <  3:328302‑M1/300‑207 (MQ=255)                                                                                                                                                                                                               cccgcgccctggatgcgCTAGTCGTCACCGAACCGACCCGCACCAACGTCAAGCACTTCCGAGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataacgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccact                                                                                         <  4:10776‑M1/293‑207 (MQ=255)                                                                                                                                                                                                                cgccgcGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactg                                                                                        <  4:286635‑M1/301‑208 (MQ=255)                                                                                                                                                                                                                                 tGATCGTCACCGACCCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacacc                                                                       <  1:7407‑M1/301‑225 (MQ=255)                                                                                                                                                                                                                                             aacccaCCCGCCCCAACGTCAACGACCTCCGCGCCTTCCTGATCCTTGAGACTGCCCAATCCGGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtc                                                           <  4:22897‑M1/301‑237 (MQ=255)                                                                                                                                                                                                                                                        aCCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccc                                                >  3:536644‑M1/1‑54 (MQ=255)                                                                                                                                                                                                                                                                  aCGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaatcccgctctacccgaccgcggaaaaccgaggcccacagcggccccccgcccccccatcaccctccactgactgccccatgtaccccaaaactttcttcaaaacgtgccccgacgcggcc                                      >  2:158337‑M1/1‑44 (MQ=255)                                                                                                                                                                                                                                                                               cATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactca                         >  4:380305‑M1/1‑31 (MQ=255)                                                                                                                                                                                                                                                                                        ccTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaataggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacagcgttactcaacggcttca                >  4:199282‑M1/1‑22 (MQ=255)                                                                                                                                                                                                                                                                                                cTGCCAAATCATGGgcggccgcgagggccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactcaacggcttcagccaggtg        <  2:380229‑M1/301‑288 (MQ=255)                                                                                                                                                                                                                                                                                                   cccAATCATGAgcggccgcgatggccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactcaacggcttcagccaggtgctg     <  4:143746‑M1/301‑291 (MQ=255)                                                                                                                                                                                                                                                                                                      aaTCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggccgaccgctcgcccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggctatagtcaacggcttcagccaggtgctgctg  >  3:352661‑M1/1‑8 (MQ=255)                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                      CGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGACGCCTGATAAAAAAGTTAAGATCCTTGCCACCCTTGGCCCTGCGATCAACGGCATCGACGATATCCGCCAGCTGGTCGAAGCCGGGGTGAACATCTTCCGCCTCAACTTCAGCCACGGCGAACATGCCGACCACGCCCTGCGCTACCAGTGGATCCGCGAAGTCGAGCAACAGCTGAACTACCCGCTGGGCATCCTCATGGACCTGCAAGGGCCGAAGCTGCGCGTTGGCCGCTTCGCCGAAGGCAAGGTGCAGTTGCAACGCGGGCAGGCCCTGCGCCTGGATCTGGACAAG  >  AE015451/4891813‑4892406

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑