分别说明blastp –query query.fasta –db all –out blast.out –evalue 1e-10 –outfmt 7
每个参数的含义。
grep -v “#” blast.out > blast2.out命令中，-v和“>”代表什么,
说明awk ‘$3>50 && $11< 1e-30’ blast2.out >blast3.out 的含义
说明awk ‘{print $2}’ blast3.out的意义
说明grep -v “#” blast.out | awk ‘$3>50 && $11<1e-30 {print $2}’ | sort | uniq > protein_ids.txt 含义