export PATH=/path/to/GMAP/Install/bin:$PATH
cp /path/to/PASA/PASApipeline-pasa-v2.4.1/pasa_conf/pasa.CONFIG.template /path/to/PASA/PASApipeline-pasa-v2.4.1/pasa_conf/conf.txt vim /path/to/PASA/PASApipeline-pasa-v2.4.1/pasa_conf/conf.txt
MYSQLSERVER=localhost MYSQL_RO_USER=test MYSQL_RO_PASSWORD=123456 MYSQL_RW_USER=pasa MYSQL_RW_PASSWORD=123456
cd /path/to/PASA/PASApipeline-pasa-v2.4.1/sample_data
cp /path/to/PASA/PASApipeline-pasa-v2.4.1/pasa_conf/pasa.alignAssembly.Template.txt alignAssembly.config vi alignAssembly.config DATABASE=sample_data
修改完成后输入:wq!保存退出。
gzip -d genome_sample.fasta.gz { time Launch_PASA_pipeline.pl -c alignAssembly.config -C -R -g genome_sample.fasta -t all_transcripts.fasta.clean -T -u all_transcripts.fasta -f FL_accs.txt --ALIGNERS blat,gmap --CPU 96; } 2>&1 |tee PASA.log
需要查看PASA.log日志,日志中打印Finished表示PASA正常结束。
输出的结果样例如图1所示。