Fig. 4

Normalized abundance of eukaryotic viral species. The heatmap shows the normalized reads counts by metagenomeSeq on log2 scale. The hierarchical clustering is based on the Euclidean distance matrix calculated from the normalized reads count. The viral species names shown in the heatmap are from the taxonomic annotation by DIAMOND and KronaTools. For each of the contigs assigned to a particular species, the ORF with the highest BLASTx identity to a reference sequence was taken, and the average identity of these different ORFs is shown in the shaded blue boxes. The red-shaded viruses were selected for qRT-PCR analysis and the names of novel viruses are shown between brackets. The samples marked with red dots are pools containing five mosquitoes and the one with star is the special sample Ab-AAF-1-3