Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	AE015451	2,032,505	1680 bp→GCGGCCGC		[pgi‑I]	[pgi‑I]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	AE015451	2032505	2034184	1680	63 [1]	[0] 57	[pgi‑I]	[pgi‑I]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 2032504	1 (0.010)	57 (0.830) +GCGGCCGC	48/584	0.2	99.2%	intergenic (‑260/‑16)	kdsA‑II/pgi‑I	2‑dehydro‑3‑deoxyphosphooctonate aldolase/glucose‑6‑phosphate isomerase
	?	AE015451	2034185 =	0 (0.000)					intergenic (+1/+1440)	pgi‑I/PP_5477	glucose‑6‑phosphate isomerase/conserved protein of unknown function

ACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACATTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAACTGGGCAAGAGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAATTACTTTCGTGGTAAACACCGTGGCTGACCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGTTGTGACTAGAGGATCTGAATGTACAATGCACTTCTGCTAATCGCTGCAATTGCCCTTACCGGATGCGC  >  AE015451/2033889‑2034444                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                    attcgattgtggagttgatcttgatcattttttactcgataggtccttagggggcaacccgcaggacatggcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTg                                                                                                                                                                                                                                                                 >  1:116610‑M2/297‑301 (MQ=255)    cgattgtggagttgatcttgatcattttttactcgataggtccttagggggcaacccgcaggacatggcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCtt                                                                                                                                                                                                                                                              >  3:405451‑M2/294‑301 (MQ=255)     gattgtggagttgatcttgatcattttttactcgataggtccttagggggcaacccgcaggacatggcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTg                                                                                                                                                                                                                                                             >  4:564281‑M2/293‑301 (MQ=255)                     ttgatcattttttactcgataggtccttagggggcaacccgcaggacatggcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCggg                                                                                                                                                                                                                                             >  3:204409‑M2/277‑301 (MQ=255)                      tgatcattttttactcgataggtccttagggggcaacccgcaggacatggcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtaacgcggccgcCCCTGCTTGATACTGGCCCGGcgggg                                                                                                                                                                                                                                            >  2:585483‑M2/276‑300 (MQ=255)                               tttactcgataggtccttagggggcaacccgcaggacatggcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTcc                                                                                                                                                                                                                                   <  1:405450‑M2/35‑1 (MQ=255)                               tttactcgataggtccgtagggggcaacccgcaggacatggcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTcc                                                                                                                                                                                                                                   <  3:116612‑M2/35‑1 (MQ=255)                                     cgataggtccttagggggcaacccgcaggacatggcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGTCCGGCGGGAGGAGCGTCCGAGCgg                                                                                                                                                                                                                             >  4:547575‑M2/261‑301 (MQ=255)                                            tccttagggggcaacccgcaggacatggcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCtggt                                                                                                                                                                                                                      <  1:204411‑M2/48‑1 (MQ=255)                                                tagggggcaacccgcaggacatggcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCtggtggtg                                                                                                                                                                                                                  >  3:158899‑M2/250‑300 (MQ=255)                                                          cccgcaggacatggcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAaca                                                                                                                                                                                                        >  1:618572‑M2/240‑301 (MQ=255)                                                           ccgcaggacatggcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAacac                                                                                                                                                                                                       <  3:618574‑M2/63‑1 (MQ=255)                                                                       gcaaagggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATg                                                                                                                                                                                           >  2:569111‑M2/227‑301 (MQ=255)                                                                          ggtgctcttccgatctcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCg                                                                                                                                                                                        <  4:567585‑M2/78‑1 (MQ=255)                                                                          cttccgatctagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCg                                                                                                                                                                                        <  1:449644‑M2/78‑1 (MQ=255)                                                                           agggataatagtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttcccgttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGGGGAGTGTCCGGGCGGAGCTGGTGGTCGCCGGCAGCACTAGCGCCGCATggctg                                                                                                                                                                                       >  2:347953‑M2/223‑299 (MQ=255)                                                                                    agtcagcgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCgcatcatggcgagatcggaag                                                                                                                                                                              >  3:449646‑M2/214‑291 (MQ=255)                                                                                      tcagcgcgcccttgctaaccggggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAgg                                                                                                                                                                            <  2:564111‑M2/90‑1 (MQ=255)                                                                                          cgcgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACactagcgcatcagggcgaggtcggaaagggcgc                                                                                                                                                                        >  2:567415‑M2/208‑285 (MQ=255)                                                                                            cgcccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCAt                                                                                                                                                                      <  4:585668‑M2/96‑1 (MQ=255)                                                                                              cccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCAttt                                                                                                                                                                    <  4:569281‑M2/98‑1 (MQ=255)                                                                                              cccttgctaaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCAttt                                                                                                                                                                    <  3:111324‑M2/98‑1 (MQ=255)                                                                                               gtgtgctcttccgatctgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTc                                                                                                                                                                   <  4:711200‑M2/99‑1 (MQ=255)                                                                                                  tgctaaccggggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGAc                                                                                                                                                                <  2:55414‑M2/102‑1 (MQ=255)                                                                                                      aaccgtggttgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGAAGCag                                                                                                                                                            >  4:258945‑M2/196‑301 (MQ=255)                                                                                                          gtggttgtgccgattctacaatttgtctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAggt                                                                                                                                                        <  4:650059‑M2/110‑1 (MQ=255)                                                                                                               tgtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTAc                                                                                                                                                   <  1:158897‑M2/115‑1 (MQ=255)                                                                                                                gtgccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGGCGGAAACACTAGCGCATCATGGCGGGGaatagcaggccgcatttcagatcggaagagcgtcg                                                                                                                                                  >  2:710955‑M2/186‑284 (MQ=255)                                                                                                                 gtccttttttacaatttttctgctgcctttcgagcatgtcgaggctgccaggcgtttgttagctgtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaactattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAg                                                                                                                                                 <  2:258887‑M2/117‑1 (MQ=255)                                                                                                                   ccgattctacaatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTc                                                                                                                                               <  4:284147‑M2/119‑1 (MQ=255)                                                                                                                             aatttttctgctgcctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGggtag                                                                                                                                     <  3:337435‑M2/129‑1 (MQ=255)                                                                                                                                 ccctacacgacgctcttccgatctggggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTaa                                                                                                                                 <  2:269969‑M2/133‑1 (MQ=255)                                                                                                                                   ctacacgacgctcttccgatctgaggcttccaggcgtttttgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAAtt                                                                                                                               <  1:219168‑M2/135‑1 (MQ=255)                                                                                                                                          cctttcgagcatggcgaggctgccaggcgtttgtgggcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGTCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGAt                                                                                                                        <  2:547418‑M2/142‑1 (MQ=255)                                                                                                                                           ctttcgagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGAtt                                                                                                                       >  4:330233‑M2/159‑301 (MQ=255)                                                                                                                                                gagcatggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATTGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTgtt                                                                                                                  >  3:307735‑M2/154‑300 (MQ=255)                                                                                                                                                   catggcgaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTcatc                                                                                                               >  3:189527‑M2/151‑301 (MQ=255)                                                                                                                                                         gaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAggtggtacagtcaacagggtaggtaaagatcggaagagcacacgtctgaa                                                                                                         >  4:270035‑M2/145‑277 (MQ=255)                                                                                                                                                         gaggctgccaggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGAAGAGGTGgtacagtcaacagggtaggtaattagatcggaagagcacacgtctg                                                                                                         >  3:219166‑M2/145‑279 (MQ=255)                                                                                                                                                               gccgggcgttggtgagcggtattttgcccgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCt                                                                                                   <  2:330154‑M2/163‑1 (MQ=255)                                                                                                                                                                  aggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAAt                                                                                                >  3:629665‑M2/136‑301 (MQ=255)                                                                                                                                                                   ggcgtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATg                                                                                               >  4:564502‑M2/135‑301 (MQ=255)                                                                                                                                                                      gtttgtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGctc                                                                                            <  4:348037‑M2/170‑1 (MQ=255)                                                                                                                                                                      gtttgtgagcggtattttgcgcgcatgaatgcggcttttcccgttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGctc                                                                                            <  1:189529‑M2/170‑1 (MQ=255)                                                                                                                                                                          gtgagcggtattttgcgcgcatgaatgcggcttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTgg                                                                                        >  3:184342‑M2/128‑301 (MQ=255)                                                                                                                                                                                                        tttttccctttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGTTGTGACTAGAgg                                                          <  2:564332‑M2/204‑1 (MQ=255)                                                                                                                                                                                                                 ttcttgcgcatgggcattgcctataatgtcgccgacgtggaacctgacaacgatcgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGTTGTGACTAGAGGATCTGAatt                                                 >  1:24581‑M2/89‑300 (MQ=255)                                                                                                                                                                                                                          ctctttccctacacgacgctcttccgatctgaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGTTGTGACTAGAGGATCTGAATGTACAATGCa                                        <  2:120956‑M2/222‑1 (MQ=255)                                                                                                                                                                                                                                   gcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGTCTCAGTCCCAGTTGTGACTCGAGGCTCTGCATGTCCAATGCCcttcttctc                               >  2:325508‑M2/71‑297 (MQ=255)                                                                                                                                                                                                                                   gcctataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGTTGTGACTAGAGGATCTGAATGTACAATGCACTTCTGCTa                               >  2:584767‑M2/71‑301 (MQ=255)                                                                                                                                                                                                                                       ataatgtcgccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGTTGTGACTAGAGGATCTGAATGTACAATGCACTTCTGCTAATCg                           <  1:629663‑M2/235‑1 (MQ=255)                                                                                                                                                                                                                                               gccgacgtggaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGTTGTGACTAGAGGATCTGAATGTACAATGCACTTCTGCTAATCGCTGCAAtt                   <  4:584952‑M2/243‑1 (MQ=255)                                                                                                                                                                                                                                                        gaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGttgtgactagaggatctgaatgtacaatgcaagatcggaagagcacacgtctgaactccag          >  4:120981‑M2/50‑271 (MQ=255)                                                                                                                                                                                                                                                        gaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGTTGTGACTAGAGGATCTGAATGTACAATGCACTTCTGCTAATCGCTGCAATTGCCCTTAcc          >  3:141241‑M2/50‑301 (MQ=255)                                                                                                                                                                                                                                                        aaacctgacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGTTGTGACTAGAGGATCTGAATGTACAATGCACTTCTGCTAATCGCTGCAATTGCCCTTAcc          <  1:184344‑M2/252‑1 (MQ=255)                                                                                                                                                                                                                                                              gacaacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGTTGTGACTAGAGGATCTGAATGTACAATGCACTTCTGCTAATCGCTGCAATTGCCCTTACCGGATgc    <  1:141239‑M2/258‑1 (MQ=255)                                                                                                                                                                                                                                                                caacgattgataggttctggctttacccgtatcgcggccgcCCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGTTGTGACTAGAGGATCTGAATGTACAATGCACTTCTGCTAATCGCTGCAATTGCCCTTACCGGCTgcgc  >  4:533144‑M2/42‑301 (MQ=255)                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                    ACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACATTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAACTGGGCAAGAGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGATCAATTACTTTCGTGGTAAACACCGTGGCTGACCCTGCTTGATACTGGCCCGGCGGGAGGAGCGTCCGAGCGGAGCTGGTGGTCGCCGGAAACACTAGCGCATCATGGCGGTGAATAGCAGGCCGCATTTCGACGCAGAGGTGGTACAGTCAACAGGGTAGGTAATTCGCAGATTTTGTCATCCATCCAAGTGCTAATGCTCTTGGATGGGGCTCAGTCCCAGTTGTGACTAGAGGATCTGAATGTACAATGCACTTCTGCTAATCGCTGCAATTGCCCTTACCGGATGCGC  >  AE015451/2033889‑2034444

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