max_hsps_1_dict[bn.qseqid][bn.sseqid] = bnelse: max_hsps_1_dict[bn.qseqid] = {bn.sseqid: bn}# 第一步:检查normal.txt去重完成后是否只剩244行了normal_cnt =0forqseqidinnormal_dict: normal_cnt +=len(normal_dict[qseqid])print('去除query ID和subject ID完全相同的行后,%s还剩下%s...
Set maximum number of HSPs per subject sequence to save for each query 意思是对于每个subject,保留多少HSP。 2.2 测试 运行命令: $ blastn -query otu.fa -outfmt '6 qseqid sseqid qstart qend sstart send evalue' -db indexes/ref -out ret/max-hsps-1.txt -max_hsps 1 max-hsps-1.txt共有...