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	55 [0]	[0] 53	pyk	Pyruvate kinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 4892113	0 (0.000)	51 (0.610) +GCGGCCGC	45/536	1.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

CAACTACAGCGAGTCGGGTGCCTCGACCCTGCGCGCGGCACGTGAACGGCCACGGGCACCGATCCTCAACCTGACGCCGAACCTGAACACCGCGCGCCGCCTGAGCGTGGCTTGGGGTGTGCATTCGGTGGTCAACGACCGCCTGCGTCAGGTCGACGAGGTGGTTTCCACCGCCCTGGAGATTGCCCAGGCGCAAGGCATGGCCAGCCGTGGCGACACGCTGCTGATCACCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGTGCTGCTGCTGCGCAACCCCCT  >  AE015451/4893231‑4893824                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                           agttcctgctcagcctcaccgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTg                                                                                                                                                                                                                                                                                                        >  1:736961‑M2/296‑300 (MQ=255)     cctgctcagcctcaccgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCggg                                                                                                                                                                                                                                                                                                   <  1:961667‑M2/10‑1 (MQ=255)                 caccgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCgggg                                                                                                                                                                                                                                                                                                  <  2:556066‑M2/11‑1 (MQ=255)                 caccgaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCgggg                                                                                                                                                                                                                                                                                                  >  1:556087‑M2/280‑290 (MQ=255)                     gaaagcctgaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGcccc                                                                                                                                                                                                                                                                                    >  1:940916‑M2/276‑300 (MQ=255)                             gaaaggcctgccgggcgtgtacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGac                                                                                                                                                                                                                                                                           >  2:563085‑M2/268‑301 (MQ=255)                             gaaaggccggccgggcgtgtacgccctggccggggaaaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTcg                                                                                                                                                                                                                                                                               <  2:105954‑M2/30‑1 (MQ=255)                                                tacgccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGACGCGTTGCGCCCCTTTc                                                                                                                                                                                                                                                                                >  2:424041‑M2/249‑277 (MQ=255)                                                    ccctggccggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGcc                                                                                                                                                                                                                                                     <  2:736935‑M2/56‑1 (MQ=255)                                                            ggtgacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcattagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACCCAAGGCCGCTCCCACAGGGcc                                                                                                                                                                                                                                                     >  1:820962‑M2/237‑292 (MQ=255)                                                               gacaccgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGcac                                                                                                                                                                                                                                           >  2:675225‑M2/234‑299 (MQ=255)                                                                    cgacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAggg                                                                                                                                                                                                                                      >  1:391265‑M2/229‑299 (MQ=255)                                                                     gacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTTCGCCCCTTTCGCCACACAAGGcc                                                                                                                                                                                                                                                                   >  1:476797‑M2/228‑269 (MQ=255)                                                                     gacggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGaa                                                                                                                                                                                                                                    <  1:563106‑M2/73‑1 (MQ=255)                                                                       cggcatcgatggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGa                                                                                                                                                                                                                                     <  2:940890‑M2/72‑1 (MQ=255)                                                                           atcgatggctcggaagaaaacgccggtgccgtcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCa                                                                                                                                                                                                                                <  2:622473‑M2/77‑1 (MQ=255)                                                                                 ggctcggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCtgt                                                                                                                                                                                                                        <  1:414696‑M2/85‑1 (MQ=255)                                                                                     cggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCtg                                                                                                                                                                                                                         <  2:967799‑M2/84‑1 (MQ=255)                                                                                      ggaagaaaacgccggtgccttcatgaccccggccagctacgcctgcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGTGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGa                                                                                                                                                                                                                                     <  2:561639‑M2/72‑1 (MQ=255)                                                                                      ggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCCGGc                                                                                                                                                                                                                                             >  2:20602‑M2/211‑274 (MQ=255)                                                                                      ggaagaaaacgccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcacaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGATGGCACAagggc                                                                                                                                                                                                                                     >  1:561660‑M2/211‑281 (MQ=255)                                                                                      ggaagaaaacgccggtgccttcatgaccccggcaagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGa                                                                                                                                                                                                                    <  2:820936‑M2/89‑1 (MQ=255)                                                                                        aagaaaacgccggtcccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCgg                                                                                                                                                                                                                <  2:264638‑M2/93‑1 (MQ=255)                                                                                                 ccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTcg                                                                                                                                                                                                                                                                               >  1:329444‑M2/200‑229 (MQ=255)                                                                                                 ccggtgccttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTcg                                                                                                                                                                                                                                                                               <  2:329433‑M2/30‑1 (MQ=255)                                                                                                      gccttcatgaccccggccagctacgccagcgccgaagccttgggcccgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTcgc                                                                                                                                                                                                     >  1:737903‑M2/195‑298 (MQ=255)                                                                                                         ttcatgaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGc                                                                                                                                                                                                               >  1:182129‑M2/192‑285 (MQ=255)                                                                                                              gaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcttcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTcc                                                                                                                                                                                                                            <  2:373712‑M2/81‑1 (MQ=255)                                                                                                              gaccccggccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGc                                                                                                                                                                                                               <  2:182125‑M2/94‑1 (MQ=255)                                                                                                                      ccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCtgt                                                                                                                                                                                                                        <  1:162397‑M2/85‑1 (MQ=255)                                                                                                                      ccagctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCtgt                                                                                                                                                                                                                        >  2:162393‑M2/179‑263 (MQ=255)                                                                                                                         gctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACCCAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTgg                                                                                                                                                                                                                      >  1:418738‑M2/176‑262 (MQ=255)                                                                                                                         gctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTgg                                                                                                                                                                                                                      <  2:418718‑M2/87‑1 (MQ=255)                                                                                                                          ctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTtgt                                                                                                                                                                                                          >  1:66582‑M2/175‑273 (MQ=255)                                                                                                                          ctacgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGccc                                                                                                                                                                                   >  1:969520‑M2/175‑296 (MQ=255)                                                                                                                             cgccagcgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGgcg                                                                                                                                                                                                                                                                                               >  2:922492‑M2/172‑185 (MQ=255)                                                                                                                              gccagcgccgaagccttgggcctgtcagccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTccc                                                                                                                                                                                                                           <  2:795334‑M2/82‑1 (MQ=255)                                                                                                                                   cgccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTtgt                                                                                                                                                                                                          <  2:66580‑M2/99‑1 (MQ=255)                                                                                                                                    gccgaagccttgggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTtg                                                                                                                                                                                                           <  2:391253‑M2/98‑1 (MQ=255)                                                                                                                                               gggcctgtcggccagcgacgagctggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTtgt                                                                                                                                                                                                          <  2:476776‑M2/99‑1 (MQ=255)                                                                                                                                                                       ggacaacaacaacggctacggctacttcgccgcgctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCa                                                                                                                                        <  2:737877‑M2/165‑1 (MQ=255)                                                                                                                                                                                   cggctacggctacttcgccgctctggatgcgctgatcgtcaccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTcccc                                                                                                                     <  1:675250‑M2/184‑1 (MQ=255)                                                                                                                                                                                                                            ccgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCAAAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGAACGCGGAAAACCGAGGCCCAAAGAGGCCCACCg                                                                                                                           >  1:351748‑M2/77‑254 (MQ=255)                                                                                                                                                                                                                             cgaaccgacccgcaccaacgtcaacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCCGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc                                                                            <  1:20604‑M2/225‑1 (MQ=255)                                                                                                                                                                                                                                                   aacgacttccgcgccatcctgatccttgagactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTggg                                                      >  2:581656‑M2/54‑300 (MQ=255)                                                                                                                                                                                                                                                          tccgcgccggcctgatccttgcgactgcccaataatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGaa                                                                                                                                                                           <  1:922518‑M2/130‑1 (MQ=255)                                                                                                                                                                                                                                                                                gactgcccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCg                                                                                                                           <  2:351737‑M2/178‑1 (MQ=255)                                                                                                                                                                                                                                                                                     cccaatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGTGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGc                                                                                                                                           <  1:122968‑M2/162‑1 (MQ=255)                                                                                                                                                                                                                                                                                       caatcatgagcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGCCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGCCCTTGTGTCGCGGTTGGAGGGCAAAtccctccc                                                                                                                                                                               >  2:677512‑M2/18‑143 (MQ=255)                                                                                                                                                                                                                                                                                       caatcatcagcggccgcGATTGCCGGGGGCGCTTTGCGCGCCTTTCGCGACACAAGGCCGCTCCCACAGGGCAACGCTGGCACAAGTGAATGCAGTCCCTGTGGTAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGtgctg                   <  2:969494‑M2/282‑1 (MQ=255)                                                                                                                                                                                                                                                                                               agcggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAACGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc                                                                            >  2:787266‑M2/10‑234 (MQ=255)                                                                                                                                                                                                                                                                                                  ggccgcGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc                                                                            >  2:351544‑M2/7‑231 (MQ=255)                                                                                                                                                                                                                                                                                                        gaTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAACAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGTTCCGCTGCTGCGCAACCCCCt  >  2:406438/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                           CAACTACAGCGAGTCGGGTGCCTCGACCCTGCGCGCGGCACGTGAACGGCCACGGGCACCGATCCTCAACCTGACGCCGAACCTGAACACCGCGCGCCGCCTGAGCGTGGCTTGGGGTGTGCATTCGGTGGTCAACGACCGCCTGCGTCAGGTCGACGAGGTGGTTTCCACCGCCCTGGAGATTGCCCAGGCGCAAGGCATGGCCAGCCGTGGCGACACGCTGCTGATCACCGCTGGCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACGGCTTCAGCCAGGTGCTGCTGCTGCGCAACCCCCT  >  AE015451/4893231‑4893824

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