Fig. 2

Rarefaction curves of observed virus family and virus distribution pattern in the 215 birds’ cloacal swab pools based on birds’ habitat. A In the rarefaction curves the horizontal ordinate represents the number of viral sequence reads obtained from the NGS data (where the “break” axis is used because some libraries contain too large a number of viral reads), the longitudinal axis represents the number of virus families observed in the virome of birds’ cloaca swab samples. B The horizontal ordinate represents different virus genomes which are further divided into 24 different virus families or groups, while the longitudinal axis represents the 215 libraries that are arranged based on birds’ habitat (breeding or wild). Heatmap representing the read number (in exponential form) in the mapping analysis using the 707 genomes against NGS data of the 215 libraries (see color legend)