Expression atlas

The merged dataset of the short-read and long-read structural annotations of six tissues was filtered for UniRef50 hits and an expression atlas created based on support of pipeline/tissue.The TAMA tool tama_mergy.py creates a gene report containing each features's source (pipeline, tissue). This gene report was parsed with tama_merge_report_parser.pl, filtered for UniRef50 hits and loaded in the R function UpSetR to visualise tissue-specific evidence for expressed genes.Parse gene report, prepare for UpSetR: ${tama_merge_report_parser.pl supermerge_hisat2_stringtie_and_thresholds_ends_gene_report.txt gene_report.csv}Get genes with UniRef50 hits in final annotation: ${awk 'BEGIN {FS="[\t;]"} $6 ~ /UniRef50/ {print $4}' supermerge_hisat2_stringtie_and_thresholds_ends.bed | sed 's/$/,/' | sort -u > gene_support_uniprot.txt}Filter UpSetR list for UniRef50-supported entries: ${grep -Ff gene_support_uniprot.txt gene_report.csv > gene_report_supported_by_uniprot.csv}