Comments (11)
from quickmerge.
Sure. First I sourced the .quickmergerc, then:
for the wrapper:
merge_wrapper.py hybrid_assembly.fasta self_assembly.fasta
which finished with no errors.
Then with individual steps as:
nucmer -l 100 -prefix out self_assembly.fasta hybrid_assembly.fasta
delta-filter -i 95 -r -q out.delta > out.rq.delta
quickmerge -d out.rq.delta -q hybrid_assembly.fasta -r self_assembly.fasta -hco 5.0 -c 1.5 -l 400000 -ml 5000
which also didn't show any errors. Both resulted in the same output content described above.
from quickmerge.
from quickmerge.
Yes. Would it be useful to upload this file?
from quickmerge.
from quickmerge.
You can download the file using the following link:
https://www.dropbox.com/s/5qlld3k2m0n4da4/out.rq.delta.gz
And the 400kb n50 is from the self-assembly file.
from quickmerge.
from quickmerge.
The reads almost certainly derive from the same species. One point that could be important is that the self_assembly.fasta is almost 50% larger than the expected genome size, which is probably the result of redundancy in the assembly. Do you think this redundancy could result in the type of whole genome alignment results being observed here?
from quickmerge.
If there is high homology between the two assemblies then it should show HCO >5. Will you be able to share the fasta files so that we can investigate this a little more. We will use your files only to understand the problem and will not share it with anyone else. Let us know.
from quickmerge.
Do you have a preferred email address where the files could be sent to?
from quickmerge.
from quickmerge.
Related Issues (20)
- Unexpected sequence in merged fasta HOT 2
- quickmerge ( multiple long-read only ) + ( multiple hybrid ) assemblies
- ERROR: mummer and/or mgaps returned non-zero
- Does running iterations of quickmerge on quickmerge output create artifacts? HOT 3
- Quickmerge before or after haplotype purging for high heter rate plant assemblies? HOT 1
- quickmerge: terminate called after throwing an instance of 'std::out_of_range' HOT 6
- Is there any problem to use quickmerge and give in input a previous output of quickmerge?
- Problem with Reference and Query
- segment fault and Core dump HOT 3
- Confusion about argument ```-l``` and ```-ml``` HOT 3
- [Mac] Running merge_wrapper.py in conda : No such file or directory: 'nucmer': 'nucmer' HOT 3
- quickmerge aborted HOT 8
- Is there a history of merges not included in anchor_summary.txt? HOT 1
- Output problems HOT 2
- sometimes quick merge does not merge contigs? HOT 3
- terminate called after throwing an instance of 'std::out_of_range' HOT 7
- Quickmerge fails to set its PATH internally HOT 6
- core dump at final step HOT 3
- Aborted (core dumped) quickmerge HOT 5
- Program runs, but doesn't provide correct output 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 quickmerge.