Skip to content

Index for manual {species}_pangenome.tsv #37

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!

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

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