Comments (2)
Hi, can you please address my inquiry above?
Thanks!
from hg-color.
Hello,
Sorry for not answering before, I was in Christmas break and wanted to take a real work-free break before I start writing my thesis.
PE libraries don't matter to much to HG-CoLoR, as it does not make use of that information. However, experiences I ran shown that using smaller short reads (125bp over 250-300 bp) provided slightly better results. This results are shown in the Supplementary Material of the paper (https://oup.silverchair-cdn.com/oup/backfile/Content_public/Journal/bioinformatics/34/24/10.1093_bioinformatics_bty521/1/bty521_supplementary_data.pdf?Expires=1547204857&Signature=tmdkoxqJvTi84m82mOwoLHRaNEg4M5sfoqPVF48xQAsUVs5d10DZQB2qAWjLFXmKMu5DYif6LfZ65p69fPHLhSAU81ygTrravdxft2GQJXwZXj7fg~sdNUd5BxuK8EcTfttc2dkCmQNysecrRStshrT5TAZweMo-n22DAEym9RDqlNAdMJf0B2A9LaUs-o-l24Nscy5rH-icSa9nsUoZCMpuSChp8Ttfm28YeWgXx~x2m4Q-Hexs~0rfopRRk9MnWGJ~AIHLnRX7YAF~GQtUmX8YDE-EkKrY6Mj~eZUkV9tDvr-4PRsH611nzV73x1WUwXlyICtrVZ3ND4a2n9ihCA__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA) Tables S6, S7, and S8. Only tested the tool up to the C. elegans genome in these experiments however, so picking the larger 260 bp SRs might be helpful to better cover repeats in your case.
No need to trim the SR, nor remove the overlapping parts or anything. The HG-CoLoR pipeline will itself correct the SR (with QuoRUM, which is fast) to achieve high quality correction.
As for the coverage, I usually use about 50x coverage of SR, using more than that did not shown significantly better results in my experiments. If you still wish to use a higher SR coverage however, I highly recommand you to lower the --bestn parameter, and increase the --solid parameter to avoid prohibitive runtimes and miscorrections.
For the SR/LR alignment step, a known issue is also that BLASR does not allow the reference file (in the case, the LR file) to be higher than 4Go. You will thus have to split your LR file into separate 4Go files, and run separate HG-CoLoR instances on each file. This will not affect the correction results, as each LR is processed independently during correction. I know that this is however quite impractical, and investigating to find a better aligner is on my TODOlist.
Moreover, another known issue is that the graph construction with PgSA tends to take very, very long time when the size of the SR file grows bigger, as it does not support multithreaded construction. So using a 50x coverage for your plant genome might take quite some time. Again, this is a known issue, and replacing PgSA with a proper FM-index, allowing parallel construction, is also on my TODOlist.
Can't promise you when the update will be done, as writing my thesis currently takes quite a lot of my time. Demands for running HG-CoLoR on large genomes however become frequent, so I might take so time to do it the next few weeks, as this is a blocking point for every large experiment.
Best,
Pierre
from hg-color.
Related Issues (17)
- Error messages HOT 12
- runtime estimation HOT 3
- Error HOT 2
- Error (core dumped) HOT 4
- What is the proper output file HOT 2
- Error While Running QuorUM HOT 5
- compile error: fatal error: test/testdata.h: No such file or directory HOT 1
- tmp files deleted before job is finished HOT 2
- Running time? HOT 4
- Run fails and only prints 'Correcting the short reads' HOT 2
- Serious safety problem, --tmpdir parameter is deleted with no checking HOT 2
- Fail due to erroneous string in alignment file? HOT 1
- avoid reverse complement files HOT 1
- blasr doesn't use multiple threads? HOT 1
- Error HOT 1
- Disk quota exceeded when there is enough space on disk HOT 6
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 hg-color.