Aligned and sorted bam files were loaded into IGV [50] to perform visual inspection of target variants. Samples carrying unexpected variants or missing expected variants were labeled as potential swapped samples. In addition, variants were called from RNA-seq following BWA/GATK pipeline [51, 52]. The identity of the samples was later verified by performing IBD analysis against genomic typing from genome-wide association study chipsets.
