evotools / nf-lo Goto Github PK
View Code? Open in Web Editor NEWA Nextflow workflow to generate lift over files for any pair of genomes
Home Page: https://nf-lo.readthedocs.io/
License: MIT License
A Nextflow workflow to generate lift over files for any pair of genomes
Home Page: https://nf-lo.readthedocs.io/
License: MIT License
Dear all,
I can't build the nflo singularity/docker image using the latest continuumio/miniconda3
base image. So I've pinned the base image and added the missing requirement. I will open a PR
Hi I am trying to use nf-LO pipeline and I am running into some errors. On my server I had specified the nflo pipeline to use 8 cpus and 32gb of RAM but it fails with the error " Process requirement exceed available CPUs -- req: 2; avail: 1" I was wondering why is that when I am giving it 8 cpus in both my lsf command and in the nextflow command.
nextflow run evotools/nf-LO
--source /scratch1/fs1/allegra.petti/khan.saad/Glass_synapse/reference/human_g1k_v37.fasta
--target /scratch1/fs1/allegra.petti/khan.saad/liftover_chain/hg19_igenomes/genome.fa
--aligner lastz
--tgtSize 10000000
--tgtOvlp 100000
--srcSize 20000000
--liftover_algorithm crossmap
--outdir ./my_liftover
--publish_dir_mode copy
--max_cpus 8
--max_memory 120.GB \
JAVA_HOME='/venv' export PATH=/opt/conda/bin:$PATH PATH=/venv/bin:$PATH LSF_DOCKER_VOLUMES='/home/khan.saad/:/home/khan.saad/ /storage1/fs1/allegra.petti/Active/:/storage1/fs1/al legra.petti/Active/ /scratch1/fs1/allegra.petti/:/scratch1/fs1/allegra.petti/' bsub -oo logs/nexflow_nflo.%J -G compute-allegra .petti -g /khan.saad/R_seurat -q general -M 128000000 -n 8 -R 'rusage[mem=128000]' -a 'docker(smk5g5/nf-lo:1.0.0)' /bin/bash ./ nflo.sh
Below is the error I am getting.
Command output:
(empty)
Work dir:
/scratch1/fs1/allegra.petti/khan.saad/liftover_chain/work/ae/5dbc7b0dff31f4dbfda486663a4019
Tip: you can replicate the issue by changing to the process work dir and entering the command bash .command.run
WARN: Killing pending tasks (1)
executor > local (7044)
[99/b8dfba] process > PREPROC:src2bit (src2Bit) [100%] 1 of 1 ✔
[77/b6fe1e] process > PREPROC:tgt2bit (tgt2Bit) [100%] 1 of 1 ✔
[87/500b9a] process > PREPROC:splitsrc (splitsrc) [100%] 1 of 1 ✔
[5d/0611b0] process > PREPROC:groupsrc (groupsrc) [100%] 1 of 1 ✔
[ca/a47e09] process > PREPROC:splittgt (splittgt) [100%] 1 of 1 ✔
[c4/622932] process > PREPROC:grouptgt (grouptgt) [100%] 1 of 1 ✔
[03/40d546] process > PREPROC:pairs (mkpairs) [100%] 1 of 1 ✔
[e1/c85215] process > ALIGNER:lastz (lastz_med.sr... [ 15%] 7009 of 46610, fa...
[8f/e0eb1a] process > ALIGNER:axtChain (axtchain_m) [ 0%] 28 of 7007
[- ] process > ALIGNER:chainMerge -
[- ] process > ALIGNER:chainNet -
[- ] process > ALIGNER:netSynt -
[- ] process > ALIGNER:chainsubset -
[- ] process > ALIGNER:chain2maf -
[- ] process > ALIGNER:name_maf_seq -
[- ] process > ALIGNER:mafstats -
[d1/3b95a6] NOTE: Process ALIGNER:lastz (lastz_med.src93.tgt68)
terminated with an error exit status (1) -- Execution is retried (1)
Error executing process > 'ALIGNER:lastz (lastz_med.src93.tgt68)'
Caused by:
Process requirement exceed available CPUs -- req: 2; avail: 1
Command executed:
echo B=0 C=0 E=30 H=0 K=3000 L=3000 M=50 O=400 T=1 Y=9400
lastz /scratch1/fs1/allegra.petti/khan.saad/liftover_chain/work/5d/0611b088a66d4b63a1adb6961c96d7/CLUST_src/src93.fa /scratch1/fs1/allegra.petti/khan.saad/liftover_chain/work/c4/622932c827e301415a2268ee960622/CLUST_tgt/tgt68.fa B=0 C=0 E=30 H=0 K=3000 L=3000 M=50 O=400 T=1 Y=9400 --ambiguous=iupac --format=lav |
lavToPsl stdin stdout |
liftUp -type=.psl stdout source.lift warn stdin |
liftUp -type=.psl -pslQ src93.tgt68.psl target.lift warn stdin
Command output:
(empty)
Work dir:
/scratch1/fs1/allegra.petti/khan.saad/liftover_chain/work/ae/5dbc7b0dff31f4dbfda486663a4019
Hello,
I am having similar issues with creating environment.yml during the workflow run. I am running it using docker using the following command,
nextflow run evotools/nf-LO -with-docker nflo:latest --source source.fa --target target.fa --distance near --aligner lastz --outdir Output_nfLo --profile conda
But, it runs all the way to make_report, and then stalls while trying to create the conda environment, and complains there are lot of versions of the packages in conflict. I thought running it via docker would not cause any issue with dependency.
Any idea what I may be doing wrong ?
Dear nd-LO developers,
I'm trying to align two local genome assemblies using nfLO, and get the following error at the bottom of the 'reports/execution_trace.txt'
file. Could you please help me figure out why?
701855 fc/df1139 144312 ALIGNER:axtChain (axtchain_m) COMPLETED 0 2023-02-06 05:52:10.789 926ms 733ms 101.9% 133.3 MB 415.4 MB 90 MB 3.6 KB
701856 ca/843384 160553 ALIGNER:chainMerge (chainmerge) FAILED 139 2023-02-06 06:22:34.216 21.1s 21s - - - - -
701857 dd/25d4a1 160837 ALIGNER:chainMerge (chainmerge) FAILED 139 2023-02-06 06:23:02.349 21.1s 21.1s - - - - -
701858 76/995157 161063 ALIGNER:chainMerge (chainmerge) FAILED 139 2023-02-06 06:23:30.339 21.5s 21.5s - - - - -
```````
The code I used was:
nextflow run evotools/nf-LO --source ../GCA_021130815.1_PanTigT.MC.v3_genomic.fna.gz --target ~/data/hg38.fa.masked.gz --outdir ../NFLOW
When the code failed, it had created three directories: genome2bit, singlechains, and reports.
Please let me know if you need more information from me,
Shweta
I got the error information when runing the n-Lo,this log showing in blow:
WARNING: Could not load nf-core/config profiles: https://raw.githubusercontent.com/nf-core/configs/master/nfcore_custom.config
Launching https://github.com/evotools/nf-LO
[pedantic_ptolemy] DSL2 - revision: 312fb56 [main]
Nextflow LiftOver v 1.6.0
======================================
source : chrY.oldVersin.fa
target : DRC.chrY.fa
aligner : minimap2
distance : medium
custom align : -cx asm5 -l 10000
custom chain : false
source chunk : 20000000
source overlap : 0
target chunk : 10000000
target overlap : 100000
output folder : ./liftover_Minimap2
liftover name : liftover
annot : false
annot type : false
liftover meth. : liftover
igenomes_base : s3://ngi-igenomes/igenomes/
igenomes_ignore : false
no_maf : false
skip netsynt : true
max cpu : 143
max mem : 26 GB
max rt : 240.h
Using CrossMap
executor > local (5)
[13/61a6e5] process > PREPROC:src2bit (src2Bit) [100%] 1 of 1 ✔
executor > local (10)
[13/61a6e5] process > PREPROC:src2bit (src2Bit) [100%] 1 of 1 ✔
[cb/e667b4] process > PREPROC:tgt2bit (tgt2Bit) [100%] 1 of 1 ✔
[3e/c206c5] process > PREPROC:splitsrc (splitsrc) [100%] 1 of 1 ✔
[83/de485d] process > PREPROC:splittgt (splittgt) [100%] 1 of 1 ✔
[73/c5a594] process > PREPROC:pairs (mkpairs) [100%] 1 of 1 ✔
[58/cefc2c] process > ALIGNER:minimap2 (minimap2.medium.src.tgt) [100%] 1 of 1 ✔
[26/7d4da4] process > ALIGNER:axtChain (axtchain_m) [100%] 1 of 1 ✔
[33/4e27ca] process > ALIGNER:chainMerge (chainmerge) [100%] 3 of 3, failed: 3, retries: 2 ✔
[- ] process > ALIGNER:chainNet -
[- ] process > ALIGNER:chainsubset -
[- ] process > ALIGNER:chain2maf -
[- ] process > ALIGNER:name_maf_seq -
[- ] process > ALIGNER:mafstats -
[c7/cb6d2d] NOTE: Process ALIGNER:chainMerge (chainmerge)
terminated with an error exit status (1) -- Execution is retried (1)
[8d/9085f1] NOTE: Process ALIGNER:chainMerge (chainmerge)
terminated with an error exit status (1) -- Execution is retried (2)
[33/4e27ca] NOTE: Process ALIGNER:chainMerge (chainmerge)
terminated with an error exit status (1) -- Error is ignored
I split the genome into 20 chromosomes and then ran nf-Lo with blat alignments for each chromosome. I don't know how to merge all chain files into one chain, can you help me ?
I've tried install nf-lo manually with this command:
wget https://raw.githubusercontent.com/evotools/nf-LO/main/install.sh
but it returned a 404 error. the log is displayed below:
--2023-06-07 16:14:11-- https://raw.githubusercontent.com/evotools/nf-LO/main/install.sh
Resolving raw.githubusercontent.com (raw.githubusercontent.com)... 185.199.109.133, 185.199.111.133, 185.199.108.133
Connecting to raw.githubusercontent.com (raw.githubusercontent.com)|185.199.109.133|:443... connected.
HTTP request sent, awaiting response... 404 Not Found
2023-06-07 16:14:11 ERROR 404: Not Found.
I'm looking forward to a reply. thanks.
Hi,
I met a problem about using nf-LO to liftover my genomes.
About Annotation param, If i support a gff file, this file is source genome's gff or target genome's gff?
Can you help me out?
Much appreciate.
Hi,
When I using nf-LO to get chain file, I met a problem.
Below is my code.
get_chain_file(){
source ~/miniconda3/bin/activate nf-LO
nextflow run evotools/nf-LO --source $1 \
--target $2 \
--annotation $3 \
--annotation_format gff \
--distance near \
--aligner $4 \
--liftover_algorithm liftover \
--outdir $5 \
--max_cpus $6 \
--max_memory 8.GB \
-profile conda
}
get_chain_file $dhl92_ref $ref_13C $dhl92_gff GSAlign DHL92_Lift_To_13C 48
Below is error log:
executor > local (18)
[0a/0b9198] process > PREPROC:src2bit (src2Bit) [100%] 1 of 1 ✔
[ca/375eee] process > PREPROC:tgt2bit (tgt2Bit) [100%] 1 of 1 ✔
[d6/67237f] process > PREPROC:splitsrc (splitsrc) [100%] 1 of 1 ✔
[06/67ee7a] process > PREPROC:splittgt (splittgt) [100%] 1 of 1 ✔
[3c/cd32ff] process > PREPROC:pairs (mkpairs) [100%] 1 of 1 ✔
[2f/9671a6] process > ALIGNER:bwt_index (bwt_index) [100%] 1 of 1 ✔
[fc/c8c3bd] process > ALIGNER:gsalign (gsalign_ne... [100%] 1 of 1 ✔
[46/9b560b] process > ALIGNER:axtChain (axtchain_n) [100%] 1 of 1 ✔
[1d/f4c0bd] process > ALIGNER:chainMerge (chainme... [100%] 1 of 1 ✔
[2e/e1b000] process > ALIGNER:chainNet (chainnet) [100%] 1 of 1 ✔
[0c/2ea9e2] process > ALIGNER:netSynt (netSyntenic) [100%] 1 of 1 ✔
[98/9f1882] process > ALIGNER:chainsubset (chains... [100%] 1 of 1 ✔
[48/1eb5ae] process > ALIGNER:chain2maf (chainmaf) [100%] 1 of 1 ✔
[4e/b212f0] process > ALIGNER:name_maf_seq (namemaf) [100%] 1 of 1 ✔
[eb/1065eb] process > ALIGNER:mafstats (mafstats) [100%] 1 of 1 ✔
[1e/428505] process > LIFTOVER:lifter (crossmap) [100%] 3 of 3, failed: 3...
[- ] process > LIFTOVER:features_stats -
[- ] process > make_report -
[5d/627e71] NOTE: Process LIFTOVER:lifter (crossmap)
terminated with an error exit status (127) -- Execution is retried (1)
[43/56ce02] NOTE: Process LIFTOVER:lifter (crossmap)
terminated with an error exit status (127) -- Execution is retried (2)
[1e/428505] NOTE: Process LIFTOVER:lifter (crossmap)
terminated with an error exit status (127) -- Error is ignored
executor > local (18)
[0a/0b9198] process > PREPROC:src2bit (src2Bit) [100%] 1 of 1 ✔
[ca/375eee] process > PREPROC:tgt2bit (tgt2Bit) [100%] 1 of 1 ✔
[d6/67237f] process > PREPROC:splitsrc (splitsrc) [100%] 1 of 1 ✔
[06/67ee7a] process > PREPROC:splittgt (splittgt) [100%] 1 of 1 ✔
[3c/cd32ff] process > PREPROC:pairs (mkpairs) [100%] 1 of 1 ✔
[2f/9671a6] process > ALIGNER:bwt_index (bwt_index) [100%] 1 of 1 ✔
[fc/c8c3bd] process > ALIGNER:gsalign (gsalign_ne... [100%] 1 of 1 ✔
[46/9b560b] process > ALIGNER:axtChain (axtchain_n) [100%] 1 of 1 ✔
[1d/f4c0bd] process > ALIGNER:chainMerge (chainme... [100%] 1 of 1 ✔
[2e/e1b000] process > ALIGNER:chainNet (chainnet) [100%] 1 of 1 ✔
[0c/2ea9e2] process > ALIGNER:netSynt (netSyntenic) [100%] 1 of 1 ✔
[98/9f1882] process > ALIGNER:chainsubset (chains... [100%] 1 of 1 ✔
[48/1eb5ae] process > ALIGNER:chain2maf (chainmaf) [100%] 1 of 1 ✔
[4e/b212f0] process > ALIGNER:name_maf_seq (namemaf) [100%] 1 of 1 ✔
[eb/1065eb] process > ALIGNER:mafstats (mafstats) [100%] 1 of 1 ✔
[1e/428505] process > LIFTOVER:lifter (crossmap) [100%] 3 of 3, failed: 3...
[- ] process > LIFTOVER:features_stats -
[- ] process > make_report -
Completed at: 20-Feb-2024 17:20:07
Duration : 12m 32s
CPU hours : 2.5 (0.1% failed)
Succeeded : 15
Ignored : 1
Failed : 3
I don't know what happend because I can't get a detail error log. How to solve this problem?
Much appreciate!
Hi,
I have a failure in running the pipline, couldnt understand exactly why.
Maybe it is looking for chain.null file?
Appreciate your help.
Thank you!
Here is the nextflow.log file (was to long to paste): https://tauex-my.sharepoint.com/:u:/g/personal/udiland_tauex_tau_ac_il/EVrh6dGhb4BOieTEvjfMbo0BcKncSEU42bVw8KarMUbEHg?e=yKs6Nu
Hello, thank you for developing this tool, I tried to run the following code :
nextflow run evotools/nf-LO --source /home/sealight1999/references/GRCH38/GRCh38_no_alt.fna --target /home/sealight1999/references/CHM13/ chm13v2.0_1.fa --aligner minimap2 –outdir /home/sealight1999/references -profile local
But encountered following warnings:
[c3/9fda46] NOTE: Process PREPROC:tgt2bit (tgt2Bit)
terminated with an error exit status (127) -- Execution is retried (1)
[ae/47e4dc] NOTE: Process PREPROC:splitsrc (splitsrc)
terminated with an error exit status (127) -- Execution is retried (1)
[e7/fdb913] NOTE: Process PREPROC:src2bit (src2Bit)
terminated with an error exit status (127) -- Execution is retried (1)
[fd/3eadf0] NOTE: Process PREPROC:splittgt (splittgt)
terminated with an error exit status (127) -- Execution is retried (1)
[03/acbfe7] NOTE: Process PREPROC:tgt2bit (tgt2Bit)
terminated with an error exit status (127) -- Execution is retried (2)
[8e/ecf45f] NOTE: Process PREPROC:splitsrc (splitsrc)
terminated with an error exit status (127) -- Execution is retried (2)
[22/e11aa5] NOTE: Process PREPROC:src2bit (src2Bit)
terminated with an error exit status (127) -- Execution is retried (2)
[62/97c512] NOTE: Process PREPROC:splittgt (splittgt)
terminated with an error exit status (127) -- Execution is retried (2)
[83/bd21fa] NOTE: Process PREPROC:tgt2bit (tgt2Bit)
terminated with an error exit status (127) -- Error is ignored
[3f/930365] NOTE: Process PREPROC:splitsrc (splitsrc)
terminated with an error exit status (127) -- Error is ignored
[9e/fdca06] NOTE: Process PREPROC:src2bit (src2Bit)
terminated with an error exit status (127) -- Error is ignored
[77/94ba10] NOTE: Process PREPROC:splittgt (splittgt)
terminated with an error exit status (127) -- Error is ignored
WARN: Graphviz is required to render the execution DAG in the given format -- See http://www.graphviz.org for more info.
What could be wrong? Thank you!
Hello.
Thank you for no-LO.
I am running it with —annotation and —annotation_format. I am using the conda profile.
It runs all the way to make_report, and then stalls while trying to create the conda environment.
I have tried to run conda env create directly to troubleshoot. Essentially, it gets stuck in the dependency solver section. By using the newer lmamba solver, it appears that some of the versions of the packages conflict. I am wondering if the versions could be relaxed or removed all together.
Thanks again.
PS: It is the dependency hell problem. 😣
Hi,
I'm running nf-LO between hg38 and pantro6 using the minimap2 preset. The docker container used was built as suggested in the docs. I'm running this on a SLURM cluster launched via a login node. It looks like minimap2 fails but the .exitcode
is still 0 so the downstream processes fail. I can confirm that the same files and settings work with lastz but not minimap2.
I've also noticed that only one minimap2 process spawns whereas a few hundred lastz processes typically spawn. It seems that splitsrc
and splittgt
each return a large single multifasta.
run script:
nextflow run evotools/nf-LO \
-latest \
-c phoenix.config \
--outdir pantro6_to_hg38_nflo \
--distance near \
--aligner minimap2 \
--max_cpus 32 \
--max_memory 128.GB \
-w work3 \
--source panTro6_cleaned.fa \
--target hg38_cleaned.fa \
-with-tower \
-profile docker -with-docker vpeddu/nflo:latest
process .command.error
:
Got 30 lifts in source.lift
Lifting stdin
Got 328 lifts in target.lift
Lifting stdin
[M::mm_idx_gen::42.331*1.45] collected minimizers
[M::mm_idx_gen::68.360*2.79] sorted minimizers
[M::main::68.360*2.79] loaded/built the index for 30 target sequence(s)
[M::mm_mapopt_update::174.858*1.28] mid_occ = 121
[M::mm_idx_stat] kmer size: 19; skip: 19; is_hpc: 0; #seq: 30
[M::mm_idx_stat::259.800*0.94] distinct minimizers: 213075345 (92.73% are singletons); average occurrences: 1.328; average spacing: 10.079; total length: 2851453592
stdin is empty
Hello @RenzoTale88.
My machine has 16 CPUs, and by default no-LO sets the max to 15 cpus. I think that is reasonable.
But, I noticed that both minimap2 and gsalign, which share the cpus definition in the config files, end up with only 7 cpus while running. Is there a particular reason for halving the cpus for these two tools?
I noticed that the function check_max returns the minimum of max cpus / 2 * the number of tries and the max number of cpus. So, it will always return at most half the number of max cpus. Is that the intended behaviour?
Thanks again.
Is there an option for nf-LO to simply generate the pipeline commands into a text file, so that I can run every step one-by-one.
nf-LO Keeps randomly crashing, is has been a highly frustrating experience.
Hi,
I am getting this general error and cant figure out why.
This is the output I get.
Appreciate your help
Thank you
I run this in the termimal in conda env nf-LO
./nextflow run evotools/nf-LO -profile test,conda
and returns the following info:
N E X T F L O W ~ version 22.10.0
Launching https://github.com/evotools/nf-LO
[sick_kimura] DSL2 - revision: 5583780 [main]
Nextflow LiftOver v 1.7.0
======================================
source : https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/001/549/955/GCF_001549955.1_ASM154995v1/GCF_001549955.1_ASM154995v1_genomic.fna.gz
target : https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/011/751/205/GCF_011751205.1_ASM1175120v1/GCF_011751205.1_ASM1175120v1_genomic.fna.gz
aligner : lastz
distance : medium
custom align : false
custom chain : false
source chunk : 1000000
source overlap : 0
target chunk : 2000000
target overlap : 100000
output folder : /media/ggj/WHY/package/nextflow/outputs
liftover name : liftover
annot : https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/001/549/955/GCF_001549955.1_ASM154995v1/GCF_001549955.1_ASM154995v1_genomic.gff.gz
annot type : gff
liftover meth. : liftover
igenomes_base : s3://ngi-igenomes/igenomes/
igenomes_ignore : false
mamba : false
no_maf : false
skip netsynt : false
max cpu : 48
max mem : 4.GB
max rt : 240.h
Using CrossMap
[- ] process > PREPROC:src2bit -
[- ] process > PREPROC:src2bit -
[- ] process > PREPROC:tgt2bit -
executor > local (1)
[1d/9a29d4] process > PREPROC:src2bit (src2Bit) [ 0%] 0 of 1
executor > local (8)
[f5/18c1ae] process > PREPROC:src2bit (src2Bit) [100%] 3 of 3, failed: 3, retries: 2 ✔
[b2/334e03] process > PREPROC:tgt2bit (tgt2Bit) [100%] 1 of 1, failed: 1, retries: 1
[c8/9af8a0] process > PREPROC:splitsrc (splitsrc) [100%] 3 of 3, failed: 3, retries: 2 ✔
[- ] process > PREPROC:groupsrc -
[5a/486c00] process > PREPROC:splittgt (splittgt) [ 0%] 0 of 1
[- ] process > PREPROC:grouptgt -
[- ] process > PREPROC:pairs -
[- ] process > ALIGNER:lastz -
executor > local (12)
[f5/18c1ae] process > PREPROC:src2bit (src2Bit) [100%] 3 of 3, failed: 3, retries: 2 ✔
[b0/1f643c] process > PREPROC:tgt2bit (tgt2Bit) [100%] 3 of 3, failed: 3, retries: 2 ✔
[c8/9af8a0] process > PREPROC:splitsrc (splitsrc) [100%] 3 of 3, failed: 3, retries: 2 ✔
[- ] process > PREPROC:groupsrc -
[b7/e1aa20] process > PREPROC:splittgt (splittgt) [100%] 3 of 3, failed: 3, retries: 2 ✔
[- ] process > PREPROC:grouptgt -
[- ] process > PREPROC:pairs -
[- ] process > ALIGNER:lastz -
[- ] process > ALIGNER:axtChain -
[- ] process > ALIGNER:chainMerge -
[- ] process > ALIGNER:chainNet -
[- ] process > ALIGNER:netSynt -
[- ] process > ALIGNER:chainsubset -
[- ] process > ALIGNER:chain2maf -
[- ] process > ALIGNER:name_maf_seq -
executor > local (12)
[f5/18c1ae] process > PREPROC:src2bit (src2Bit) [100%] 3 of 3, failed: 3, retries: 2 ✔
[b0/1f643c] process > PREPROC:tgt2bit (tgt2Bit) [100%] 3 of 3, failed: 3, retries: 2 ✔
[c8/9af8a0] process > PREPROC:splitsrc (splitsrc) [100%] 3 of 3, failed: 3, retries: 2 ✔
[- ] process > PREPROC:groupsrc -
[b7/e1aa20] process > PREPROC:splittgt (splittgt) [100%] 3 of 3, failed: 3, retries: 2 ✔
[- ] process > PREPROC:grouptgt -
[- ] process > PREPROC:pairs -
[- ] process > ALIGNER:lastz -
[- ] process > ALIGNER:axtChain -
[- ] process > ALIGNER:chainMerge -
[- ] process > ALIGNER:chainNet -
[- ] process > ALIGNER:netSynt -
[- ] process > ALIGNER:chainsubset -
[- ] process > ALIGNER:chain2maf -
[- ] process > ALIGNER:name_maf_seq -
[- ] process > ALIGNER:mafstats -
[- ] process > LIFTOVER:lifter -
[- ] process > LIFTOVER:features_stats -
[- ] process > make_report -
Staging foreign file: https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/001/549/955/GCF_001549955.1_ASM154995v1/GCF_001549955.1_ASM154995v1_genomic.fna.gz
[1d/9a29d4] NOTE: Process PREPROC:src2bit (src2Bit)
terminated with an error exit status (1) -- Execution is retried (1)
[d5/a639d1] NOTE: Process PREPROC:splitsrc (splitsrc)
terminated with an error exit status (1) -- Execution is retried (1)
[49/3aef03] NOTE: Process PREPROC:src2bit (src2Bit)
terminated with an error exit status (1) -- Execution is retried (2)
[81/d5e694] NOTE: Process PREPROC:splitsrc (splitsrc)
terminated with an error exit status (1) -- Execution is retried (2)
[f5/18c1ae] NOTE: Process PREPROC:src2bit (src2Bit)
terminated with an error exit status (1) -- Error is ignored
[c8/9af8a0] NOTE: Process PREPROC:splitsrc (splitsrc)
terminated with an error exit status (1) -- Error is ignored
[b2/334e03] NOTE: Process PREPROC:tgt2bit (tgt2Bit)
terminated with an error exit status (1) -- Execution is retried (1)
[5a/486c00] NOTE: Process PREPROC:splittgt (splittgt)
terminated with an error exit status (1) -- Execution is retried (1)
[d3/70a534] NOTE: Process PREPROC:tgt2bit (tgt2Bit)
terminated with an error exit status (1) -- Execution is retried (2)
[e0/caf685] NOTE: Process PREPROC:splittgt (splittgt)
terminated with an error exit status (1) -- Execution is retried (2)
[b0/1f643c] NOTE: Process PREPROC:tgt2bit (tgt2Bit)
terminated with an error exit status (1) -- Error is ignored
[b7/e1aa20] NOTE: Process PREPROC:splittgt (splittgt)
terminated with an error exit status (1) -- Error is ignored
looking forward to your reply!Thanks!
Hello! I am using nf-LO to bulid liftover chains for Facs6toFacs5, but I met the following errs. Does it mean the disk space is not enough? (I am not sure) Looking forward to your reply. Thanks!
BTW,I use conda environmet and the example data runs successfully(-profile test,conda).
Here is the err info:
executor > local (50)
[37/3b6a75] process > PREPROC:src2bit (src2Bit) [100%] 1 of 1 ✔
[08/9a11ef] process > PREPROC:tgt2bit (tgt2Bit) [100%] 1 of 1 ✔
[88/1c52c2] process > PREPROC:splitsrc (splitsrc) [100%] 1 of 1 ✔
[6f/def6e1] process > PREPROC:groupsrc (groupsrc) [100%] 1 of 1 ✔
[ba/1794c8] process > PREPROC:splittgt (splittgt) [100%] 1 of 1 ✔
[17/acab5e] process > PREPROC:grouptgt (grouptgt) [100%] 1 of 1 ✔
[87/a62aee] process > PREPROC:pairs (mkpairs) [100%] 1 of 1 ✔
[90/f5326a] process > ALIGNER:lastz (lastz_near.2... [ 3%] 10 of 265, failed...
[6a/8f2413] process > ALIGNER:axtChain (axtchain_n) [100%] 7 of 7
[- ] process > ALIGNER:chainMerge -
[- ] process > ALIGNER:chainNet -
[- ] process > ALIGNER:chainsubset -
[- ] process > ALIGNER:chain2maf -
[- ] process > ALIGNER:name_maf_seq -
[- ] process > ALIGNER:mafstats -
[db/556bfd] NOTE: Process ALIGNER:lastz (lastz_near.11.tgt8)
failed -- Execution is retried (1)
[ac/e063b1] NOTE: Process ALIGNER:lastz (lastz_near.6.tgt10)
failed -- Execution is retried (1)
[26/33681b] NOTE: Process ALIGNER:lastz (lastz_near.6.tgt1)
failed -- Execution is retried (1)
executor > local (54)
[37/3b6a75] process > PREPROC:src2bit (src2Bit) [100%] 1 of 1 ✔
[08/9a11ef] process > PREPROC:tgt2bit (tgt2Bit) [100%] 1 of 1 ✔
[88/1c52c2] process > PREPROC:splitsrc (splitsrc) [100%] 1 of 1 ✔
[6f/def6e1] process > PREPROC:groupsrc (groupsrc) [100%] 1 of 1 ✔
[ba/1794c8] process > PREPROC:splittgt (splittgt) [100%] 1 of 1 ✔
[17/acab5e] process > PREPROC:grouptgt (grouptgt) [100%] 1 of 1 ✔
[87/a62aee] process > PREPROC:pairs (mkpairs) [100%] 1 of 1 ✔
[a6/1e196a] process > ALIGNER:lastz (lastz_near.9... [ 5%] 14 of 270, failed...
[6a/8f2413] process > ALIGNER:axtChain (axtchain_n) [100%] 7 of 7
[- ] process > ALIGNER:chainMerge -
[- ] process > ALIGNER:chainNet -
[- ] process > ALIGNER:chainsubset -
[- ] process > ALIGNER:chain2maf -
[- ] process > ALIGNER:name_maf_seq -
[- ] process > ALIGNER:mafstats -
[db/556bfd] NOTE: Process ALIGNER:lastz (lastz_near.11.tgt8)
failed -- Execution is retried (1)
[ac/e063b1] NOTE: Process ALIGNER:lastz (lastz_near.6.tgt10)
failed -- Execution is retried (1)
[26/33681b] NOTE: Process ALIGNER:lastz (lastz_near.6.tgt1)
failed -- Execution is retried (1)
[13/b53784] NOTE: Process ALIGNER:lastz (lastz_near.7.tgt21)
failed -- Execution is retried (1)
[4e/94bcd1] NOTE: Process ALIGNER:lastz (lastz_near.11.tgt22)
failed -- Execution is retried (1)
[37/566b5a] NOTE: Process ALIGNER:lastz (lastz_near.11.tgt4)
failed -- Execution is retried (1)
[b9/3c93d9] NOTE: Process ALIGNER:lastz (lastz_near.11.tgt1)
failed -- Execution is retried (1)
Hi,
I'm running a job on to mammalian genomes.
nextflow run ../bosTauTomBal/nf-LO --source ../bosTau6.fa --target /project/mchaisso_100/projects/Whales/mPhosSin1.pri/assembly.orig.fasta --max-cpus 16 --max_memory 64.GB -with-trace -with-report -resume -profile slurm,conda
At the step:
[54/e0d255] process > ALIGNER:chainMerge (chainmerge) [100%] 3 of 3, failed: 3, retries: 2
[- ] process > ALIGNER:chainNet -
[- ] process > ALIGNER:netSynt -
[- ] process > ALIGNER:chainsubset -
[- ] process > ALIGNER:chain2maf -
[- ] process > ALIGNER:name_maf_seq -
[- ] process > ALIGNER:mafstats -
[9d/5d42ff] NOTE: Error submitting process 'ALIGNER:chainMerge (chainmerge)' for execution -- Execution is retried (1)
[55/e09c1b] NOTE: Error submitting process 'ALIGNER:chainMerge (chainmerge)' for execution -- Execution is retried (2)
[54/e0d255] NOTE: Error submitting process 'ALIGNER:chainMerge (chainmerge)' for execution -- Error is ignored
slurm is not able to submit the command:
sbatch $PWD/work/54/e0d255bdfcaa440fbfe84e1e09cee3/.command.run
sbatch: error: Batch job submission failed: Pathname of a file, directory or other parameter too long
The problem is the file is too large (5.3M):
-rw-rw---- 1 mchaisso mchaisso_100 5.3M Dec 29 19:09 work/54/e0d255bdfcaa440fbfe84e1e09cee3/.command.run
I think the default slurm max script is 4.5 M, and I don't have permissions to increase this. The problem is in the zillions of rm and ln lines in nxf_stage(). Ideally nextflow could do the work of creating a script that runs the rm and ln's, and to submit that.
[0f/a27c7f] NOTE: Process PREPROC:grouptgt (grouptgt)
terminated with an error exit status (127) -- Execution is retried (1)
[6c/14fce7] NOTE: Process PREPROC:grouptgt (grouptgt)
terminated with an error exit status (127) -- Execution is retried (2)
[7e/a18dee] NOTE: Process PREPROC:grouptgt (grouptgt)
terminated with an error exit status (127) -- Error is ignored
WARN: Graphviz is required to render the execution DAG in the given format -- See http://www.graphviz.org for more info.
Completed at: 20-Nov-2023 20:40:11
A declarative, efficient, and flexible JavaScript library for building user interfaces.
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google ❤️ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.