Comments (2)
--n-hap
is used to determine the coverage of heterozygous nodes or contigs. For your sample, hifiasm thinks the homozygous coverage is 26, and the heterozygous coverages are 26/2 = 13 and 26/4 = 6 using --n-hap 2
and --n-hap 4
, respectively. Hifiasm keeps any node in the assembly graph with coverage above the heterozygous coverage threshold as a real node, instead of sequencing errors. This is why --n-hap 4
leads to a larger graph. Could you please have a try with --hom-cov 55
and --n-hap 2
? Since bv looking at the k-mer plot, there are only two peaks and the homozygous coverage should be 55.
from hifiasm.
Dear author,
I tried your suggestions, and here are the results:
hifiasm --n-hap 2 -t 120 -l 0 --hom-cov 55 the generated .p_ctg.gfa file is of size 56G, and the .p_utg.gfa file is of size 66G.
Since mine is a homologous tetraploid, which form should I choose for assembly, p_ctg or p_utg?
p_ctg N50: 100MB
p_utg N50: 1MB
I believe that increasing the depth of HiFi data will not increase the length of N50.
For the above question, could you provide some debugging suggestions? Thank you for your valuable time and assistance. I sincerely look forward to your response!
from hifiasm.
Related Issues (20)
- Spend too long times to run hifiasm HOT 4
- Switch error on X and Y chromosome HOT 2
- *.ovlp.bin file HOT 1
- Resolving switching error (?)
- Interchromosomal misjoin HOT 1
- Read error correction does not reduce the number of kmers present once, twice or three times
- Recreate p_ctg from p_utg HOT 3
- In the diploid assembly, hfiiasm identified a value that did not exist in the k-mer plot as the "homozygous read coverage threshold".
- fungi diploid assemble phasing errors
- What is the strategy for Ultra-long ONT integration? HOT 3
- No further progress after generating *bin files HOT 1
- Issues with Hexaploid Assembly
- Numerous spurious contig overlaps in assemblies HOT 1
- output GFA with J lines in scaffolding mode
- Output that maps repeat shrapnels to unitigs and to scaffolded chromsomes
- High Duplication in my Assembly...l3 (purge-dup) doesn't help much HOT 1
- How to classify hifi reads into haplotype genomes? minimap2 did not work well HOT 2
- False inserts over duplicated regions (TRBV12-3/4) HOT 2
- Whether unitig will be filtered based on the hic signal of unitig HOT 2
- Hifiasm Assembly Purging 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 hifiasm.