Summary of the bioinformatics workflow. After de novo assembly of sequencing data, individual reads were mapped to the obtained assembly i.e. transcript abundance data were calculated per individual. Subsequently expression difference analysis and BLAST searches against Uniprot-Swissprot and NCBI non-redundant database were preformed. BLAST searches were also run against both chicken and zebra finch coding DNA and genome sequences (retrieved from Ensembl).