Statistic of pre-filter data. A total of 24 GB of raw data were generated and used for assembly with total sequence depth of 214.29 and total physical depth of 3203.63. Table S2. Statistic post filtering. A total of 18.76 GB data were generated after filtering process with total sequence depth of 167.5 and total physical depth of 2290.07. Table S3. Transposable elements (TEs) content in the Assembled B. pahangi Genome. The highest TEs predicted at the DNA level was found by RepBase (0.75 %). Table S4. General statistics of predicted protein-coding genes. Three de novo based methods used: AUGUSTUS, SNAP and GLIMMERHMM. Table S5. Statistics of function annotation. The method for functional annotation is divided into two types that are automated and manual. The automated method, TrEMBL database shows the highest number of annotation among all methods used with percentage of 96.64 % while InterPro database produces the most number of annotations with percentage of 71.32 % among manually curated databases. Table S6. B. pahangi genes mapped to Swissprot via blast. Table S7. B. pahangi genes with intrepro annotations. Table S8. B. pahangi genes with Gene ontology annotations. Table S9. B. pahangi genes mapped to KEGG via blast. Table S10. List of 569 B. pahangi unique genes. Table S11. The 403 Wolbachia genes in B. pahangi unique genes. Table S12. The 26 B. Pahangi unique proteins with their respective KEGG pathway class. Table S13. The 803 Wolbachia genes in B. pahangi genome. Table S14. General statistics of repeats in genome. From the table, TRF shows the biggest repeat size of 2.9 Mbp, which represent 3.34 % of the genome size. (XLSX 2292 kb)

Statistic of pre-filter data. A total of 24 GB of raw data were generated and used for assembly with total sequence depth of 214.29 and total physical depth of 3203.63. Table S2. Statistic post filtering. A total of 18.76 GB data were generated after filtering process with total sequence depth of 167.5 and total physical depth of 2290.07. Table S3. Transposable elements (TEs) content in the Assembled B. pahangi Genome. The highest TEs predicted at the DNA level was found by RepBase (0.75 %). Table S4. General statistics of predicted protein-coding genes. Three de novo based methods used: AUGUSTUS, SNAP and GLIMMERHMM. Table S5. Statistics of function annotation. The method for functional annotation is divided into two types that are automated and manual. The automated method, TrEMBL database shows the highest number of annotation among all methods used with percentage of 96.64 % while InterPro database produces the most number of annotations with percentage of 71.32 % among manually curated databases. Table S6. B. pahangi genes mapped to Swissprot via blast. Table S7. B. pahangi genes with intrepro annotations. Table S8. B. pahangi genes with Gene ontology annotations. Table S9. B. pahangi genes mapped to KEGG via blast. Table S10. List of 569 B. pahangi unique genes. Table S11. The 403 Wolbachia genes in B. pahangi unique genes. Table S12. The 26 B. Pahangi unique proteins with their respective KEGG pathway class. Table S13. The 803 Wolbachia genes in B. pahangi genome. Table S14. General statistics of repeats in genome. From the table, TRF shows the biggest repeat size of 2.9 Mbp, which represent 3.34 % of the genome size. (XLSX 2292 kb)