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	40 [0]	[0] 26	pyk	Pyruvate kinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 4892113	0 (0.000)	26 (0.660) +GCGGCCGC	26/534	0.8	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

GTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGACGCCTGATAAAAAAGTTAAGATCCTTGCCACCCTTGGCCCTGCGATCAACGGCATCGACGATATCCGCCAGCTGGTCGAAGCCGGGGTGAACATCTTCCGCCTCAACTTCAGCCACGGCGAACATGCCGACCACGCCCTGCGCTACCAGTGGATCCGCGAAGTCGAGCAACAGCTGAACTACCCGCTGGGCATCCTCATGGACCTGCAAGGGCCGAAGCTGCGCGTTGGCCGCTTCGCCGAAGGCAAGGTGCAGTT  >  AE015451/4891818‑4892369                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                 gTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAACCTGCCCAATCATGAgcgg                                                                                                                                                                                                                                                              >  2:334325‑M1/1‑296 (MQ=255)           gTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggcc                                                                                                                                                                                                                                                            <  2:303768‑M1/292‑7 (MQ=255)                 gCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCACCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgc                                                                                                                                                                                                                                                    <  2:179208‑M1/294‑15 (MQ=255)                       ccTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttg                                                                                                                                                                                                                                       <  1:157632‑M1/301‑28 (MQ=255)                            ccGAAAGCCTGAAAGGCCTGCCGGGCGTTTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgtt                                                                                                                                                                                                                                        <  2:305690‑M1/295‑27 (MQ=255)                               aaaGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGGCTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccct                                                                                                                                                                                                                                >  2:105728‑M1/1‑266 (MQ=255)                                aaGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCCCCAACGTCAACGACTTCCGCCCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggc                                                                                                                                                                                                                                             >  2:110810‑M1/1‑265 (MQ=255)                                      gAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccc                                                                                                                                                                                                                                 <  1:110814‑M1/292‑34 (MQ=255)                                      gAAAGGCCTGCCGGGCGGGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccggg                                                                                                                                                                                                                                                <  1:387311‑M1/277‑19 (MQ=255)                                                ccGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattggagggcaaagccctccctgaataccgctctacccga                                                                                                                                                                                                                 <  2:43468‑M1/298‑50 (MQ=255)                                                 cGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgc                                                                                                                                                                                                              >  2:171956‑M1/1‑248 (MQ=255)                                                   ggCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcc                                                                                                                                                                                                           <  1:334334‑M1/301‑56 (MQ=255)                                                          aCGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGTTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccac                                                                                                                                                                                                        <  1:105732‑M1/297‑59 (MQ=255)                                                            gCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccccccccacagggc                                                                                                                                                                                                   >  2:138873‑M1/1‑237 (MQ=255)                                                                      gTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgcccccacagggccacgctggc                                                                                                                                                                                          >  1:86828‑M1/1‑227 (MQ=255)                                                                       tGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaa                                                                                                                                                                                            >  1:42467‑M1/1‑226 (MQ=255)                                                                        gACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCGTCCTGATCCTTGAGACTGCCCAATCATTAgcggccgcgattggagggcaaagccctcactgaataccgcactacccgac                                                                                                                                                                                                                >  2:143797‑M1/1‑225 (MQ=255)                                                                            ccGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaag                                                                                                                                                                                     <  1:171962‑M1/298‑78 (MQ=255)                                                                                 ggCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaagg                                                                                                                                                                                                                  <  2:74258‑M1/264‑49 (MQ=255)                                                                                 ggCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaagg                                                                                                                                                                                                                  >  1:74262‑M1/1‑216 (MQ=255)                                                                                         tGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACAACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacccaaggccgctccc                                                                                                                                                                                                          >  2:194331‑M1/1‑208 (MQ=255)                                                                                         tGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACAACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctccc                                                                                                                                                                                                          <  1:194338‑M1/264‑57 (MQ=255)                                                                                          ggCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtg                                                                                                                                                                    <  1:138878‑M1/301‑95 (MQ=255)                                                                                                gaataaaacgccGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaaggg                                                                                                                                                                                   <  1:420677‑M1/276‑80 (MQ=255)                                                                                                gaagaaAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacaccaggccgctcccacagggccacgct                                                                                                                                                                                             >  2:420665‑M1/1‑201 (MQ=255)                                                                                                 aagaaAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctcca                                                                                                                                                              <  2:42464‑M1/300‑101 (MQ=255)                                                                                                                   tCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggctcgttgcgcccctttcgcgacacaaggccgcccccacagggccacgctggcacaagggaatgcagtccctgtgggagcggcct                                                                                                                                                          >  1:308900‑M1/1‑182 (MQ=255)                                                                                                                                     aCGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGGGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgc                                                                                                                                                                              <  2:86824‑M1/248‑85 (MQ=255)                                                                                                                                              ccGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggc                                                                                                                                                                                                <  1:143802‑M1/221‑67 (MQ=255)                                                                                                                                                    ccTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggc                                                                                                                                                                                          >  2:30777‑M1/1‑149 (MQ=255)                                                                                                                                                    ccTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggc                                                                                                                                                                                          <  1:30780‑M1/221‑73 (MQ=255)                                                                                                                                                                                        acaacGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactcaacggcttcagccaggtgc                                                                          >  2:323569‑M1/1‑113 (MQ=255)                                                                                                                                                                                                   ggctacTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactcaacggcttcagccaggtgctgctgctgcgcaac                                                            <  1:323577‑M1/300‑199 (MQ=255)                                                                                                                                                                                                               gcgcTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCAAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggcct                                                                                                                                                          <  2:308892‑M1/194‑105 (MQ=255)                                                                                                                                                                                                                                    accgaaccgaCCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccctgactgggccacggcgttactcaacggcttcagccaggtgctgctgctgcgcaaccccctgtgcggcctgtgctgcctgctggccatcc                          >  1:84574‑M1/1‑69 (MQ=255)                                                                                                                                                                                                                                            gaCCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactcaacggcttcagccaggtgctgctgctgcgcaaccccctgtgc                                                   >  2:369084‑M1/1‑61 (MQ=255)                                                                                                                                                                                                                                            gaCCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactcaacggcttcagccaggtgctgctgctgcgcaaccccctgtgc                                                   <  1:369093‑M1/268‑208 (MQ=255)                                                                                                                                                                                                                                                               gACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactcaacggcttcagccaggtgctgctgctgcgcaaccccctgtgcggcctgttctgcctgctggccatcct                         >  2:39333‑M1/1‑42 (MQ=255)                                                                                                                                                                                                                                                               gACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactcaacggcttcagccaggtgctgctgctgcgcaaccccctgtgcggcctgtgctgcctgctggccatcctgctgaccgcacccaac         <  1:39336‑M1/291‑250 (MQ=255)                                                                                                                                                                                                                                                                                                     aTGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgcgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccacagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggc  >  2:321525‑M1/1‑4 (MQ=255)                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                 GTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGACGCCTGATAAAAAAGTTAAGATCCTTGCCACCCTTGGCCCTGCGATCAACGGCATCGACGATATCCGCCAGCTGGTCGAAGCCGGGGTGAACATCTTCCGCCTCAACTTCAGCCACGGCGAACATGCCGACCACGCCCTGCGCTACCAGTGGATCCGCGAAGTCGAGCAACAGCTGAACTACCCGCTGGGCATCCTCATGGACCTGCAAGGGCCGAAGCTGCGCGTTGGCCGCTTCGCCGAAGGCAAGGTGCAGTT  >  AE015451/4891818‑4892369

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