Comments (3)
Hi,
Thank you for your interest!
The reported alignment error rate is for the alignment of raw reads on the current assembly. Even if the assembly is polished, you would expect 12-16% errors coming from the reads, so the numbers you have do not necessarily mean that the genome requires more polishing.
Still, you can make an extra iteration and check if the statistics improve (most likely, for a tiny bit): use '--iterations 3' and '--resume' options to add one polishing iteration to the existing assembly.
from flye.
Hi,
I used the newest version of Flye (2.7-b1585) to handle nanopore reads (~60X, with default parameters except --asm-coverage 30). The assembly (N50: 14.7Mb) is much better than the one produced by Canu. Its Busco reaches almost 90%. The genome can annotate all of hox genes(39 genes) using homology-based prediction. But there are a lot of frameshift mutations in 38 genes.
How to know how many iterations I need to polish? how about the Busco score?
I also check if it's necessary to do more polishes with Racon. However, racon makes the quality of genome assembly worse. Busco score decreases to 86% and more hox genes get 'lost'.
Any suggestions are welcome.
Thanks a lot.
from flye.
Polishing methods are now rapidly developing - it is hard to give a concrete advice on the best strategy. I'd suggest to look into recent literature - this might be a good start: https://genomebiology.biomedcentral.com/articles/10.1186/s13059-019-1727-y
from flye.
Related Issues (20)
- ERROR: No disjointigs were assembled HOT 1
- Question for my understanding of the tool HOT 1
- `--deterministic` flag should set a random seed for random processes HOT 1
- ERROR: Looks like the system ran out of memory (flye) HOT 4
- "Caught unhandled exception: std::bad_alloc" error, even for example file HOT 1
- what's the meaning of repeat and mult. columns in assembly_info.txt with metagenomic data? HOT 2
- small deletion errors in assemblies HOT 1
- co-assemble samples with different read depth HOT 3
- Provide potential sequence circularity in Fasta headers HOT 1
- ERROR: Error running minimap2 HOT 10
- flye-minimap2 error HOT 4
- How can I aseemble Plasmid and Chromosome separately? HOT 4
- Assembling with "high" heterozygosity HOT 3
- [Question] How to estimate -g genome-size for --meta mode? HOT 3
- OSError: AF_UNIX path too long HOT 2
- 2.9.3 docker crashes HOT 2
- Quetion about the assembly_info.txt file HOT 1
- died with <Signals.SIGILL: 4 HOT 2
- Best parameters for ONT reads corrected with Illumina before assembly HOT 1
- Assembling ONT Multiple Chromosomes Form a Cluster 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 flye.