Skip to content

Index for manual {species}_pangenome.tsv #37

@Hocnonsense

Description

@Hocnonsense

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!

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions