mkdir /path/to/PBLAT/test cp hg19.fa.gz /path/to/PBLAT/test cp NONCODEv5_human.fa.gz /path/to/PBLAT/test
cd /path/to/PBLAT/test
gzip -d hg19.fa.gz gzip -d NONCODEv5_human.fa.gz cat NONCODEv5_human.fa | grep -v "^#" | head -799996 > test-799996.fa
{ time pblat -threads=96 hg19.fa test-799996.fa out.psl ; } 2>&1 |tee pblat-kunpeng.log
需要查看“pblat-kunpeng.log”日志中的“real”数值,单位是“s”,数值越小性能越优。
输出的结果样例如图1所示。