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 &