Comments (10)
@815325223 happy to see your interest in RNACocktail.
I guess your volumes are not mounted correctly to the docker image. Can you try to mount to /work_dir/data and /work_dir/work instead.
from rnacocktail.
@msahraeian I don't think I have a mount problem. I tried using your mount path.
run_rnacocktail.py editing --alignment /work_dir/data/D8.sorted.bam --variant /work_dir/data/D8.snv.vcf --strand_pos /work_dir/data/hg19_strand_pos.bed --genes_pos /work_dir/data/hg19_genes_pos.bed --outdir /work_dir/work --workdir /work_dir/work --threads 10 --sample D8 --ref_genome /work_dir/data/Homo_sapiens.GRCh37.dna.primary_assembly.fa --knownsites /work_dir/data/common_all_20180418.vcf --giremi_dir /user/local/bin/
INFO 2019-03-04 01:47:46,022 src.utils Running RNASeqPipeline 0.2.2
INFO 2019-03-04 01:47:46,023 src.utils Command-line /usr/local/bin/run_rnacocktail.py editing --alignment /work_dir/data/D8.sorted.bam --variant /work_dir/data/D8.snv.vcf --strand_pos /work_dir/data/hg19_strand_pos.bed --genes_pos /work_dir/data/hg19_genes_pos.bed --outdir /work_dir/work --workdir /work_dir/work --threads 10 --sample D8 --ref_genome /work_dir/data/Homo_sapiens.GRCh37.dna.primary_assembly.fa --knownsites /work_dir/data/common_all_20180418.vcf --giremi_dir /user/local/bin/
INFO 2019-03-04 01:47:46,023 src.utils Arguments are Namespace(VariantAnnotator_opts='', alignment='/work_dir/data/D8.sorted.bam', editing_caller='GIREMI', gatk='GenomeAnalysisTK.jar', genes_pos='/work_dir/data/hg19_genes_pos.bed', giremi_dir='/user/local/bin/', giremi_opts='', htslib_dir='', java='java', java_opts='-Xms1g -Xmx5g', knownsites='/work_dir/data/common_all_20180418.vcf', mode='editing', outdir='/work_dir/work', ref_genome='/work_dir/data/Homo_sapiens.GRCh37.dna.primary_assembly.fa', sample='D8', samtools='samtools', start=0, strand_pos='/work_dir/data/hg19_strand_pos.bed', threads=10, timeout=10000000, variant='/work_dir/data/D8.snv.vcf', workdir='/work_dir/work')
INFO 2019-03-04 01:47:46,023 src.utils Run log will be saved in /work_dir/work/logs/run-editing-20190304-014746.log
INFO 2019-03-04 01:47:46,023 src.utils Run in mode: editing
INFO 2019-03-04 01:47:46,023 src.utils Running RNA editing calling step using GIREMI
INFO 2019-03-04 01:47:46,023 src.run_editing Running RNA editing detection (GIREMI) for D8
ERROR 2019-03-04 01:47:46,024 src.run_editing Aborting!
INFO 2019-03-04 01:47:46,024 src.run_editing GIREMI failed!
ERROR 2019-03-04 01:47:46,024 src.run_editing No alignment file /work_dir/data/D8.sorted.bam
Traceback (most recent call last):
File "/usr/local/bin/run_rnacocktail.py", line 931, in
sys.exit(run_pipeline(args,parser))
File "/usr/local/lib/python2.7/dist-packages/src/main.py", line 242, in run_pipeline
workdir=args.workdir, outdir=args.outdir, timeout=args.timeout)
File "/usr/local/lib/python2.7/dist-packages/src/run_editing.py", line 372, in run_editing
raise Exception(excp)
Exception: No alignment file /work_dir/data/D8.sorted.bam
from rnacocktail.
It seems that your bam is not visible to the code.
Can you check you have your bam in your local machine at /D8/D8.sorted.bam
and you can see that when you load it to docker using:
docker run --rm -u container_R -v /D8:/work_dir/data -v /Container:/work_dir/write rnacocktail:latest ls /work_dir/data/D8.sorted.bam
from rnacocktail.
Yes, I did see its existence.
$ docker run --rm -u container_R -v /D8:/work_dir/data -v /Container:/work_dir/write marghoob/rnacocktail:latest ls /work_dir/data/D8.sorted.bam
/work_dir/data/D8.sorted.bam
from rnacocktail.
Can you try the following command, so that I can see full details about you bam file:
docker run --rm -u container_R -v /D8:/work_dir/data -v /Container:/work_dir/write rnacocktail:latest ls -al /work_dir/data/D8.sorted.bam
from rnacocktail.
docker run --rm -u container_R -v /D8:/work_dir/data -v /Container:/work_dir/write rnacocktail:latest ls -la /work_dir/data/D8.sorted.bam
lrwxrwxrwx 1 1012 NGS-BI-RNA 105 Feb 26 06:22 /work_dir/data/D8.sorted.bam -> /D8/D8.sorted.bam
id container_R
uid=3001(container_R) gid=1001(NGS) groups=1001(NGS),1305(NGS-BI-RNA)
from rnacocktail.
Thank you for your continued help
from rnacocktail.
I think the problem is the symlink. You should have absolute path in docker.
Can you check ls -al /D8/D8.sorted.bam
in your host system?
from rnacocktail.
It's works. The data is indeed a soft link. But I have a new problem again. “Error: Unable to access jarfile GenomeAnalysisTK.jar” I checked the dockerfile and did not install GATK, right?
from rnacocktail.
Sounds good. Yes, GATK is not part of our docker.
You can download GATK 3.5 at:
https://software.broadinstitute.org/gatk/download/archive
https://software.broadinstitute.org/gatk/download/auth?package=GATK-archive&version=3.5-0-g36282e4
and provide the jar file as input like --gatk /path/to/GenomeAnalysisTK.jar
from rnacocktail.
Related Issues (12)
- Exception: Unable to detect format from ['SNV;ENSG00000225630', '1', '+', '50', '0'] HOT 10
- Variant Calling HOT 2
- Error in run_editing STEP 10 HOT 3
- Keep moving HOT 4
- docker test fails on fusioncatcher in 0.3.0 HOT 1
- Transcript quantification of pacbio long reads (ccs)
- option '-k' is only recognised by salmon quant of salmon version 0.11.0 but not above
- samtools version 1.3 or higher incompatible HOT 1
- IOError: [Errno 36] File name too long: HOT 2
- hisat2_jun2bed.py error HOT 1
- Cannot start run_rnacocktail.py HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from rnacocktail.