chaopower
/
pipeline
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

splicing 位点 exon 展示

master
chaopower 2023-12-11 13:49:59 +08:00
parent 6c1b428d45
commit 388296bffa
1 changed files with 9 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
script/filter_snpindel.pl
View File

@ -123,8 +123,16 @@ while (<IN>) {
}
$hgvs =~ /(\S+):exon(\d+):c\.(\S+)$/;
my $spl = $3;
my $exon = $2;
# splicing 位点 的前后 1-2bp 的 或者 MET的基因
if ($spl =~ /\d+[\+|\-][1|2]\D+/) {
if ($spl =~ /\d+\+[1|2]\D+/) {
my $intron = $exon;
$hgvs =~ s/exon(\d+)/exon$exon;intron$intron/;
$line[9] = join(":", ($gene, $hgvs));
}
elsif ($spl =~ /\d+\-[1|2]\D+/) {
my $intron = $exon-1;
$hgvs =~ s/exon(\d+)/intron$intron;exon$exon/;
$line[9] = join(":", ($gene, $hgvs));
}
elsif ($gene eq "MET") {