Comments (6)
Hi, the function would need to be modified since it currently requires the .scaffolds.tdt file right now. Here's an example below of what that looks like.
Qry Scaffold Description QryLen QryStart QryEnd Strand Ref RefLen RefStart RefEnd RefMap Identity Length Coverage N Rank Inv Trans
query1.01 1 1 len=51.46 Mb 1.03% 1(+) 74,333:77,337,987; 0.90% 1(-); 8.98% other; 51464084 7953 51441722 + 1 77392008 74333 77337987 Anchored 485780 541756 531399.0 446 1 465116.0 4619909.0
query1.02 1 1 len=43.32 Mb 8.98% 1(+) 5,070:77,338,662; 1.21% 1(-); 12.20% other; 43319783 18524 43288930 + 1 77392008 5070 77338662 Placed 3613688 4005849 3890364.0 571 1 525171.0 5283668.0
It might be a while before we could work on that, so in the meantime if you think you can just spoof that format that might be the fastest way.
from clinehelpr.
Thank you very much. It looks tricky, but I'll try.
from clinehelpr.
Hi,
If you have a fasta file for the reference genome and another fasta file for your query scaffolds, You can run your fasta files through minimap2 to generate a PAF output file, and then use the PAF file as input into PAFScaff. PAFScaff then generates the scaffolds.tdt file that you can use with the Ideograms.
Is that what you had in mind, or am I perhaps misunderstanding what you were wanting to do?
-Bradley
from clinehelpr.
Thank you for your response, Bradley.
Maybe I didn't make myself clear; I mapped the reads for bgc to the draft sequence of the study species with a gff file (almost at the chromosome level) and use this information instead of the reference genome of the other relative species. Perhaps I should also look around other packages to draw an ideogram.
from clinehelpr.
Oh ok. Now I understand.
I'm going to mark this as a requested feature that we will try to get to in the future.
In the meantime, if you don't want to use the PAFScaff files to make the ideogram, you can try to write your own script to use the RIdeogram package (https://cran.r-project.org/web/packages/RIdeogram/vignettes/RIdeogram.html), which is what we used under the hood to make the ideogram plots. In the script you'd just have to get your data into a format that ideogram accepts.
-Bradley
from clinehelpr.
Thank you very much, Bradley. I was thinking about the rideogram. The ClineHelpR helped me a lot. Thank you again. -Makiko
from clinehelpr.
Related Issues (17)
- Hybrid Index HOT 1
- Error with parsing alleles depth HOT 7
- Error HOT 9
- get_bgc_outliers error HOT 15
- Rename snps$crazy.a and snps$crazy.b to something more intuitive
- Rebuild docker image with updates to vcf2bgc.py HOT 1
- Fix run_bgc.sh estpost LnL file naming HOT 1
- run bgc error in docker HOT 3
- ClineHelpR R package conda install? HOT 1
- jinja2 HOT 2
- R function: combine_bgc_output HOT 2
- GSL error. HOT 2
- phiPlot visualization issue HOT 1
- Assessing the convergence of runs HOT 3
- Add vcf2bgc script
- vcf2bgc.py file HOT 6
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 clinehelpr.