Hi, I have a problem annotating genome of my rabies samples. Although I have a file with gene specific proteins in .faa format, at the time of annotation it marks me that some are either smaller or larger. The code I am using is the following, I don't know if I am making a mistake.
Thank you very much in advance.
for assembly in /home/secuenciacion_cenasa/Analisis_corridas/SPAdes_viral/*.fa; do
ID=$(basename ${assembly} | cut -d '-' -f '1')
ename=$(basename ${assembly} | cut -d '.' -f '1')
Hi, I have a problem annotating genome of my rabies samples. Although I have a file with gene specific proteins in .faa format, at the time of annotation it marks me that some are either smaller or larger. The code I am using is the following, I don't know if I am making a mistake. Thank you very much in advance.
for assembly in /home/secuenciacion_cenasa/Analisis_corridas/SPAdes_viral/*.fa; do ID=$(basename ${assembly} | cut -d '-' -f '1') ename=$(basename ${assembly} | cut -d '.' -f '1')
prokka ${assembly} \ --proteins /home/secuenciacion_cenasa/Programas_Bioinformaticos/prokka/db/genus/Rabia_prokka_db/tmp/Rabies_proteins.faa --outdir /home/secuenciacion_cenasa/Analisis_corridas/Resultados_allvirus/prokka/${ID}${gene}_anotado \ --evalue 1e-9 \ --coverage 100 \ --rawproduct