-
Notifications
You must be signed in to change notification settings - Fork 6
Open
Description
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
Labels
No labels