使用Launch_PASA_pipeline.pl 脚本报错

使用Launch_PASA_pipeline.pl 脚本报错

报错信息:

(base) admin1@admin1-PowerEdge-R730xd:/mnt/6aa0485e-9dc6-4840-b61a-0798134475bf/jiegou$  /usr/local/src/PASApipeline/Launch_PASA_pipeline.pl  -c alignAssembly.config -C -R -t Trinity-GG.fasta --ALIGNERS minimap2 -g ./genome.fasta --CPU 68 1>pasa.log
-connecting to SQLite db: /tmp/Bra-pasa.sqlite
-*** Running PASA pipeine:
-- Skipping CMD: /usr/local/src/PASApipeline/scripts/create_sqlite_cdnaassembly_db.dbi -c alignAssembly.config -S '/usr/local/src/PASApipeline/schema/cdna_alignment_sqliteschema', checkpoint [/mnt/6aa0485e-9dc6-4840-b61a-0798134475bf/jiegou/__pasa_Bra-pasa.sqlite_SQLite_chkpts/create_db.ok] exists.
-- Skipping CMD: /usr/local/src/PASApipeline/scripts/upload_transcript_data.dbi -M '/tmp/Bra-pasa.sqlite' -t Trinity-GG.fasta  -f NULL , checkpoint [/mnt/6aa0485e-9dc6-4840-b61a-0798134475bf/jiegou/__pasa_Bra-pasa.sqlite_SQLite_chkpts/upload_transcripts.ok] exists.
* [Thu Nov 13 21:21:22 2025] Running CMD: /usr/local/src/PASApipeline/scripts/run_spliced_aligners.pl --aligners minimap2 --genome ./genome.fasta --transcripts Trinity-GG.fasta -I 500000 -N 1 --CPU 68
-running minimap2
CMD: /usr/local/src/PASApipeline/scripts/process_minimap2_alignments.pl --genome ./genome.fasta --transcripts Trinity-GG.fasta --CPU 68 -I 500000 -o Trinity-GG.fasta.mm2.bam
-- Skipping CMD: minimap2 -d ./genome.fasta.mm2/genome.fasta.mmi -t 68 ./genome.fasta, checkpoint [/mnt/6aa0485e-9dc6-4840-b61a-0798134475bf/jiegou/./genome.fasta.mm2/genome.fasta.mmi.ok] exists.
CMD: bash -c 'set -o pipefail && minimap2 -ax splice  --secondary=no -O6,24 -B4 -L -t 68 -cs -ub -G 500000   ./genome.fasta.mm2/genome.fasta.mmi Trinity-GG.fasta | samtools view -Sbt genome.fasta.fai | samtools sort -o Trinity-GG.fasta.mm2.bam && samtools index Trinity-GG.fasta.mm2.bam '
[WARNING] Indexing parameters (-k, -w or -H) overridden by parameters used in the prebuilt index.
[M::main::6.314*1.00] loaded/built the index for 152 target sequence(s)
[M::mm_mapopt_update::7.619*1.00] mid_occ = 834
[M::mm_idx_stat] kmer size: 15; skip: 10; is_hpc: 0; #seq: 152
[M::mm_idx_stat::8.519*1.00] distinct minimizers: 73053642 (43.22% are singletons); average occurrences: 4.613; average spacing: 5.341; total length: 1799778430
[M::worker_pipeline::63.164*51.25] mapped 227394 sequences
[M::main] Version: 2.22-r1101
[M::main] CMD: minimap2 -ax splice --secondary=no -O6,24 -B4 -L -t 68 -cs -ub -G 500000 ./genome.fasta.mm2/genome.fasta.mmi Trinity-GG.fasta
[M::main] Real time: 63.522 sec; CPU: 3237.524 sec; Peak RSS: 48.108 GB
[E::hts_idx_check_range] Region 536864449..536878569 cannot be stored in a bai index. Try using a csi index with min_shift = 14, n_lvls >= 6
[E::sam_index] Read 'TRINITY_GG_7592_c1_g1_i1' with ref_name='chr1', ref_length=559816472, flags=16, pos=536864450 cannot be indexed
samtools index: failed to create index for "Trinity-GG.fasta.mm2.bam": Numerical result out of range
Error, cmd: bash -c 'set -o pipefail && minimap2 -ax splice  --secondary=no -O6,24 -B4 -L -t 68 -cs -ub -G 500000   ./genome.fasta.mm2/genome.fasta.mmi Trinity-GG.fasta | samtools view -Sbt genome.fasta.fai | samtools sort -o Trinity-GG.fasta.mm2.bam && samtools index Trinity-GG.fasta.mm2.bam ' died with ret (256) at /usr/local/src/PASApipeline/scripts/process_minimap2_alignments.pl line 130.
Error, cmd:
/usr/local/src/PASApipeline/scripts/process_minimap2_alignments.pl --genome ./genome.fasta --transcripts Trinity-GG.fasta --CPU 68 -I 500000 -o Trinity-GG.fasta.mm2.bam
 died with ret (512) at /usr/local/src/PASApipeline/scripts/../PerlLib/Process_cmd.pm line 18.
    Process_cmd::process_cmd("/usr/local/src/PASApipeline/scripts/process_minimap2_alignmen"...) called at /usr/local/src/PASApipeline/scripts/run_spliced_aligners.pl line 149
    main::run_minimap2() called at /usr/local/src/PASApipeline/scripts/run_spliced_aligners.pl line 88
Error, cmd: /usr/local/src/PASApipeline/scripts/run_spliced_aligners.pl --aligners minimap2 --genome ./genome.fasta --transcripts Trinity-GG.fasta -I 500000 -N 1 --CPU 68 died with ret 512 No such file or directory at /usr/local/src/PASApipeline/PerlLib/Pipeliner.pm line 187.
    Pipeliner::run(Pipeliner=HASH(0x629238bb3210)) called at /usr/local/src/PASApipeline/Launch_PASA_pipeline.pl line 1061


请先 登录 后评论