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] 38	pyk	Pyruvate kinase

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

TAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGACGCCTGATAAAAAAGTTAAGATCCTTGCCACCCTTGGCCCTGCGATCAACGGCATCGACGATATCCGCCAGCTGGTCGAAGCCGGGGTGAACATCTTCCGCCTCAACTTCAGCCACGGCGAACATGCCGACCACGCCCTGCGCTACCAGTGGATCCGCGAAGTCGAGCAACAGCTGAACTACCCGCTGGGCATCCTCATGGACCTGCAAGGGCCGAAGCTGCGCGTTGGCCGCTTCGCCGAAGGCAAGGTGCAGTTGCAACGCGGGCAGGCCCTGCGCC  >  AE015451/4891819‑4892392                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                        ctgtccctaccggacgctattccgttcggaaagcctgaaaggcctgccgggcgtgtaCGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggcc                                                                                                                                                                                                                                                                                   <  2:330453‑M1/273‑7 (MQ=255)      tctttccctacacgacgctcttccgatctctgaaaggcctgccgggcgtgtacgccctgGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgat                                                                                                                                                                                                                                                                              <  1:42065‑M1/272‑12 (MQ=255)             ccTGCTCAGCCTAACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccggg                                                                                                                                                                                                                                                                       <  2:751185‑M1/301‑19 (MQ=255)                             gAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccagatcggaagagcacacgtctgaactcc                                                                                                                                                                                                                                                       >  4:330532‑M1/1‑267 (MQ=255)                               gtttcggaaaggcCTGCCGGGGGCGTACGGCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgccccttt                                                                                                                                                                                                                                                     <  2:269308‑M1/295‑37 (MQ=255)                                aGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTTACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttc                                                                                                                                                                                                                                                    <  4:530993‑M1/301‑38 (MQ=255)                                aGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttc                                                                                                                                                                                                                                                    <  4:406255‑M1/301‑38 (MQ=255)                                  ccTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcccccctttcgc                                                                                                                                                                                                                                                  >  1:227973‑M1/1‑262 (MQ=255)                                   cTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttagat                                                                                                                                                                                                                                                 >  4:269374‑M1/1‑261 (MQ=255)                                   cTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgatagatcggaagagcacacgtctgaactcca                                                                                                                                                                                                                                                 >  3:42066‑M1/1‑261 (MQ=255)                                                       cttccgatctgcaggaatcgcggacgcagtcgagggctcggaagaaaacgccgatgccgTCATGACCCCTGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacag                                                                                                                                                                                                                             <  1:333266‑M1/272‑61 (MQ=255)                                                                 gCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgccacaagcacgctcccccagagatagggag                                                                                                                                                                                                                   >  3:333267‑M1/1‑231 (MQ=255)                                                                    ggTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggca                                                                                                                                                                                                                <  4:347385‑M1/301‑74 (MQ=255)                                                                                                aagaaAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcgg                                                                                                                                                                                    <  3:469393‑M1/301‑102 (MQ=255)                                                                                                         ccGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtc                                                                                                                                                                           >  4:58678‑M1/1‑191 (MQ=255)                                                                                                                 ttCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccccgctggcccaagggaatgcagcccctgtgcgaccggcctggtgtccggattgg                                                                                                                                                                   >  2:221457‑M1/1‑183 (MQ=255)                                                                                                                     tGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcggtgcgcccctttcgcgacacaaggccgcccccccagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgttgcgagtgggggg                                                                                                                                                               >  3:116948‑M1/1‑179 (MQ=255)                                                                                                                        ccccGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcgtcgctctgtcgagcggagggggaca                                                                                                                                                            >  3:125307‑M1/1‑176 (MQ=255)                                                                                                                          ccGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaag                                                                                                                                                          <  1:116946‑M1/301‑128 (MQ=255)                                                                                                                                   gttcagacgtgtgctcttccgatcttgtcggccagcgacgagctggacaacAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctg                                                                                                                                                 <  3:144774‑M1/276‑137 (MQ=255)                                                                                                                                                  gCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctaccc                                                                                                                                  <  4:221507‑M1/301‑152 (MQ=255)                                                                                                                                                           cTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcgga                                                                                                                         >  4:28995‑M1/1‑141 (MQ=255)                                                                                                                                                            tGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgagatcggaagagcgtcgggtagggg                                                                                                                        >  1:144772‑M1/1‑140 (MQ=255)                                                                                                                                                             gTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaca                                                                                                                       >  4:142406‑M1/1‑139 (MQ=255)                                                                                                                                                                  ccAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctacctgaataccgctctacccgaccgcggaacaccga                                                                                                                  >  1:12874‑M1/1‑134 (MQ=255)                                                                                                                                                                            ctggaaaaacaac‑acggccacggCTACTTCGCGGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagag                                                                                                       <  2:58661‑M1/294‑179 (MQ=255)                                                                                                                                                                             cTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagag                                                                                                       <  3:12874‑M1/301‑179 (MQ=255)                                                                                                                                                                                          ccggctacggctacTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctccc                                                                                          <  3:227971‑M1/300‑192 (MQ=255)                                                                                                                                                                                                                            gATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccc                                                        >  1:359765‑M1/1‑76 (MQ=255)                                                                                                                                                                                                                                  cACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatt                                                  <  3:359766‑M1/301‑232 (MQ=255)                                                                                                                                                                                                                                     cgaaccgaCCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgccaccccatcaacctccactgactgccccatgtacaccaaaaatttct                                               >  3:280800‑M1/1‑67 (MQ=255)                                                                                                                                                                                                                                      aacaccacccacgacaagctcaacgacctccgcgcaatcCTGATCCTTGAGACGGCCAGAACACGAgcggccgcgcttgccgggggcgcgttgctcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgt                                              <  2:142377‑M1/282‑236 (MQ=255)                                                                                                                                                                                                                                             ccgcaaccaactgcaccgacttccgcgccatccTGTTCCTTGATACTGCCCAAGCATCGgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccg                                       <  1:280801‑M1/285‑243 (MQ=255)                                                                                                                                                                                                                                               cGCACCAACGTCAACGACTTCGGGGCAATCCTGATCCTGGAGACTGCCCAATCATGGgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtg                                     <  2:28985‑M1/301‑245 (MQ=255)                                                                                                                                                                                                                                                   ccgcgcttttccgacctcgcgccagccttctccttgagACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgcccc                                 <  2:748169‑M1/284‑249 (MQ=255)                                                                                                                                                                                                                                                     aaCGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgc                               >  3:218115‑M1/1‑51 (MQ=255)                                                                                                                                                                                                                                                              gACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgcacccccatcaacctccactgcctgccccatgtacaccaaaaatttcgtcaacccgtgccccgaccgggccac                      >  1:157565‑M1/1‑42 (MQ=255)                                                                                                                                                                                                                                                                    cgcgCCATCCTGATCCTTGAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccagatcggaagagcacac                >  4:748448‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                                                                             cTGATCCTTAAGACTGCCCAATCATGAgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgtcccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactcaacgg       >  3:98681‑M1/1‑27 (MQ=255)                                                                                                                                                                                                                                                                                  ccTTGAGACTGCCCAATCATGGgcggccgcgattgccgggggcgcgttgcgcccctttcgcgacacaaggccgctcccacagggccacgctggcacaagggaatgcagtccctgtgggagcggccttgtgtcgcgattggagggcaaagccctccctgaataccgctctacccgaccgcggaaaaccgaggcccaaagaggcccaccgctcccccatcaacctccactgactgccccatgtacaccaaaaatttcgtcaacccgtgccccgactgggccacggcgttactcaacggcttca  <  3:157567‑M1/301‑280 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                        TAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCCGCGCTGGATGCGCTGATCGTCACCGAACCGACCCGCACCAACGTCAACGACTTCCGCGCCATCCTGATCCTTGAGACTGCCCAATCATGACGCCTGATAAAAAAGTTAAGATCCTTGCCACCCTTGGCCCTGCGATCAACGGCATCGACGATATCCGCCAGCTGGTCGAAGCCGGGGTGAACATCTTCCGCCTCAACTTCAGCCACGGCGAACATGCCGACCACGCCCTGCGCTACCAGTGGATCCGCGAAGTCGAGCAACAGCTGAACTACCCGCTGGGCATCCTCATGGACCTGCAAGGGCCGAAGCTGCGCGTTGGCCGCTTCGCCGAAGGCAAGGTGCAGTTGCAACGCGGGCAGGCCCTGCGCC  >  AE015451/4891819‑4892392

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