通常做blast,我都是输出asn格式,然后用blast_formatter抽出需要的信息写成表格格式。

序列的description行当然不会有重复,但是ID (第一个空白前面的字符)是有重复的,从qseqid来看,就分不清是那个序列,于是准备让blast_formatter把query title(整个description line)给输出来,结果仔细阅读了blastn -help发现不支持呀不支持!然而它支持stitle和salltitle,也就是subject的可以,但query的不行,这绝逼是一个大坑!_| ̄|○

Continue reading

Author's picture

Guangchuang Yu

a senior-in-age-but-not-senior-in-knowledge bioinformatician

Postdoc researcher

Hong Kong