Raw reads were pre-processed by removing the adaptors and low-quality reads using BBMap. The filtered reads were normalized for depth based on kmer counts using BBNorm function. De novo transcriptomes were generated using both Trinity and velvet-oases. CD-HIT-EST was used to merge the two de novo transcriptomes and reduce the transcript redundancy to 98% similarity and generate unique genes. Transcriptome assembly completeness was evaluated with BUSCO (Benchmarking Universal Single Copy Orthologs) database. Functional annotation was done using blastp function of ncbi-blast using the nr database with evalue 1E-20 and num_alignments 3.