查看原文
其他

blast一大坑:无法输出query sequence title

2016-11-17 Y叔 biobabble

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


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



您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存