Comments (14)
change -t 8
to -t 64
perhaps
from haslr.
Thanks, i have runned the 55 threads shell and not break the 8 threads process. How long do you think it will take to run the both scripts
from haslr.
Did you change the output directory? I have no idea how long it might take? Depends on coverage of long and short reads?
from haslr.
sure,i set a new output dir
from haslr.
it is always running minia for over 24 hours ,is it normal?
from haslr.
Minia is very fast, but genome size and coverage influence its runtime also probably choice of k-mer length and other similar types of settings.
from haslr.
So,is there any recommended parameters for running human genome assembly?
from haslr.
Are you trying out the assembler with someone else's data or do you have a new human genome assembly that you would like to make with your own data? I would think that it would have finished by now (~5 days running). Again, you haven't specified the coverage of the Illumina or I guess Oxford Nanopore data that you are using. You can also read the paper describing HASLR for perhaps more information on the program.
from haslr.
Sorry,im trying to assembly a human genome, The coverage of both short reads and long reads is 120X
from haslr.
I would recommend you try either GraphAligner
(https://github.com/maickrau/GraphAligner) or Ratatosk
(https://github.com/DecodeGenetics/Ratatosk) to error correct your Nanopore reads with your Illumina reads then assemble with Flye
(https://github.com/fenderglass/Flye) using the --nano-corr
option. Ratatosk
even has a faster reference based method whereby to correct the reads (I haven't used this method, so I don't know the details). For Flye
I really don't think you need 120x Nanopore coverage, especially if you can correct the reads. See here for running Ratatosk
or here for running GraphAligner
.
Edit: I guess you could use 120x Nanopore reads for a Human assembly (https://github.com/fenderglass/Flye/blob/3ee5b3390a5f88c36d0869d0382c75aba3b1f5cc/README.md#flye-benchmarks), although these data come from CHM13 (homozygous cell line). Also note the 4000 CPU hours (divide 4000 by number of available cores and you get approximately how many wall hours the assembly would take).
from haslr.
Thanks,jelber2.
so haslr is not advised ?why?
from haslr.
In my experience, HASLR will generate very good statistics (N50, etc) for assembly using raw long reads and accurate short reads, but the error rate (indels and substitutions) of the final assembly is similar to the error rate of the long reads and not the short reads. One can improve the error rates by using long reads corrected by the short reads, and using the corrected long reads as input, but then the assembly statistics suffer. This is based off of simulation of course, and simulations are sometimes useful but can never fully capture the intricacies of real data.
from haslr.
Hi @koujiaodahan and thanks for trying HASLR.
I'm surprised that Minia is taking so long to finish. In my experience, on short read datasets from human genome with about 40x coverage, it takes about 5 hours to finish. Are you sure that Minia assembly was the step that took a long time to finish?
If yes, one solution could be subsampling short reads to about 40-50x coverage. You can use fastutils command that comes with HASLR for that purpose. So assuming you have a paired end dataset, you do the following:
fastutils interleave -q sr_1.fastq sr_2.fastq | fastutils subsample -q -g 3g -d 40 > sr_40x.fastq
With regards to the error rate of the final assembly that was raised by @jelber2, if you eventually want to perform polishing for your assembly, our results show that polished HASLR assemblies are as accurate as polished assemblies from other tools.
from haslr.
yeah, i agree that the coverage is too high,so i downsampled and i got error which i released at #20 .
and i want to know whether your polishing method means running wtdbg2.pl after running haslr?
from haslr.
Related Issues (20)
- assembling long reads using HASLR... failed ERROR HOT 9
- Debug Message in log File
- Bug while running script: "could not find file"
- Short Homo Sapien Assembly from Genome in a Bottle Data HOT 10
- Multiple long read files in the latest version
- Changing K-mer size of Minia leads to final fasta of zero bits
- Failing while checking minia
- question multiple libraries?
- Does this tool do phasing?
- kmer determination
- Much Shorter Assembly than Expected HOT 1
- AttributeError: 'module' object has no attribute 'run' HOT 1
- lr0x.fasta... [ERROR] option -d/--depth is required HOT 1
- new release for bioconda
- Assembling pacbio, nanopore and illumina reads HOT 2
- ERROR: "haslr_assemble" returned non-zero exit status
- Clean raw data prior to assembly
- Assembling the test data does not generate an assembly HOT 2
- final.fa is 0 kb HOT 2
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 haslr.