flex 使用多个序列比对工具的共识序列管道

要求：用法： find_mapped_reads.py返回一个仅包含映射指定区域内读取的fasta文件。 python3 find_mapped_reads.py [fasta file with reads] [m5 mapping output from blasr] [start index] [end index] [output file name]示例： python3 find_mapped_reads.py reads.fa mapping.m5 100 5100 reads_100_5100.fa flex.py从多序列比对输出创建一致序列。 python3 flex.py [fasta file output from Muscle] [minimum number of non-gap nucleotides] [nucleotide pe