Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Index for manual {species}_pangenome.tsv #37

Open
Hocnonsense opened this issue Aug 10, 2022 · 0 comments
Open

Index for manual {species}_pangenome.tsv #37

Hocnonsense opened this issue Aug 10, 2022 · 0 comments

Comments

@Hocnonsense
Copy link

panphlan/panphlan_map.py

Lines 361 to 362 in a7846e4

for gene, (begin, end) in contig2gene[contig].items():
if begin <= position and position <= end +1:

Hello and thanks for your code.
I'me tring to apply this software to my metagenome analysis, and I've already annote all my MAGs with prodigal followed by mmseqs and eggNOG. So I'd like to generate a {species}_pangenome.tsv table without PanPhlAn_pangenome_exporter.
However, I'me not sure the style of genome start and stop in this table. For example, in gff table, a gene "starts at 3 and ends at 11" refer to a 6 bp sequence:

>contig_1
AATCGTCGTCGA
  ^       ^

>contig1|gene1
TCGTCGTCG

However, what does the number in "start and stop" of {species}_pangenome.tsv should be?, can I just add the raw number from gff file to this table?
Thanks for your advices!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant