问题描述：遇到一位同学他的SEQ数据是13年的，illumina双端测序，而且是无参物种，分别是病毒和昆虫拼接好的转录组序列。而且还不是同一批样本。然而我比较熟悉RNAseq这个DENOVo，从来没有做过，不过同学比较着急。他需要两个fasta中的交叉序列。因此我觉得用BLAST应该可以达到他的要求。以下为解决步骤方便以后再遇到类似情况：
1.安装blast+
sudo apt install blast
2.用参考序列建库
makeblastdb -in *.fasta -dbtype nucl -parse_seqids -out outdatabase
3.将需要比对的序列与上一步的核酸库进行比对，并输出结果
blastn -query /media/bai/final.need2.fa -out LAST.fasta.blast -db outdatabase -outfmt
2021-03-04 blast+ 比对两个denovo RNAseq序列
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 一：安装blast 1. 下载blast: $wget ftp://ftp.ncbi.nlm.nih.gov/bl...
- 1. Blast （1）格式化数据库formatdb -i db.seq -p T -o T -l logfile...
