Closed luckysardar closed 1 year ago
Hi, you don't need to rebuild the docker image for using the pipeline. The pipeline will download it for you from dockerhub. Sometimes rebuilding the images some dependencies are moved from repositories etc. So better to use the image we generated.
Thanks for the quick response,
I Have tried to run the pipeline
nextflow run make_anno.nf -with-docker -bg --vcffile mgp.v5.merged.snps_all.dbSNP142.vcf.gz --speciesA CAST_EiJ --speciesB 129S1_SvImJ --genome GRCm38_68.fa --outvcf CAST_EiJ-129S1_SvImJ.vcf > log
log file
N E X T F L O W ~ version 22.10.0 Launching
make_anno.nf[mighty_ekeblad] DSL2 - revision: ea836071c0 Cannot enable more than one container engine -- Choose either one of: docker, singularity
Uh, you are right I set singularity as default. I changed it. Can you do a git pull and try again? Thanks for spotting the error.
Okay Thanks i will do that and get back to you
Hello I have run the same command which was mentioned above gave me warning and no output was generated
log file
N E X T F L O W ~ version 22.10.0
Launching `make_anno.nf` [tiny_raman] DSL2 - revision: ea836071c0
BIOCORE@CRG RNAseq - N F ~ version 1.0
╔═╗┬ ┬ ┌─┐┬ ┌─┐ ╔═╗┌─┐┌─┐┌─┐┬┌─┐┬┌─┐ ╦═╗╔╗╔╔═╗┌─┐┌─┐┌─┐
╠═╣│ │ ├┤ │ ├┤ ╚═╗├─┘├┤ │ │├┤ ││ ╠╦╝║║║╠═╣└─┐├┤ │─┼┐
╩ ╩┴─┘┴─┘└─┘┴─┘└─┘ ╚═╝┴ └─┘└─┘┴└ ┴└─┘ ╩╚═╝╚╝╩ ╩└─┘└─┘└─┘└
====================================================
BIOCORE@CRG Allele Specific RNAseq - N F ~ version 1.0
====================================================
vcffile : mgp.v5.merged.snps_all.dbSNP142.vcf.gz
speciesA : CAST_EiJ
speciesB : 129S1_SvImJ
genome : GRCm38_68.fa
outvcf : CAST_EiJ-129S1_SvImJ.vcf
email for notification : mymail@mydomain.eu
warning
=============================================================================
= WARNING =
= You are running this script using DSL2 syntax, however it does not =
= contain any 'workflow' definition so there's nothing for Nextflow to run. =
= =
= If this script was written using Nextflow DSL1 syntax, please add the =
= setting 'nextflow.enable.dsl=1' to the nextflow.config file or use the =
= command-line option '-dsl1' when running the pipeline. =
= =
= More details at this link: https://www.nextflow.io/docs/latest/dsl2.html =
=============================================================================
Yes, the newer version of nextflow assumes is DSL2. So... you have to add -dsl1 option when running nextflow. I'll add a note to the documentation
Thank you Its running now i will let you know if i face anything
hello again I got a error some permission error i have given all permission to nextflow and try to run as root user also but process did not completed
N E X T F L O W ~ version 22.10.0
Launching `make_anno.nf` [evil_mandelbrot] DSL1 - revision: ea836071c0
BIOCORE@CRG RNAseq - N F ~ version 1.0
╔═╗┬ ┬ ┌─┐┬ ┌─┐ ╔═╗┌─┐┌─┐┌─┐┬┌─┐┬┌─┐ ╦═╗╔╗╔╔═╗┌─┐┌─┐┌─┐
╠═╣│ │ ├┤ │ ├┤ ╚═╗├─┘├┤ │ │├┤ ││ ╠╦╝║║║╠═╣└─┐├┤ │─┼┐
╩ ╩┴─┘┴─┘└─┘┴─┘└─┘ ╚═╝┴ └─┘└─┘┴└ ┴└─┘ ╩╚═╝╚╝╩ ╩└─┘└─┘└─┘└
====================================================
BIOCORE@CRG Allele Specific RNAseq - N F ~ version 1.0
====================================================
vcffile : mgp.v5.merged.snps_all.dbSNP142.vcf.gz
speciesA : CAST_EiJ
speciesB : 129S1_SvImJ
genome : GRCm38_68.fa
outvcf : CAST_EiJ-129S1_SvImJ.vcf
email for notification : mymail@mydomain.eu
Error executing process > 'parseVCF (mgp.v5.merged.snps_all.dbSNP142.vcf.gz)'
Caused by:
Unable to create directory=/Allelic_analysis_new/allele_specific_RNAseq/makeAnno/work/c8/93e60c21c12211cd17b010d93426c5 -- check file system permissions
Failed to invoke `workflow.onComplete` event handler
-- Check script 'make_anno.nf' at line: 102 or see '.nextflow.log' file for more details
Pipeline BIOCORE@CRG Master of Pore completed!
Started at 2022-10-21T18:01:15.783383+05:30
Finished at 2022-10-21T18:01:16.614244+05:30
Time elapsed: 830ms
Execution status: failed
is it anything to do with system config I have 64GB ram with 12 core
this is weird. Can you please add -w $PWD/work to the command line? It looks like is trying to write to a place where you have no privileges... Are you using Linux or mac?
I am using linux
Let me know if you manage
Hello I have manage to ran as sudo user it ran succesfully, Got the output for the one of the sample SRR14004348 when I am checking allelic proportion chromsome wide data is biased towords one of the allele also same for other sample can you please suggest me anything i am missing. I just took mean across chromsomes. followed same pipeline with input file mentioned
what we expect from autosome should give almost 50:50 ratio most of the chromosome showing biased average ratio is it because of some aneuploidy or what please suggest
Thank you
chromosome A/(A+B) B/(A+B)
X 0.855416115526518 0.144583884473482
13 0.504994567271183 0.495005432728817
9 0.436021416100318 0.563978583899682
7 0.395250986658452 0.604749013341547
5 0.373462992941653 0.626537007058346
8 0.366074201526175 0.633925798473825
2 0.230499212359011 0.76950078764099
11 0.204298039835776 0.795701960164223
6 0.11859084240173 0.88140915759827
4 0.037873537098041 0.962126462901958
3 0.030748955771279 0.96925104422872
16 0.026845202408981 0.973154797591019
18 0.023344626780433 0.976655373219567
14 0.022645613522724 0.977354386477276
17 0.021418953610684 0.978581046389316
1 0.021037927064929 0.978962072935071
19 0.01496138275413 0.985038617245871
15 0.013040629549786 0.986959370450214
10 0.011114480140262 0.988885519859738
12 0.008684188037739 0.991315811962261
MT 1.58665051548522E-05 0.999984133494845
Hey @luckysardar,
if you're saying that you're running this pipeline on this PGCLC day 5 XGFP- rep 3 sample than I guess the results you get are pretty good.
If I remember correctly only chr X and 13 are fully hybrid with good distribution of the Cas and Mus SNPs. You can check a distribution of the SNPs here. Anyway I'm pinging @Jacq3lin3 and @MoritzBauer (the paper first authors) as they know better.
So the allelic ratio of chr 13 is perfectly 50% and chr9 also very close, while chr X is not reflecting the X chromosome Inactivation (XCI) happening in XGFP- PGLCs.
Nicco
@lucacozzuto I created a while back a small sh script that counts the sum allele counts per file. I'd like to push it to the development branch but I don't have write access. I'd avoid creating a fork if you agree.
Yes. Please go ahead
Thank you @Ni-Ar @lucacozzuto for the clarification
Yes. Please go ahead
@lucacozzuto I think that biocoreCRG needs to grant me write access. I'm reading info here.
Cause from my GitHub Desktop I see this message when trying to push to the development branch.
Strange. You should be able to fork and to make a pull request
yes I could fork + PR. I wanted (if possible) to commit to the "development" branch directly.
Thank for the newflow pipeline
I am running
docker build .
getting bellow error