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

CGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGACGCCTGATAAAAAAGTTAAGATCCTTGCCACCCTTGGCCCTGCGATCAACGGCATCGACGATATCCGCCAGCTGGTCGAAGCCGGGGTGAACATCTTCCGCCTCAACTTCAGCCACGGCGAACATGCCGACCACGCCCTGCGCTACCAGTGGATCCGCGAAGTCGAGCAACAGCTGAACTACCCGCTGGGCATCCTCATGGACCTGCAAGGGCCGAAGCTGCGCGTTGGCCGCTTCGCCGAAGGCAAGGTGCAGTTGCAACGCGGGCAGGCCCTGCGCCTGGATCTGGAC  >  AE015451/4891822‑4892403                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                   cGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgc                                                                                                                                                                                                                                                                                            <  1:925962‑M1/300‑9 (MQ=255)      aGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattg                                                                                                                                                                                                                                                                                       >  1:736961‑M1/1‑287 (MQ=255)          ccTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccggg                                                                                                                                                                                                                                                                                  <  1:961667‑M1/301‑19 (MQ=255)                      cACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggg                                                                                                                                                                                                                                                                                 <  2:556066‑M1/290‑20 (MQ=255)                      cACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggg                                                                                                                                                                                                                                                                                 >  1:556087‑M1/1‑271 (MQ=255)                          gAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccc                                                                                                                                                                                                                                                                   >  1:940916‑M1/1‑267 (MQ=255)                                  gAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgac                                                                                                                                                                                                                                                          >  2:563085‑M1/1‑259 (MQ=255)                                  gAAAGGCCGGCCGGGCGTGTACGCCCTGGCCGGGGAAACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcg                                                                                                                                                                                                                                                              <  2:105954‑M1/297‑39 (MQ=255)                                                     tACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccggggacgcgttgcgcccctttc                                                                                                                                                                                                                                                               >  2:424041‑M1/1‑240 (MQ=255)                                                         cccTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggcc                                                                                                                                                                                                                                    <  2:736935‑M1/300‑65 (MQ=255)                                                                 ggTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATTAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacccaaggccgctcccacagggcc                                                                                                                                                                                                                                    >  1:820962‑M1/1‑228 (MQ=255)                                                                    gACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcac                                                                                                                                                                                                                          >  2:675225‑M1/1‑225 (MQ=255)                                                                    gACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcgg                                                                                                                                                                                                                                                                                                >  2:894587‑M1/1‑225 (MQ=255)                                                                    gACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcgg                                                                                                                                                                                                                                                                                                <  1:894613‑M1/229‑5 (MQ=255)                                                                         cGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaaggg                                                                                                                                                                                                                     >  1:391265‑M1/1‑220 (MQ=255)                                                                          gACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgtttcgcccctttcgccacacaaggcc                                                                                                                                                                                                                                                  >  1:476797‑M1/1‑219 (MQ=255)                                                                          gACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaa                                                                                                                                                                                                                   <  1:563106‑M1/300‑82 (MQ=255)                                                                            cGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaaggga                                                                                                                                                                                                                    <  2:940890‑M1/297‑81 (MQ=255)                                                                                aTCGATGGCTCGGAAGAAAACGCCGGTGCCGTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgca                                                                                                                                                                                                               <  2:622473‑M1/298‑86 (MQ=255)                                                                                      ggCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgt                                                                                                                                                                                                       <  1:414696‑M1/300‑94 (MQ=255)                                                                                          cGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctg                                                                                                                                                                                                        <  2:967799‑M1/295‑93 (MQ=255)                                                                                           ggAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCTGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgtgggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaaggga                                                                                                                                                                                                                    <  2:561639‑M1/282‑81 (MQ=255)                                                                                           ggAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgccggc                                                                                                                                                                                                                            >  2:20602‑M1/1‑202 (MQ=255)                                                                                           ggAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCACAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgatggcacaagggc                                                                                                                                                                                                                    >  1:561660‑M1/1‑202 (MQ=255)                                                                                           ggAAGAAAACGCCGGTGCCTTCATGACCCCGGCAAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtggga                                                                                                                                                                                                   <  2:820936‑M1/299‑98 (MQ=255)                                                                                             aagaaAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGCTCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCCTCCTGATCCTTGAGCCTGCCCAATCATTAgcggccgc                                                                                                                                                                                                                                                                                            >  2:232670‑M1/1‑200 (MQ=255)                                                                                             aagaaAACGCCGGTCCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcgg                                                                                                                                                                                               <  2:264638‑M1/301‑102 (MQ=255)                                                                                                      ccGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcg                                                                                                                                                                                                                                                              >  1:329444‑M1/1‑191 (MQ=255)                                                                                                      ccGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcg                                                                                                                                                                                                                                                              <  2:329433‑M1/229‑39 (MQ=255)                                                                                                           gCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCCGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgc                                                                                                                                                                                    >  1:737903‑M1/1‑186 (MQ=255)                                                                                                              ttCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggc                                                                                                                                                                                              >  1:182129‑M1/1‑183 (MQ=255)                                                                                                                   gACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCTTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtcc                                                                                                                                                                                                           <  2:373712‑M1/267‑90 (MQ=255)                                                                                                                   gACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggc                                                                                                                                                                                              <  2:182125‑M1/280‑103 (MQ=255)                                                                                                                           ccAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgt                                                                                                                                                                                                       <  1:162397‑M1/263‑94 (MQ=255)                                                                                                                           ccAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgt                                                                                                                                                                                                       >  2:162393‑M1/1‑170 (MQ=255)                                                                                                                              gCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacccaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgg                                                                                                                                                                                                     >  1:418738‑M1/1‑167 (MQ=255)                                                                                                                              gCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgg                                                                                                                                                                                                     <  2:418718‑M1/262‑96 (MQ=255)                                                                                                                               cTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccc                                                                                                                                                                  >  1:969520‑M1/1‑166 (MQ=255)                                                                                                                               cTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgt                                                                                                                                                                                         >  1:66582‑M1/1‑166 (MQ=255)                                                                                                                                  cGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcg                                                                                                                                                                                                                                                                              >  2:922492‑M1/1‑163 (MQ=255)                                                                                                                                   gCCAGCGCCGAAGCCTTGGGCCTGTCAGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccc                                                                                                                                                                                                          <  2:795334‑M1/252‑91 (MQ=255)                                                                                                                                        cgcCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgt                                                                                                                                                                                         <  2:66580‑M1/264‑108 (MQ=255)                                                                                                                                         gcCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttg                                                                                                                                                                                          <  2:391253‑M1/262‑107 (MQ=255)                                                                                                                                                    gggCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgt                                                                                                                                                                                         <  2:476776‑M1/252‑108 (MQ=255)                                                                                                                                                                            ggACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggccca                                                                                                                       <  2:737877‑M1/294‑174 (MQ=255)                                                                                                                                                                                        cggctacggctaCTTCGCCGCTCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccc                                                                                                    <  1:675250‑M1/301‑193 (MQ=255)                                                                                                                                                                                                                                 ccgaaccgaCCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcaaaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaacgcggaaaaccgaggcccaaagaggcccaccg                                                                                                          >  1:351748‑M1/1‑68 (MQ=255)                                                                                                                                                                                                                                  cgaaccgaCCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgccggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttc                                                           <  1:20604‑M1/300‑234 (MQ=255)                                                                                                                                                                                                                                                        aaCGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactggg                                     >  2:581656‑M1/1‑45 (MQ=255)                                                                                                                                                                                                                                                               tCCGCGCCGGCCTGATCCTTGCGACTGCCCAATAATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaa                                                                                                                                                          <  1:922518‑M1/176‑139 (MQ=255)                                                                                                                                                                                                                                                                                     gaCTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccg                                                                                                          <  2:351737‑M1/202‑187 (MQ=255)                                                                                                                                                                                                                                                                                          cccAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgtgagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggc                                                                                                                          <  1:122968‑M1/181‑171 (MQ=255)                                                                                                                                                                                                                                                                                            cAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacaggcccacgctggcacaagggaatgcagtccctgtgggagcgcccttgtgtcgcggttggagggcaaatccctccc                                                                                                                                                              >  2:677512‑M1/1‑9 (MQ=255)                                                                                                                                                                                                                                                                                            cAATCATCAgcggccgcgattgccgggggcgctttgcgcgcctttcgcgacacaaggccgctcccacagggcaacgctggcacaagtgaatgcagtccctgtggtagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactcaacggcttcagccaggtgctg  <  2:969494‑M1/299‑291 (MQ=255)                                                                                                                                                                                                                                                                                                    agcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaacgaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttc                                                           >  2:787266‑M1/1‑1 (MQ=255)                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                   CGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGACGCCTGATAAAAAAGTTAAGATCCTTGCCACCCTTGGCCCTGCGATCAACGGCATCGACGATATCCGCCAGCTGGTCGAAGCCGGGGTGAACATCTTCCGCCTCAACTTCAGCCACGGCGAACATGCCGACCACGCCCTGCGCTACCAGTGGATCCGCGAAGTCGAGCAACAGCTGAACTACCCGCTGGGCATCCTCATGGACCTGCAAGGGCCGAAGCTGCGCGTTGGCCGCTTCGCCGAAGGCAAGGTGCAGTTGCAACGCGGGCAGGCCCTGCGCCTGGATCTGGAC  >  AE015451/4891822‑4892403

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: ‑