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

GTCAACTACAGCGAGTCGGGTGCCTCGACCCTGCGCGCGGCACGTGAACGGCCACGGGCACCGATCCTCAACCTGACGCCGAACCTGAACACCGCGCGCCGCCTGAGCGTGGCTTGGGGTGTGCATTCGGTGGTCAACGACCGCCTGCGTCAGGTCGACGAGGTGGTTTCCACCGCCCTGGAGATTGCCCAGGCGCAAGGCATGGCCAGCCGTGGCGACACGCTGCTGATCACCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGT  >  AE015451/4893229‑4893803                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                       cgagttcctgctcagcctaaccgaaagcctgaaaggccggccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGAtt                                                                                                                                                                                                                                                                                     <  4:610365‑M2/4‑1 (MQ=255)        ctgctcagcctcaccgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCgggg                                                                                                                                                                                                                                                                              >  1:215869‑M2/291‑301 (MQ=255)                 ctcaccgaaagcctgaaaggcctgccgggcgtgtacgcccgggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTgc                                                                                                                                                                                                                                                                     <  1:296075‑M2/20‑1 (MQ=255)                         aagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCttt                                                                                                                                                                                                                                                             <  1:396566‑M2/28‑1 (MQ=255)                          agcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTc                                                                                                                                                                                                                                                            <  1:458921‑M2/29‑1 (MQ=255)                                    gcgtgccgggggtgtaggcccgggcaggtgaaaccgccggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAg                                                                                                                                                                                                                                                 <  2:618710‑M2/39‑1 (MQ=255)                                     cctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAgg                                                                                                                                                                                                                                                <  4:243554‑M2/40‑1 (MQ=255)                                         ccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgccccatcctgatccttgaacctgcccaatcatgaccggccgcCATTGCCGTGGGCGCGTTGCGCCCCTTTCGCG‑CCACACggcccc                                                                                                                                                                                                                                            >  2:446835‑M2/258‑299 (MQ=255)                                           gggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagcctgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTTCGCGCCTTTCGCG‑ACCCAAGGCcgccc                                                                                                                                                                                                                                          >  1:194290‑M2/256‑299 (MQ=255)                                                       cctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCAc                                                                                                                                                                                                                              <  3:215875‑M2/58‑1 (MQ=255)                                                       ccgcgccggtgacaccgagcgcttcgatggctccgaagaaaacgccggtgccttcatgaccccggccatctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCAc                                                                                                                                                                                                                              <  2:462888‑M2/58‑1 (MQ=255)                                                              ggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGgcc                                                                                                                                                                                                                       >  4:489435‑M2/237‑300 (MQ=255)                                                                        acggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagcctgcccaatcatgaacggccggGATTGCCGGGGGCGAGTTgcgcccctttc‑cgcacacaagcccgctccccagtggccccgcgggcaacggggagcg                                                                                                                                                                                                             >  2:589225‑M2/227‑274 (MQ=255)                                                                        acggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATg                                                                                                                                                                                                             >  3:365375‑M2/227‑301 (MQ=255)                                                                            catcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGt                                                                                                                                                                                                         <  3:444734‑M2/79‑1 (MQ=255)                                                                                  tggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcctgagcggccgcGATAGCCGGGGGCGCGCTGCGCCCCTTTCGCG‑ACCCAAGGCCGCCCCCACCGGGCCACGCTGGCACAAGGGAATGGCGTCccttt                                                                                                                                                                                                   >  2:64924‑M2/217‑299 (MQ=255)                                                                                       gggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAg                                                                                                                                                                                              <  1:135805‑M2/90‑1 (MQ=255)                                                                                           ggagacagccggggccgtcctgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGc                                                                                                                                                                                          <  2:448775‑M2/94‑1 (MQ=255)                                                                                                acgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATAGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTtgt                                                                                                                                                                                     <  4:64939‑M2/99‑1 (MQ=255)                                                                                                   ccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGGCCCTGTGGGAGCGGCCTTGCGTc                                                                                                                                                                                  >  2:183466‑M2/200‑301 (MQ=255)                                                                                                   ccggggcctccatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTc                                                                                                                                                                                  <  2:489284‑M2/102‑1 (MQ=255)                                                                                                               tgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCCCCCACAGGGCCACGCTGGGCCAAGGGAATGCAGTCCCTGTTGCAGCGGCCTTGTGTCCCGATTgggggg                                                                                                                                                                      >  2:603941‑M2/188‑301 (MQ=255)                                                                                                                                             cctggggccggcgggccagcgaccagctggacaacaacaacggctacggctgcttcgcgccgctggatgggcttgtcgtcaccgaaccgacccgcaccaccgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCg                                                                                                                                        <  2:478815‑M2/144‑1 (MQ=255)                                                                                                                                                           cccttaccctacccgccgccctgccgatcttacggctacttagccgcgctggatacgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCg                                                                                                                          <  1:200278‑M2/158‑1 (MQ=255)                                                                                                                                                             cagcgacgcgcttgaccacaacccctgcgacggctacttcgccgcgctggatgcgctgagcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAg                                                                                                                        <  1:365365‑M2/160‑1 (MQ=255)                                                                                                                                                                             aacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCAc                                                                                                        <  4:589427‑M2/176‑1 (MQ=255)                                                                                                                                                                                      ggctacggctacttcgccgcgctggatgcgctgatcttcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTccccc                                                                                               <  4:183516‑M2/185‑1 (MQ=255)                                                                                                                                                                                         tacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATaccgctctacccgaccgcggaaaaccgagatcggaagagcacacgtctgaactccag                                                                                            >  3:200284‑M2/114‑273 (MQ=255)                                                                                                                                                                                              ctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTAACCGACCGCGGAAAACCGAGGCCCAAAGAGGCccacccctcccccccccaccc                                                                                       >  2:373151‑M2/109‑293 (MQ=255)                                                                                                                                                                                                 cttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaagcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGGCACGCTGGCACAAGGGAATGAAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATGGGAGGGCAAAGCCCGCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATAAACCGCCa                                                                                    >  4:690063‑M2/106‑301 (MQ=255)                                                                                                                                                                                                        gcgctggaggcgctgagcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCactgactg                                                                             <  3:194296‑M2/203‑1 (MQ=255)                                                                                                                                                                                                          gctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTCCCCGACCGCGGAAAACCGAGGCCCACAGCGTCCCACCGCTCCCccatccacctcccactcccgcc                                                                           >  2:527835‑M2/97‑291 (MQ=255)                                                                                                                                                                                                              cctcttccgatcttcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCAt                                                                       <  3:300846‑M2/209‑1 (MQ=255)                                                                                                                                                                                                               atgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGGACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATg                                                                      <  4:604143‑M2/210‑1 (MQ=255)                                                                                                                                                                                                                 gcgctgatcgtcaccgaaccgacccgcaccaacgtcagcgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCACCCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                                    >  1:303688‑M2/90‑301 (MQ=255)                                                                                                                                                                                                                           tcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAACGAGGCCCACCGCTCCCCCATCAACCTCCcctgactgccccatagatcggcagcgg                                                          >  1:300839‑M2/80‑288 (MQ=255)                                                                                                                                                                                                                              ccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc                                                       <  4:373254‑M2/225‑1 (MQ=255)                                                                                                                                                                                                                                         cgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTg                                            <  4:446973‑M2/236‑1 (MQ=255)                                                                                                                                                                                                                                                        aacttgcgcgccatcctgatccttgagactgcccaatcatgagcggccgaGATTGCCGGGGGCGGGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCAc                             <  4:528002‑M2/251‑1 (MQ=255)                                                                                                                                                                                                                                                                               tgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGCAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCACCCTCCACTCACTGCCCCCTGTACACCAAACCTTTCGTAACCCCGCTCCCCGTCTCGGCCCAGCCGTTACTCAACGgcttcgccc      >  2:17393‑M2/28‑297 (MQ=255)                                                                                                                                                                                                                                                                                   actgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGt  >  1:360902‑M2/24‑301 (MQ=255)                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                       GTCAACTACAGCGAGTCGGGTGCCTCGACCCTGCGCGCGGCACGTGAACGGCCACGGGCACCGATCCTCAACCTGACGCCGAACCTGAACACCGCGCGCCGCCTGAGCGTGGCTTGGGGTGTGCATTCGGTGGTCAACGACCGCCTGCGTCAGGTCGACGAGGTGGTTTCCACCGCCCTGGAGATTGCCCAGGCGCAAGGCATGGCCAGCCGTGGCGACACGCTGCTGATCACCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCG‑ACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGT  >  AE015451/4893229‑4893803

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