Act as an experienced bioinformatician proficient in viral genomic assembly and annotation task. I have 3 input files, (1) 10142.1.149555.ATGTC.subset_500k.fastq.gz (2) adapters.fa (3) Ref_database folder with reference genome fa files. For the interleaving fastq.gz file split into read1 and read2, then trim the reads. After trim, I want to use the trim files to do metaviralSPAdes, which is a specific mode of SPAdes for metaviral assembly with "--metaviral" option in command. After assembly please check the contig quality using trim files. Then, builds a custom database from reference sequences, then annotates contigs and filtered read pairs using the database. make sure downloads taxonomy data and adds reference sequences from the 'Ref_database' folder to the database. use this database to annotate assembled contigs, and annotate filtered reads, put results in separate folder. Please load the following modules/tools before using them: trimmomatic/0.39 , spades/3.15 , quast/5.0 , kraken2. Use 2 threads and write code in bash.  In your entire response, the letter f should appear less than 40.
user