I am analysing soil samples for two different marker genes (bacteria and fungi) and my sequencing service provider has some unusual technical procedures: Instead of running bacteria and fungi samples separately, they combined - for each sample - fungi and bacteria first-stage PCR products and then used both together for Illumina library preparation.
In a nut-shell, I ended up with fasta files that contain two different marker genes. Has anybody tips how I can deal with that during my taxonomic classification (using two data bases, SILVA and UNITE)?
Would it work to first extract reference sequences for both silva and unite for respective primers and then merge reference sequences (using merge-seqs?) and reference taxonomy files (any idea how?) and then train the classifer on the joint files?
Thanks for any help!!