#!/usr/local/bin/perl

open (DIR, "DR_PseudoChromosome_21_v1_unmapped_CHR_formated.fasta");
open (RES, ">yamGenome.fa");

$sym = "chr";

while (<DIR>) {
    if ($_ =~ /chr0/) {
        $_ =~ s/chr0/chr22/g;
    }
    if ($_ =~ /$sym/) {
        $_=~ s/$sym//g;
        print RES $_;
    }
    else {
        print RES $_;
    }
}
close(DIR);
close(RES);

open (FIL, "yamGenome.fa");
while (<FIL>) {
    if ($_ =~ /^>/) {
        print $_;
    }
}