Comments (3)
Just started exploring consensus.pl myself. Looks to be suspect. Help documentation is wrong. From what I can tell the program expects (1) reference file, (2) bed file, (3) bam file which explains tallnuttcsiro's errors. Program does not catch missing files. Even when files are given it can not open gzipped file. Has not been updated since 2013 thus my confidence in the program is low,
from fasta_utilities.
No need to be inflammatory @westerman.
As you say the script needs three input files, I apologize for the documentation inconsistency.
It can open a gzipped or bzipped file if you have pigz or pbzip on your PATH. (see scripts/consensus.pl#L25-L29).
It also requires that wrap.pl
is in ~/fu
which it will not be most cases. (scripts/consensus.pl#L87) and is likely the reason you are both having difficulty. That line should be changed to look for wrap.pl
on the PATH rather than hard-coded.
This also assumes mpileup output from ~2013, which is pre-HTSlib, so I cannot be sure the format has remained the same.
I would probably not recommend using this script directly, but instead use it as a reference on parsing the output and adapt to your needs.
from fasta_utilities.
Wasn't trying to be inflammatory. Just pointing out that the help is not correct, that I did not (and cannot) get a gzip file to open, and that the program hasn't been updated since 2013. Thus, for all I know, the program may be unusable and abandoned. Heavens knows I have some programs of my own like that -- still in my overall path but ones I never use because they are broken.
As for the gzip file, below is my command line, the response, and where pigz is located:
]> consensus.pl Carpophilus.fa 017150_c-yuc-GATK-HaplotypeCaller.vcf.gz 017150_c-yuc.sorted.bam
No such file or directory:Could not open pigz -dc < 017150_c-yuc-GATK-HaplotypeCaller.vcf.gz|
which pigz
/usr/bin/pigz
Running the command 'pigz -dc < 017150_c-yuc-GATK-HaplotypeCaller.vcf.gz' without the pipe outside of 'consensus.pl' works fine. Also works with 'more' or 'cat' after the pipe. But is not working within consensus.pl.
Thanks for your work on the rest of the utilities.
from fasta_utilities.
Related Issues (4)
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 fasta_utilities.