1.使用命令随机从fastq文件中提取序列:
python /share/work/huangls/piplines/01.script/random_fastq_to_fa.py -f C5_R1.fastq.gz -o ~ -p its5 -n 100
2.将序列提交到NCBI中用blast做NR/NT库中进行比对:https://blast.ncbi.nlm.nih.gov/Blast.cgi
生物信息
702 篇文章
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!