Fig. 2

Analyses of heterozygosity of MAC genome (A), intron structure (B) in Pseudourostyla cristata and nine representative ciliates. A Log plot of a Kmer spectral genome composition analysis of P. cristata. Heterozygosity was estimated by jellyfish and GenomeScope2, based on 21-mers in Illumina sequence reads of the MAC genome. The len = inferred haploid genome length, uniq = percentage non-repetitive sequence, aa = Homozygous, ab = Heterozygous, kcov = mean kmer coverage for heterozygous bases, err = error rate of the reads, dup = average rate of read duplications, k = Kmer, p = ploid. The observed 21-mers frequency distribution is depicted in purple. The black lines represent the modeled distribution of 21-mers in the full genome. The orange lines represent the modeled distribution of the unique fraction of the genome. We find ~ 104× and ~ 208× coverage for heterozygous and homozygous peaks in our dataset, respectively. B Motif sequences of intron with most abundant size category in ten representative ciliates. The purple shadow and polyline show the percentage of base A at each position in the second half of intron sequences. Asterisks represent the conserved A nucleotide (The percentage of base A at the corresponding position is greater than 50%) which most likely represents a branch site