Comments (5)
The new version of HATCHet includes a new step called SNPCalling
and it has a slightly different structure. We are working on updating the documentation and in the meanwhile could you please check out the complete demo available at:
from hatchet.
The new version of HATCHet includes a new step called
SNPCalling
and it has a slightly different structure. We are working on updating the documentation and in the meanwhile could you please check out the complete demo available at:
Started with the demo. Downloaded the files and now execute runUnphased.sh.
In the SNPcaller-step I get a couple of errors all like this:
Process Caller-5: Traceback (most recent call last): File "/home/zyto/unger/anaconda3/envs/Hatchet/lib/python3.7/multiprocessing/process.py", line 297, in _bootstrap self.run() File "/home/zyto/unger/anaconda3/envs/Hatchet/lib/python3.7/site-packages/hatchet/utils/SNPCaller.py", line 121, in run snps = self.callSNPs(bamfile=next_task[0], samplename=next_task[1], chromosome=next_task[2]) File "/home/zyto/unger/anaconda3/envs/Hatchet/lib/python3.7/site-packages/hatchet/utils/SNPCaller.py", line 138, in callSNPs raise ValueError(sp.error('SNP Calling failed on {} of {}, please check errors in {}!').format(chromosome, samplename, errname)) ValueError: SNP Calling failed on chr5 of Normal, please check errors in /home/zyto/unger/WES_Rezidiv_HNSCC_Clonality/Hatchet/hatchet/examples/demo-complete/snps/Normal_chr5_bcftools.log!
Here the Normal_chr5_bcftools.log:
E::bgzf_read] Read block operation failed with error 4 after 44546 of 56208 bytes Failed to open https://ftp.ncbi.nih.gov/snp/organisms/human_9606_b151_GRCh37p13/VCF/GATK/00-All.vcf.gz: could not load index
Any suggestion?
from hatchet.
The new version of HATCHet includes a new step called
SNPCalling
and it has a slightly different structure. We are working on updating the documentation and in the meanwhile could you please check out the complete demo available at:
complete demoStarted with the demo. Downloaded the files and now execute runUnphased.sh.
In the SNPcaller-step I get a couple of errors all like this:
Process Caller-5: Traceback (most recent call last): File "/home/zyto/unger/anaconda3/envs/Hatchet/lib/python3.7/multiprocessing/process.py", line 297, in _bootstrap self.run() File "/home/zyto/unger/anaconda3/envs/Hatchet/lib/python3.7/site-packages/hatchet/utils/SNPCaller.py", line 121, in run snps = self.callSNPs(bamfile=next_task[0], samplename=next_task[1], chromosome=next_task[2]) File "/home/zyto/unger/anaconda3/envs/Hatchet/lib/python3.7/site-packages/hatchet/utils/SNPCaller.py", line 138, in callSNPs raise ValueError(sp.error('SNP Calling failed on {} of {}, please check errors in {}!').format(chromosome, samplename, errname)) ValueError: SNP Calling failed on chr5 of Normal, please check errors in /home/zyto/unger/WES_Rezidiv_HNSCC_Clonality/Hatchet/hatchet/examples/demo-complete/snps/Normal_chr5_bcftools.log!
Here the Normal_chr5_bcftools.log:
E::bgzf_read] Read block operation failed with error 4 after 44546 of 56208 bytes Failed to open https://ftp.ncbi.nih.gov/snp/organisms/human_9606_b151_GRCh37p13/VCF/GATK/00-All.vcf.gz: could not load index
Any suggestion?
Manually downloading did the trick...
from hatchet.
@brian-arnold please close the issue after updating the documentation with downloading and local usage instructions
from hatchet.
The config.sh file now allows users to specify the LIST variable, directing HATCHet to use a locally stored VCF of known germline SNP positions instead of HATCHet automatically downloading one from the internet. If LIST is left blank ("") in the config.sh file, then HATCHet attempts to download a list of SNP positions based on the reference genome version and chromosome annotation (which are specified by REF_VERS and CHR_NOTATION in config.sh).
from hatchet.
Related Issues (20)
- Running HATCHet with phasing
- Failed to liftover chromosomes with picard
- Failed to parse the output of the following command because the final objective was not found: {cmd} HOT 2
- Full documentation is not in fact available from cs.brown.edu HOT 3
- Error in read_bb() in cluster_bins.py due to change in pandas df.groupby() HOT 4
- Error in plot_clusters() in plot_cn_1d2d.py due to change in pandas df.groupby() HOT 1
- Error when no reads are recorded in tumor.bed before/after centromere location
- Specifying size for WES samples in count-reads HOT 1
- Set default value for ILP time limit
- can i ues this method for checking WGD of other species
- HATCHet2 RDR/BAF locally-aware clustering error with sparse bins
- String "None" is not converted to Python `None` object in `parse_combine_counts_args`
- Length of values does not match length of index
- resolve issue with CI testing python 3.7
- explore using more recent SHAPEIT version
- AssertionError at clustering step HOT 1
- add flag and make GC correction optional
- Reference panel requirement HOT 1
- About "download_panel"
- Issue runing Hatchet2 with hg38 reference_version
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 hatchet.