my ($name,$locus_id,$chrom,$transcript) = split(/\|/,$gene_header);
$gene = $transcript; #???
}elsif ($species eq "Cucumis_sativus") {
- my ($name,$gene_id,$scaffold,$unknown) = split(/\|/,$gene_header);
- print "$name\t$gene_id\t$scaffold\t$unknown\n";
- $gene = $gene_id;
+ my ($name,$main_isomer,$scaffold,$isomer) = split(/\|/,$gene_header);
+ $gene = $isomer;
}elsif ($species eq "Danio_rerio") {
my ($gene_id,$transcript) = split(/\|/,$gene_header);
$gene = $gene_id;