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

ACAGCGAGTCGGGTGCCTCGACCCTGCGCGCGGCACGTGAACGGCCACGGGCACCGATCCTCAACCTGACGCCGAACCTGAACACCGCGCGCCGCCTGAGCGTGGCTTGGGGTGTGCATTCGGTGGTCAACGACCGCCTGCGTCAGGTCGACGAGGTGGTTTCCACCGCCCTGGAGATTGCCCAGGCGCAAGGCATGGCCAGCCGTGGCGACACGCTGCTGATCACCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGTGCTGCTGCTGCGCAACCCCCT  >  AE015451/4893236‑4893824                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                             ctgcccagccgaaccgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCgggg                                                                                                                                                                                                                                                                                                    <  1:365067‑M2/11‑1 (MQ=255)  tgctcagcctcaccgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCggggg                                                                                                                                                                                                                                                                                                   >  3:193335‑M2/290‑301 (MQ=255)         ccccaccgaaagccggaacggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTg                                                                                                                                                                                                                                                                                            <  2:440973‑M2/19‑1 (MQ=255)                  aagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgcgcggccgcGATTGCCGGGGGCGCCTTGCGCCCCttt                                                                                                                                                                                                                                                                                   >  2:559343‑M2/274‑301 (MQ=255)                  aagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCCCGTTGCGCCCCttt                                                                                                                                                                                                                                                                                   >  1:133354‑M2/274‑301 (MQ=255)                          aaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGacac                                                                                                                                                                                                                                                                           >  1:568338‑M2/266‑301 (MQ=255)                          aaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCCCGTTGCGCCCCTTTCGCGacac                                                                                                                                                                                                                                                                           >  1:368407‑M2/266‑301 (MQ=255)                            ggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACaa                                                                                                                                                                                                                                                                         <  3:568353‑M2/38‑1 (MQ=255)                             ccctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAg                                                                                                                                                                                                                                                                        <  3:424528‑M2/39‑1 (MQ=255)                                ggccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGcc                                                                                                                                                                                                                                                                     <  1:193333‑M2/42‑1 (MQ=255)                                agacgtgtgctcttccgatctccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGcc                                                                                                                                                                                                                                                                     <  3:250302‑M2/42‑1 (MQ=255)                                                cctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CAc                                                                                                                                                                                                                                                    >  1:316756‑M2/244‑301 (MQ=255)                                                 ctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACg                                                                                                                                                                                                                                                   <  4:434891‑M2/59‑1 (MQ=255)                                                    gccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTg                                                                                                                                                                                                                                                <  4:716990‑M2/62‑1 (MQ=255)                                                     ccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGcgccccgttcgcgacacaaggccagatcggaagagcgtcgtgtg                                                                                                                                                                                                                                                >  1:250296‑M2/239‑280 (MQ=255)                                                              ccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAagggg                                                                                                                                                                                                                                      >  3:109383‑M2/230‑300 (MQ=255)                                                                  cggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCACAC‑CGGGCACAAgggaaagg                                                                                                                                                                                                                                  >  2:419108‑M2/226‑298 (MQ=255)                                                                   ggcatcgatggctcggaagaaaacggcggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCa                                                                                                                                                                                                                                 <  4:392103‑M2/77‑1 (MQ=255)                                                                    ccttccgatcttcggaagaaaacgccggtgccttcatgcccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAg                                                                                                                                                                                                                                <  1:416231‑M2/78‑1 (MQ=255)                                                                      atcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGcagac                                                                                                                                                                                                                              >  1:563215‑M2/222‑299 (MQ=255)                                                                       cacacgacgctcttccgatctccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTcc                                                                                                                                                                                                                             <  1:118520‑M2/81‑1 (MQ=255)                                                                           gggctcggaagaaaacgccggggccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCtgt                                                                                                                                                                                                                         <  1:38263‑M2/85‑1 (MQ=255)                                                                               tcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACaagggaatgcagagatcggaaga                                                                                                                                                                                                                     >  3:416241‑M2/213‑290 (MQ=255)                                                                                       aaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCCCAAGGGAATGCAGTCCCTGTGGGAGCGGCCtt                                                                                                                                                                                                             >  1:146300‑M2/205‑301 (MQ=255)                                                                                       aaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCtt                                                                                                                                                                                                             <  3:631241‑M2/97‑1 (MQ=255)                                                                                        aacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTtg                                                                                                                                                                                                            <  4:419210‑M2/98‑1 (MQ=255)                                                                                            ccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACgctggcacaagggaatgcagtccagatcggaagagcacacgtct                                                                                                                                                                                                        >  3:118524‑M2/200‑280 (MQ=255)                                                                                                    ttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTgg                                                                                                                                                                                                <  3:316764‑M2/110‑1 (MQ=255)                                                                                                       atgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAgg                                                                                                                                                                                             <  4:559481‑M2/113‑1 (MQ=255)                                                                                                       atgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAgg                                                                                                                                                                                             >  3:264586‑M2/189‑301 (MQ=255)                                                                                                        gtgtgctcttccgatctgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAggg                                                                                                                                                                                            <  4:549417‑M2/114‑1 (MQ=255)                                                                                                               ggccagctaagccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGcc                                                                                                                                                                                     <  3:563230‑M2/121‑1 (MQ=255)                                                                                                                  cagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGGCGCGATTGGAGGGCAAagccccc                                                                                                                                                                                  >  1:300562‑M2/178‑299 (MQ=255)                                                                                                                   agctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAagggaatgcagtcccttggggagcggcctttgttcgcgcttgaggcgcaagcacctcc                                                                                                                                                                                 >  2:286582‑M2/177‑273 (MQ=255)                                                                                                                   agctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCCGGCACAAGGGAATGCATGCCACGTGGGAGCGGCCGTGTGTCCGGCTTtgagggcaaaagccccc                                                                                                                                                                                 >  2:143716‑M2/177‑294 (MQ=255)                                                                                                                    cctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGccctccc                                                                                                                                                                                <  3:300569‑M2/126‑1 (MQ=255)                                                                                                                     ctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGccctccct                                                                                                                                                                               >  3:7407‑M2/175‑301 (MQ=255)                                                                                                                         gccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTCGCCCCAGGGAATGCAGTCCCTGTGGGAGCGGCcttgtgtcgcgagcggagggatggtgggagagcgtgg                                                                                                                                                                           >  2:549279‑M2/171‑284 (MQ=255)                                                                                                                           cagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATAc                                                                                                                                                                         >  1:351037‑M2/169‑301 (MQ=255)                                                                                                                                 ggaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgctcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGctct                                                                                                                                                                   <  1:109381‑M2/139‑1 (MQ=255)                                                                                                                                 cgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGctct                                                                                                                                                                   <  3:146306‑M2/139‑1 (MQ=255)                                                                                                                                  gaagccgtgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTa                                                                                                                                                                  <  1:264581‑M2/140‑1 (MQ=255)                                                                                                                                    agccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTAcc                                                                                                                                                                <  3:133360‑M2/142‑1 (MQ=255)                                                                                                                                     gcctggggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTAccc                                                                                                                                                               <  3:368418‑M2/143‑1 (MQ=255)                                                                                                                                                     ccagcgacgagctggacaacaacagcggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGGGTCGCGATTGGAGGCCAAACCCCCCGCGGAAAAACGCTCGACcctaccgcgggaacccag                                                                                                                                               >  2:10776‑M2/143‑293 (MQ=255)                                                                                                                                                       agcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGGTTGGAGGGCAAGGACctccgtgaataccaccctaccacacgccggaaaaacgaag                                                                                                                                             >  2:22890‑M2/141‑282 (MQ=255)                                                                                                                                                                                       ctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAACGCCTTCCCTGACTACCGCTCTACCCGACTGCCGCAACCCGA‑GCCCCAAGTGGGCCACCCGCTCCCCcatcacccc                                                                                                            >  2:503005‑M2/109‑297 (MQ=255)                                                                                                                                                                                         acttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTcc                                                                                                          >  1:328294‑M2/107‑301 (MQ=255)                                                                                                                                                                                            ccgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTCCact                                                                                                       <  3:328302‑M2/198‑1 (MQ=255)                                                                                                                                                                                            cccgcgccctggatgcgctagtcgtcaccgaaccgacccgcaccaacgtcaagcacttccgagccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATAACGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTCCact                                                                                                       <  4:10776‑M2/198‑1 (MQ=255)                                                                                                                                                                                             cgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTCCactg                                                                                                      <  4:286635‑M2/199‑1 (MQ=255)                                                                                                                                                                                                              tgatcgtcaccgacccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACAcc                                                                                     <  1:7407‑M2/216‑1 (MQ=255)                                                                                                                                                                                                                          aacccacccgccccaacgtcaacgacctccgcgccttcctgatccttgagactgcccaatccggagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTc                                                                         <  4:22897‑M2/228‑1 (MQ=255)                                                                                                                                                                                                                                     accaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGccc                                                              >  3:536644‑M2/63‑301 (MQ=255)                                                                                                                                                                                                                                               acgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATCCCGCTCTACCCGACCGCGGAAAACCGAGGCCCACAGCGGCCCC‑CCGCCCCCCCATCACCCTCCACTGACTGCCCCATGTACCCCAAAACTTTCTTCAAAACGTGccccgacgcggcc                                                    >  2:158337‑M2/53‑295 (MQ=255)                                                                                                                                                                                                                                                            catcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCa                                       >  4:380305‑M2/40‑301 (MQ=255)                                                                                                                                                                                                                                                                     ccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAATAGGCCCA‑CCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACAGCGTTACTCAACGGCTTCa                              >  4:199282‑M2/31‑301 (MQ=255)                                                                                                                                                                                                                                                                             ctgccaaatcatgggcggccgcGAGGGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGtg                      <  2:380229‑M2/279‑1 (MQ=255)                                                                                                                                                                                                                                                                                cccaatcatgagcggccgcGATGGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGtgctg                   <  4:143746‑M2/282‑1 (MQ=255)                                                                                                                                                                                                                                                                                   aatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCGA‑CCGCTCGCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCTATAGTCAACGGCTTCAGCCAGGtgctgctg                >  3:352661‑M2/17‑301 (MQ=255)                                                                                                                                                                                                                                                                                                 gcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGTGCTGCTGCTGCGCAACCCCCt  <  1:536630‑M2/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                             ACAGCGAGTCGGGTGCCTCGACCCTGCGCGCGGCACGTGAACGGCCACGGGCACCGATCCTCAACCTGACGCCGAACCTGAACACCGCGCGCCGCCTGAGCGTGGCTTGGGGTGTGCATTCGGTGGTCAACGACCGCCTGCGTCAGGTCGACGAGGTGGTTTCCACCGCCCTGGAGATTGCCCAGGCGCAAGGCATGGCCAGCCGTGGCGACACGCTGCTGATCACCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGC‑CACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCA‑CCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGTGCTGCTGCTGCGCAACCCCCT  >  AE015451/4893236‑4893824

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