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
0 answers
No answers yet.
Log in to answer this question.