my @matrix = (
    [t1 t1 t2 t2 t2], <- transcripts
    [a1 a2 a1 a1 a2], <- alleles
    [intron intron UTR_CG UTR UTR], <- locations
);