【Linux blast比对】fastq转fasta比对到参考序列

FASTQ转换为FASTA

nohup less -S ../sample-1_L3_2.fq.gz |awk '{if(NR%4 == 1){print ">" substr($0, 2)}}{if(NR%4 == 2){print}}' > ./sample-1_L3_2.fasta &

比对

$vi run_blast1.sh 

blastn \
 -query sample-1_L3_2.fasta \
 -db artificial.blastdb \
 -out sample_results2.xls \
 -outfmt 6 \
 -num_threads 8 

$nohup sh run_blast1.sh &

参考资料：

©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成，浏览时请结合常识与多方信息审慎甄别。
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容