diff --git a/alignment/filter_genefusions.pl b/alignment/filter_genefusions.pl index 44b4453f7a598a0284327b8a9b5fc9143f3abe48..f44cf58f34f07ec90b277e1944f3e8ffd042f271 100755 --- a/alignment/filter_genefusions.pl +++ b/alignment/filter_genefusions.pl @@ -29,6 +29,8 @@ while (my $line = <OM>) { my @exonfiles = `ls */*.exons.txt`; foreach $efile (@exonfiles) { chomp($efile); + my @leftexons; + my @rightexons; my ($dir,$pair,@etc) = split(/\/|\./,$efile); open EFILE, "<$efile" or die $!; my $header = <EFILE>; diff --git a/alignment/starfusion.sh b/alignment/starfusion.sh index 73599fb76244a6a0b4ab86b1a664dd7af3122c96..d462f5ad4c7e252f7d136eac6ee7dfbb1e423469 100644 --- a/alignment/starfusion.sh +++ b/alignment/starfusion.sh @@ -60,7 +60,7 @@ else cp ${pair_id}_star_fusion/star-fusion.fusion_candidates.final.abridged ${pair_id}.starfusion.txt fi -module load singularity/2.6.1 +module load singularity/3.0.2 export PYENSEMBL_CACHE_DIR="/project/shared/bicf_workflow_ref/singularity_images" cut -f 5-8 ${pair_id}.starfusion.txt |perl -pe 's/\^|:/\t/g' | awk '{print "singularity exec /project/shared/bicf_workflow_ref/singularity_images/agfusion.simg agfusion annotate -db /project/shared/bicf_workflow_ref/singularity_images/pyensembl/GRCh38/ensembl92/agfusion.homo_sapiens.92.db -g5", $1,"-j5",$4,"-g3",$6,"-j3",$9,"-o",$1"_"$4"_"$6"_"$9}' |grep -v 'LeftGene' |sh