核酸-核酸
makeblastdb -in Sitalica.con -dbtype nucl -out Sitalicagenome
blastn -task blastn -query CCA1.txt -db Sitalica_312_v2.fa -outfmt 7 -evalue 0.00001 -max_target_seqs 1 -num_threads 20 -out CCA1_gffclean.txt
蛋白-蛋白
makeblastdb -in Nitab-v4.5_proteins_Edwards2017.fasta -parse_seqids -dbtype prot -out Nitab-v4.5_db
blastp -query cold_AT_gffclean.fa -db ../Nitab-v4.5_db -outfmt 7 -evalue 0.00001 -max_target_seqs 1 -num_threads 20 -out cold_AT_gffclean.txt
#Fields: query acc.ver, subject acc.ver, % identity, alignment length, mismatches, gap opens, q. start, q. end, s. start, s. end, evalue, bit score