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

Calculating Gene-wise Aminoacid Diversity (aa-pi) bug #7

Open
jianshu93 opened this issue Oct 15, 2021 · 3 comments
Open

Calculating Gene-wise Aminoacid Diversity (aa-pi) bug #7

jianshu93 opened this issue Oct 15, 2021 · 3 comments

Comments

@jianshu93
Copy link

Dear Author,

when running the following:

perl pogenom.pl --vcf_file pico127_pico127.002.vcf --out pico127_pico127.002_pi_fst_gene --gff_file pico127_pico127.002.gff --genome_size 2660050 --genetic_code_file standard_genetic_code.txt

I have the following error:

Use of uninitialized value $mod_contig_seq in substr at pogenom.pl line 924.
substr outside of string at pogenom.pl line 924.

Any idea? vcf file is generated by the Input_POGENOM.sh and gff is generated by prodigal

Thanks,

Jianshu

@jianshu93
Copy link
Author

If I want to run multiple samples should I change the dataset name to "dataset1,dataset2,dataset3"?

How should I do it

Thanks

Jianshu

@lfdelzam
Copy link
Collaborator

lfdelzam commented Oct 21, 2021

question 2) R/ INPUT POGENOM only runs one dataset at a time. It is not possible to use ""dataset1,dataset2,dataset3" as you suggested.

@JFsanchezherrero
Copy link

Hi there,

I guess I found a solution for question 1.

Basically, gff file might contain fasta sequence at the end but pogenom.pl does not read properly neither fasta sequence nor genome size from GFF and basically although all the information is included within GFF you need to include --genome_size too and in this case, the FASTA file.

So, the solution for the error is to include the FASTA file too.

perl pogenom.pl --vcf_file pico127_pico127.002.vcf --out pico127_pico127.002_pi_fst_gene --gff_file pico127_pico127.002.gff --genome_size 2660050 --genetic_code_file standard_genetic_code.txt --fasta_file XXXX.fna

It works fine, it produces no error and Gene-wise Aminoacid Diversity (aa-pi) is calculated.

Best regards
Jose

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

3 participants