Skip to content
New issue

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

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

Already on GitHub? # to your account

redundant gff files #49

Closed
ShresthaRima opened this issue Feb 19, 2021 · 3 comments
Closed

redundant gff files #49

ShresthaRima opened this issue Feb 19, 2021 · 3 comments

Comments

@ShresthaRima
Copy link

ShresthaRima commented Feb 19, 2021

@evolarjun @kmt @slottad @nawrockie

aa.giff.txt
Hi!
I ran the code:

amrfinder -n for_amrfinder.fna -g aa.giff -O Klebsiella

and received the below error:

*** ERROR ***
Parameter --gff is redundant

The attached .txt file is aa.giff file (glimpse only). As the original .gff was from PROKKA, I changed them into the amr .giff format as per instruction at #24 and tutorial.

I could not figure out what should not be redundant in the .giff file. I see the first column with the same name.
I, being, sufficient user of linux, look forward for the solution.
thanks

@vbrover
Copy link
Contributor

vbrover commented Feb 19, 2021

A .gff file is needed to link the protein FASTA and DNA FASTA files.
Since there is no protein FASTA file, the .gff file will not be used.

@evolarjun
Copy link
Contributor

Hi Rima,

Thanks for trying AMRFinderPlus. When you only use the assembly (nucleotide) file, as Slava mentioned above, you don't need to include a GFF because AMRFinderPlus can figure out coordinates from the nucleotide file. If you just remove -g aa.giff from the command-line AMRFinderPlus should run.

However, since you already have a Prokka annotation, we would recommend you run AMRFinderPlus in the full "combined" nucleotide plus protein mode. That enables AMRFinderPlus to use HMMs to enhance sensitivity for novel AMR genes. If you want to run in combined mode you would need to include the nucleotide and GFF files (which you have used above) as well as the protein FASTA file from Prokka's output.

You would run it something like the following (I called the protein FASTA protein.fa)

amrfinder -p protein.fa -n for_amrfinder.fna -g aa.giff -O Klebsiella

Please let us know if you have further issues.

Arjun

@ShresthaRima
Copy link
Author

Thank you @evolarjun and @vbrover. It did solved the problem.

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

No branches or pull requests

3 participants