Comments (3)
P_ctg consists of one set of haplotype but its contigs might be switched between two haplotypes. I guess it is usually used as reference. If the paternal short reads are available, hifiasm can output fully phased assembly, i.e., two fully phased haplotypes. We think fully phased assembly should be also generated with HiC or Strand-seq if paternal data is not available, but for now hifiasm itself cannot automatically do that. If you can get phased partitioned list of HiFi reads with HiC, then feeding it to hifiasm using option '-3/-4' can also output fully phased assembly.
from hifiasm.
@chhylp123 Thanks very much for the reply! I see, some contigs might be switched..
btw, P_ctg is one full set of halpotype, but a_ctg is not, I think it contains only the alternative regions, I wondered how can I add the collapsed consensus regions into a_ctg as well, then I could obtain a full set of alternative haplotype (?).
from hifiasm.
It is hard. If you'd like to do it manually, probably: 1) run hifiasm with -l0, 2) find homozygous unitigs at r_utg with homozygous coverage, 3) add these unitigs to a_ctg.
from hifiasm.
Related Issues (20)
- Mitigate Overlapping Sequence Assignments in Haplotypes HOT 3
- Help!!! Segmentation fault (core dumped) HOT 1
- Question about the depth of ONT ultra-long reads HOT 1
- Homotetraploid, super-large genome, with different parameters, the size of p_utg varies greatly? HOT 1
- setting K parameter in yak HOT 2
- how to make the correct genome size estimation for allotetraploid species? HOT 2
- Possible missing one haplotype in human assemblies HOT 2
- No haploid.gfa files output in trio-binning mode HOT 3
- Hifi + Hi-c + ONT assembly fails
- In Trio-binning, always more on hap1 despite (almost) same sequences for paternal and maternal
- discontinuous assembly with shorter pacbio hifi reads but high coverage HOT 2
- Is x20 of Hifi data enough to construct draft assembly of 6.5Gb genome? HOT 1
- line 8: 110334 Aborted(core dumped) HOT 1
- Ultra Long intergration failed: no output for UL kmer counting HOT 3
- missing 8Mb sequences in the assembly HOT 5
- Empty haplotype 2 gfa files by ONT integration HOT 1
- Basic Question About HiFi Input HOT 3
- Spend too long times to run hifiasm HOT 1
- Switch error on X and Y chromosome HOT 2
- *.ovlp.bin file 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.