Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	AE015451	5,342,745	1681 bp→GCGGCCGC		[pgi‑II]	[pgi‑II]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	AE015451	5342745	5344425	1681	54 [0]	[0] 52	[pgi‑II]	[pgi‑II]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 5342744	0 (0.000)	52 (0.760) +GCGGCCGC	43/584	0.4	100%	intergenic (+190/‑17)	panC/pgi‑II	pantothenate synthetase/glucose‑6‑phosphate isomerase
	?	AE015451	5344426 =	0 (0.000)					intergenic (+1/‑177)	pgi‑II/acsA‑II	glucose‑6‑phosphate isomerase/acetyl‑CoA synthetase

GCAGTGATCAATGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCAAGGGAACCCCGCGCAATGGCGTATTACCGTACACCCCACGATGTGACGGCCCTGCCCGCCTGGCAGGCGCTTCAGAAACACCGCGACGCCATGCAAAGCTTCAGCATGCGCGAAGCCTTCGCCGCCGATGCCAAGCGCTTCGACCAGTTTTCCCTCAGCGCTTGCGGCCTGTTCCTCGACTACTCGAAAAACCTGATCACCGAGCAAAGCCGCGATCTGCTGGTCAATCTGGCCAATGAAGTGGGCCTGCAGGACGCCATCAAGTCGATGTTCAGCGGTGAGATCAT  >  AE015451/5342450‑5343032                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                 gCAGTGATCAATGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCT‑GGCTTGCTTGCACTGCCCCCCAACTCCCCCTGCCATACCCATTCCCCGCACGTTGCCTTTGCCTATAATGGTGCCAGCCTGACCCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCCCCgccgcc                                                                                                                                                                                                                                                                                            >  4:106222‑M1/1‑295 (MQ=255) gCAGTGATCAATGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggcc                                                                                                                                                                                                                                                                                            >  4:636492‑M1/1‑295 (MQ=255)           aTGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagct                                                                                                                                                                                                                                                                                  >  3:344221‑M1/1‑285 (MQ=255)             cccctttccctaaacgacgctcttccgatctacctgggtaacacgcggttgatcgacaaccTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttg                                                                                                                                                                                                                                                                                <  2:611620‑M1/271‑19 (MQ=255)                          tggtgATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttc                                                                                                                                                                                                                                                                   <  4:186510‑M1/301‑32 (MQ=255)                            gtgATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATG‑CCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcg                                                                                                                                                                                                                                                                >  3:118460‑M1/1‑267 (MQ=255)                            gctcttccgatctcttacctgggtaacACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgc                                                                                                                                                                                                                                                                 <  1:848‑M1/288‑34 (MQ=255)                              gATTGCGGCGGCTTACCTGGGTAAAACGCGGTTGATCGAAAACCTCTACCTGAATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATG‑CCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcggg                                                                                                                                                                                                                                                              <  1:118458‑M1/301‑37 (MQ=255)                                 gccggcggctTCCCTGGGTAACACGCGGGTGACCGACAACCTCTACCTGCGTTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATGGGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTGGCCTATAATGGTGCAAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggta                                                                                                                                                                                                                                                            <  2:106198‑M1/299‑39 (MQ=255)                                       ggcTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgccccgagcttgtgccggcctcttcgcgggtaaacccc                                                                                                                                                                                                                                                      >  2:50916‑M1/1‑257 (MQ=255)                                         cTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcagatcggaagagc                                                                                                                                                                                                                                                    >  3:848‑M1/1‑255 (MQ=255)                                                  gTAACACGCGGTGGATCGACAACCTCTACCTGCATTGGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggca                                                                                                                                                                                                                                           <  1:344220‑M1/301‑56 (MQ=255)                                                              ttccgatctacctctacctGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctc                                                                                                                                                                                                                               <  4:254515‑M1/292‑68 (MQ=255)                                                              ctggagttcagacgtgtgctcttccgatctagaagacagcataacggtcttgggctgccttgcagcCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctc                                                                                                                                                                                                                               <  3:36603‑M1/271‑68 (MQ=255)                                                                  cgcgaaccgctgcctgaagttggaagagaagacgtcatgaCGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtag                                                                                                                                                                                                                           <  2:636289‑M1/282‑72 (MQ=255)                                                                   gACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtagg                                                                                                                                                                                                                          <  4:50933‑M1/301‑73 (MQ=255)                                                                   gACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgagttacctcgtagg                                                                                                                                                                                                                          >  4:345641‑M1/1‑229 (MQ=255)                                                                       ccctCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcg                                                                                                                                                                                                                      <  3:226264‑M1/300‑77 (MQ=255)                                                                       aCCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctccgatcggag                                                                                                                                                                                                                      >  2:254457‑M1/1‑225 (MQ=255)                                                                          tctACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggg                                                                                                                                                                                                                   >  4:191517‑M1/1‑222 (MQ=255)                                                                                tctatttGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttaccc                                                                                                                                                                                                             <  1:378550‑M1/298‑86 (MQ=255)                                                                                   aTTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccaga                                                                                                                                                                                                          >  3:378550‑M1/1‑213 (MQ=255)                                                                                      tGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggccccttcgcgggtaaacccgcctacccggcaacgcggtacctcgtaggagcgggcttacccgccaaa                                                                                                                                                                                                       >  3:280244‑M1/1‑210 (MQ=255)                                                                                        gagttcagacgtgtgctcttccgatctgctgccttgcagcccattcgcgggcaagCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagag                                                                                                                                                                                                     <  3:139486‑M1/274‑94 (MQ=255)                                                                                            agaagaCAGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggccacgcgttacctcagatcggaagagcgtcgtgtagggaaagag                                                                                                                                                                                                 >  1:36602‑M1/1‑204 (MQ=255)                                                                                                      aTAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttt                                                                                                                                                                                       >  1:89139‑M1/1‑194 (MQ=255)                                                                                                                   gCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagagagatcggaagcgcgtcgggtagggaaa                                                                                                                                                                          >  1:139484‑M1/1‑181 (MQ=255)                                                                                                                    gggccttgacgcccattcgcgGGCAAGCCCGTGCCTGCACCTTTACCTGTAGGGGCT‑GGCTTGCTTGCAATGCCCCCAAAACCCCCCTGCCCTACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcgtcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccag                                                                                                                                                                         <  1:280245‑M1/291‑122 (MQ=255)                                                                                                                    ccgcatgcactcccagtcccgggacagcccgctcctgcacgtttaccggtaCGGGCG‑GTCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAACGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccag                                                                                                                                                                         <  2:345557‑M1/276‑122 (MQ=255)                                                                                                                                             aGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggatcgggtttacccgctaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgt                                                                                                                                                >  1:285155‑M1/1‑155 (MQ=255)                                                                                                                                               cccGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggccccttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactggattacccaatagccgttg                                                                                                                                              >  1:416554‑M1/1‑153 (MQ=255)                                                                                                                                                  gCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcag                                                                                                                                           >  4:425985‑M1/1‑150 (MQ=255)                                                                                                                                                     ccTGCACGTTATCCTGTA‑GGTCGCGGCTCGCTTGCAATGCCCCCAAACTCCCCCTGCCATACCCATTCCCAGCACCTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtcc                                                                                                                                        <  2:425876‑M1/301‑155 (MQ=255)                                                                                                                                                        gCACGTTTCCCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatc                                                                                                                                     <  3:416556‑M1/301‑158 (MQ=255)                                                                                                                                                        cgccgtctacctgtcGGGGCG‑GTCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatc                                                                                                                                     <  2:191471‑M1/294‑158 (MQ=255)                                                                                                                                                         cACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcggttaaccccgcctacaaggcaacgcggtacctcgtcggaccgggtttccccccgcacaggccggcccaggctgaacaccgctcccagctacactgtcttatccattagcccttcccgtccctct                                                                                                                                    >  2:69960‑M1/1‑143 (MQ=255)                                                                                                                                                                          cttccctccacgcccctcttccgatctctccccctcccctaccccttcccagcacgtGGCCTTTGCCTATACTGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCACGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcag                                                                                                                    <  1:364882‑M1/273‑175 (MQ=255)                                                                                                                                                                                 ttgcttgcAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtaccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacc                                                                                                             <  4:764084‑M1/301‑182 (MQ=255)                                                                                                                                                                                   gcttgcAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggtgcggggttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtattatcgaatagccgttgcaatccatcgcccctttcacgagcagaaccccccc                                                                                                           >  2:409763‑M1/1‑118 (MQ=255)                                                                                                                                                                                       ccctacccgacgctcttccgatctcatacccattcccagcacgtggcctTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccat                                                                                                       <  1:59640‑M1/277‑188 (MQ=255)                                                                                                                                                                                           cgcCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccatgttc                                                                                                   <  3:89142‑M1/300‑192 (MQ=255)                                                                                                                                                                                            gCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccatgttcg                                                                                                  <  4:69979‑M1/301‑193 (MQ=255)                                                                                                                                                                                                     aTCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcagagatcggaagagcacacgtctgaactc                                                                                         >  3:364883‑M1/1‑100 (MQ=255)                                                                                                                                                                                                          ccTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagacgccggcccaggttgaacaccgctcccagctacactgtcttatcgaataaccgttgcagtccatcgccccttacaagagcaggaccacccgccatgttcgatatccgcaaataa                                                                                    >  1:127017‑M1/1‑95 (MQ=255)                                                                                                                                                                                                            tGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccatgttcgatatccgcaaataccc                                                                                  >  3:503117‑M1/1‑93 (MQ=255)                                                                                                                                                                                                             gCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccatgttcgatatccgcaaatacccc                                                                                 >  3:376132‑M1/1‑92 (MQ=255)                                                                                                                                                                                                               cATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccatagatcggaagagcacacgtctgaa                                                                               >  3:59642‑M1/1‑90 (MQ=255)                                                                                                                                                                                                                    ccATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccatgttcgatatccgccaatacccccaggccc                                                                          >  4:664818‑M1/1‑85 (MQ=255)                                                                                                                                                                                                                          cccgcacGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGCCGGCCGTGTCCAAGCCCCCACCCAGCACCgcggccgctccgagcttgggccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccatgttcgatatccgcaaatacccccaggccctggccg                                                                    <  2:664596‑M1/298‑223 (MQ=255)                                                                                                                                                                                                                                             tataATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccatgttcgatatccgcaaatacccccaggccctggccgtcagccagtccgccgcccc                                                 >  1:372808‑M1/1‑60 (MQ=255)                                                                                                                                                                                                                                                                                  tgtCCAAGCCCTCGCCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccatgttcgatatccgcaaatacccccaggccctggccgtcagccagtccgccgccctcacccccgaagactaccgccgcctctaccgccagtcg            <  3:372808‑M1/301‑279 (MQ=255)                                                                                                                                                                                                                                                                                      cAAGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccatgttcgatatccgcaaatacccccaggccctggccgtcagccagtccgccgccctcacccccgaagcctaccgccgcctctaccgccagtcagtcc        >  1:200159‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                                                                                        aGCCCTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccatgttcgatatccgcaaatacccccaggccctggccgtcagccagtccgccgccctcacccccgaacactaccgccgcctctcccgccagtcggtcgaa      >  3:197162‑M1/1‑17 (MQ=255)                                                                                                                                                                                                                                                                                            cTCACCACGCACCgcggccgctccgagcttgtgccggcctcttcgcgggtaaacccgcctaccaggcaacgcgttacctcgtaggagcgggtttacccgcgaagaggccggcccaggttgaacaccgctcccagctacactgtcttatcgaatagccgttgcagtccatcgccccttacaagagcaggaccacccgccatgttcgatatccgcaaatacccccaggccctggccgtcagccagtccgccgccctcacccccgaagactaccgccgcctctaccgccagtcggtcgaagacc  <  4:409872‑M1/301‑289 (MQ=255)                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                 GCAGTGATCAATGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCG‑GGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCAAGGGAACCCCGCGCAATGGCGTATTACCGTACACCCCACGATGTGACGGCCCTGCCCGCCTGGCAGGCGCTTCAGAAACACCGCGACGCCATGCAAAGCTTCAGCATGCGCGAAGCCTTCGCCGCCGATGCCAAGCGCTTCGACCAGTTTTCCCTCAGCGCTTGCGGCCTGTTCCTCGACTACTCGAAAAACCTGATCACCGAGCAAAGCCGCGATCTGCTGGTCAATCTGGCCAATGAAGTGGGCCTGCAGGACGCCATCAAGTCGATGTTCAGCGGTGAGATCAT  >  AE015451/5342450‑5343032

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