This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Missing taxonomic information from Trinotate annotation report

Hello there,

I have loaded my BLASTX and CD-HIT-EST information into a Trinotate.sqlite file, and then tried to run a trinotate annotation report. However, I am getting gene names but no taxonomic information after that (i.e., kingdom) TRINITY_DN88414_c0_g1 TRINITY_DN88414_c0_g1_i2_batch2_hall sp|Q9SYM0|VTE6_ARATH^sp|Q9SYM0|VTE6_ARATH^Q:127-606,H:67-226^75.625%ID^E:4.3e-90^.^.

Before, I have run this and have not had the problem, so I am not sure why this is happening and what can be done to fix this. I am needing to filter out my genes by fungi, and I cannot do this if I don't have the kingdom info to screen for fungal genes. Below is my job script for the blastx database:

#!/bin/bash

#SBATCH --nodes=1

#SBATCH --ntasks-per-node=96

#SBATCH --mem=1860GB

#SBATCH --time=168:00:00

#SBATCH --job-name=blastx_NYCHA_hall

#SBATCH --account=PAS1182

#SBATCH --mail-type=ALL

module load blast-plus/2.16.0

export BLASTDB=$HOME/local/src/database

blastx -query /fs/ess/PAS1182/Olivia/NYCHA/cd_hit_est_2/combined_cd/CD_HIT_EST_hall_tag.fasta -db uniprot_sprot -num_threads 78 -max_target_seqs 1 -outfmt 6 -evalue 1e-3 > /fs/ess/PAS1182/Olivia/jobs/annotations/batch2/hall/blastx.outfmt6

and trinotate annotation report: #!/bin/bash

#SBATCH --nodes=1

#SBATCH --ntasks-per-node=48

#SBATCH --time=168:00:00

#SBATCH --job-name=Trinotate_Annotation_Report_NYCHA_1

#SBATCH --account=PAS1182

#SBATCH --mail-type=ALL

export TRINOTATE_HOME=$HOME/local/src/Trinotate-Trinotate-v3.2.2

cd /fs/ess/PAS1182/Olivia/jobs/annotations/batch2/hall

$TRINOTATE_HOME/Trinotate Trinotate.sqlite report -E 0.00001 > trinotate_annotation_report_batch2_hall_NYCHA_.xlsx

Any and all help would be greatly appreciated.

Thanks,

Olivia

blastx trinotate rna-seq

0 answers

No answers yet.

Log in to answer this question.