Comments (3)
Hi Fabian,
that fungus does not really look very monoploid to me. It actually look like a very clear tetraploid. Is that even possible?
Monoploid genomes show only one coverage peak, but your genome shows 3 to 4. And smudgeplot shows lots of k-mer pairs - namely k-mer pairs with ~150x and 50x coverages (the unannotated smudge) but also 100x and 100x (what is annotated as AB smudge). Given the gneomescope model and positions of the two smudges, I think the correct annotations would be AAAB and AABB respectively. When you do the plot, you can give a coverage prior 54 and the smudges should get correct labels.
Finally, it does not look like a population variability either. If there would be a genetic diveristy of individuals in the isolate, the peaks would not be spaced so evenly and the coverage ratios would not match precise ratios (0.25 and close to 0.5). So this genomic dataset looks terribly tetraploid to me.
Sorry about the misslables, I need to revisit the annotation algorithm soon, I have seen quite a few smudgeplots that were clearly misslabeled.
from smudgeplot.
Really when using coverage before 54, the stains must receive correct labels. I still need to think about the tetraploidy of this strain:
I noticed that the case is very similar to the Saccharomyces tutorial (by the way, this is a Saccharomyces).
https://github.com/KamilSJaron/smudgeplot/wiki/tutorial-saccharomyces
from smudgeplot.
Oh cool, Sacharomyces quite wild when it comes to ploidy levels.
By the way, the last of any AB smudge would suggest this is not degenerated polyploid at all. I.e. there are no sequences that would be private only to two out of four haplotypes.
I will close it for now, feel free to reopen if anything.
from smudgeplot.
Related Issues (20)
- Haploid x Smudgeplot HOT 1
- Lower coverage cutoff higher than upper coverage cutoff HOT 1
- sharing strange smudgeplot for imporve smudgeplot HOT 1
- Smudgeplot for allotetraploidy HOT 5
- Interpretation of smudgy smudgeplot for diatom genome HOT 2
- Could I use to these AAB kmer to make evaluation the assembled phased genome? HOT 3
- failed to obtain the final plot HOT 3
- Smudgeplot estimated incorrect ploidy level HOT 6
- Outputting smudgeplot in pdf format HOT 1
- What is "1n = 25" on the smudgeplot plots? HOT 2
- Smudgeplot and genomescope HOT 4
- Problems installing smudgeplot HOT 3
- Can not find R package 'smudgeplot' HOT 1
- Interpreting smudgeplot output HOT 10
- Question about kmc parameters HOT 3
- issue about smudge_pairs HOT 3
- Smudgeplot is not suggesting an expected ploidy level HOT 2
- Salix RAD-data smudgeplots HOT 6
- estimate_1n_coverage_1d_subsets averages over two distinct high peaks
- Fix estimate_1n_coverage_highest_peak to take histogram 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 smudgeplot.