Comments (4)
I have not tested it in viruses yet.
Can you elaborate a bit more on what happened? Why do you think it is not "optimal", exactly? What were your expectations and what is the result you got?
from superpang.
Maybe the genome size of viruses is small, the core.fasta
file is empty for some viruses.
Such as Enterovirus A
from superpang.
I can reproduce with the HIV genomes that you mentioned in #6.
One problem is that by default SuperPang expects genomes to come from the same ANI>95% cluster. The ANI between the HIV genomes is lower. This can be controlled by using the -i
parameter, maybe setting up down to 0.8 or 0.7 for HIV. This greatly improves the collapsing of homologous regions, and results in a better looking assembly (when you look at the graph.fastg
file in Bandage.
Still, mOTUpan is not working correctly even when lowering the -i
parameter, and I'm not sure why. I have also played around with some other settings in SuperPang that might have helped, but got nothing so far.
I will discuss this with its developer.
For now, you can just manually filter the NBPs according to their coverage.
E.g. I ran SuperPang 0.9.4beta1 (available in conda) on 39 HIV genomes with SuperPang.py --fasta rawHIV/*fna -i 0.7
.
mOTUpan failed to predict a core genomes, but e.g. the header for one of the NBPs I got was:
>NODE_Sc0-0-noinfo_length_8924_cov_38.94_tag_noinfo;
So then I can see that the coverage of that NBP in my original genomes is 38.94, very close to the 39 original genomes (a couple kmers might be missing from one of the input genomes so it's not exactly 39). So I would definitely say that it is core.
Anyways I'll keep looking at this
from superpang.
I discussed this with the creator of mOTUpan, and tweaking some parameters makes it work for your data.
Just go with
SuperPang.py -f rawHIV/*fna -i 0.7 --force-overwrite -t 24 -x 70
The -i 0.7
part will collapse regions with an ANI>70%. By default it is set to 0.95
, which makes sense for bacterial species but not so much for viruses I guess.
The -x 70
tweaks mOTUpan's default completeness, which leads to a better core prediction. You can also try even higher values (e.g. 80 or 90).
from superpang.
Related Issues (15)
- IndexError: list index out of range HOT 2
- Read checkm file at the begining even if not yet needed
- multiple file paths for `--fasta` HOT 1
- IndexError: list index out of range HOT 8
- Error after identifying connected components HOT 15
- the pangeome size is larger than all size of input genomes HOT 7
- mOTUpan HOT 4
- errors when running the test-SuperPang.py HOT 4
- There was an error running homogenize.py HOT 12
- There was an error running homogenize.py HOT 13
- ValueError: invalid vertex: 3; condense-edges.py error. HOT 3
- assert len(kmers) == 2 AssertionError HOT 5
- input genome position information HOT 2
- NBP2origins.tsv and graph.NBP2origins.csv contain corrupted SourceContig names 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 superpang.