Comments (2)
No worky, here's the approach and outputs
make new environment
conda create -n bayescmg
activate it
conda activate bayescmg
run the run.sh with all parameters
bash ~/bin/VarBayes/run.sh -v 17-10-13_ARUP-Marth-Subependymal_Heterotopia_Final_1511985898.vcf.gz -p 17-10-13_ARUP-Marth-Subependymal_Heterotopia.ped -r /scratch/ucgd/lustre/work/u0691312/reference/ucgd_reference/GRCh37/human_g1k_v37_decoy_phix.fasta -c -g /scratch/ucgd/lustre/work/u0691312/reference/slivar.gnomad.hg37.zip -d /scratch/ucgd/lustre/work/u0691312/reference/ensembl/ -u /scratch/ucgd/lustre/work/u0691312/reference/ensembl/Plugins/ -l /scratch/ucgd/lustre/work/u0691312/reference/ensembl/Plugins/revel_all_chromosomes_vep.tsv.gz -t 0.01 -j 0.6 -y 0.1 -o 350 -e 0.1 -f ./
wget: /uufs/chpc.utah.edu/common/HIPAA/u0691312/bin/miniconda3/lib/libuuid.so.1: no version information available (required by wget)
--2020-03-15 10:50:08-- ftp://ftp.ncbi.nlm.nih.gov/pub/clinvar/vcf_GRCh37/clinvar.vcf.gz
=> ‘data/clinvar.grc37.vcf.gz’
Resolving ftp.ncbi.nlm.nih.gov (ftp.ncbi.nlm.nih.gov)... 130.14.250.10, 2607:f220:41e:250::11
Connecting to ftp.ncbi.nlm.nih.gov (ftp.ncbi.nlm.nih.gov)|130.14.250.10|:21... connected.
Logging in as anonymous ... Logged in!
==> SYST ... done. ==> PWD ... done.
==> TYPE I ... done. ==> CWD (1) /pub/clinvar/vcf_GRCh37 ... done.
==> SIZE clinvar.vcf.gz ... 25088270
==> PASV ... done. ==> RETR clinvar.vcf.gz ... done.
Length: 25088270 (24M) (unauthoritative)
100%[=========================================================================================================================================================================>] 25,088,270 21.2MB/s in 1.1s
2020-03-15 10:50:12 (21.2 MB/s) - ‘data/clinvar.grc37.vcf.gz’ saved [25088270]
wget: /uufs/chpc.utah.edu/common/HIPAA/u0691312/bin/miniconda3/lib/libuuid.so.1: no version information available (required by wget)
--2020-03-15 10:50:12-- ftp://ftp.ncbi.nlm.nih.gov/pub/clinvar/vcf_GRCh37/clinvar.vcf.gz.tbi
=> ‘data/clinvar.grc37.vcf.gz.tbi’
Resolving ftp.ncbi.nlm.nih.gov (ftp.ncbi.nlm.nih.gov)... 130.14.250.10, 2607:f220:41e:250::11
Connecting to ftp.ncbi.nlm.nih.gov (ftp.ncbi.nlm.nih.gov)|130.14.250.10|:21... connected.
Logging in as anonymous ... Logged in!
==> SYST ... done. ==> PWD ... done.
==> TYPE I ... done. ==> CWD (1) /pub/clinvar/vcf_GRCh37 ... done.
==> SIZE clinvar.vcf.gz.tbi ... 283593
==> PASV ... done. ==> RETR clinvar.vcf.gz.tbi ... done.
Length: 283593 (277K) (unauthoritative)
100%[=========================================================================================================================================================================>] 283,593 1.00MB/s in 0.3s
2020-03-15 10:50:13 (1.00 MB/s) - ‘data/clinvar.grc37.vcf.gz.tbi’ saved [283593]
/uufs/chpc.utah.edu/common/HIPAA/u0691312/bin/VarBayes/run.sh: line 131: vep: command not found
/uufs/chpc.utah.edu/common/HIPAA/u0691312/bin/VarBayes/run.sh: line 152: bgzip: command not found
/uufs/chpc.utah.edu/common/HIPAA/u0691312/bin/VarBayes/run.sh: line 153: tabix: command not found
/uufs/chpc.utah.edu/common/HIPAA/u0691312/bin/VarBayes/run.sh: line 156: externals/slivar/slivar: No such file or directory
/uufs/chpc.utah.edu/common/HIPAA/u0691312/bin/VarBayes/run.sh: line 162: vep: command not found
/uufs/chpc.utah.edu/common/HIPAA/u0691312/bin/VarBayes/run.sh: line 183: bgzip: command not found
/uufs/chpc.utah.edu/common/HIPAA/u0691312/bin/VarBayes/run.sh: line 184: tabix: command not found
python VarBayes.py -v data/slivar.tmp.vep.vcf.gz -f 17-10-13_ARUP-Marth-Subependymal_Heterotopia.ped -d ./ -c data/clinvar.grc37.vep.vcf.gz -e 0.1 -o 350 -p 0.1 -a 0.01 -r 0.6
python: can't open file 'VarBayes.py': [Errno 2] No such file or directory
from bayescmg.
We now check for these dependencies and exit if they are not present.
from bayescmg.
Related Issues (15)
- Write scripts for parsing on PostP and other useful things HOT 1
- Make sure BP7 is VEP Consequence is synonymous and NOT splice_region and check counting the number of synonymous variants HOT 1
- Update github documentation and README HOT 2
- run.sh crashes even with all necessary components in PATH HOT 3
- Rename everything VaBayes to bAyesCMG HOT 1
- Fix run.sh message HOT 1
- run.sh should not download clinvar and other data to current directory HOT 3
- Add variant quality filtering using slivar
- Check variant logic being only applied to variants present in proband
- Skip (pre) processing steps if possible
- Add a flag to keep intermediate files or not HOT 1
- VQSRTranche variants being retained HOT 1
- Need to annotate VCF before slivar comp het function HOT 1
- Getting multiple clinvar tbis HOT 1
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 bayescmg.