09-Qsub-Pair_Alignment_minimap2.pl

#!/usr/bin/perl
#Author: huilong du
#Note: align the corrected pacbios and non-scaffolded contigs selfly
use warnings;
use strict;

my $infile=shift;
my $script=shift;
my $queue=shift;
my $genome=shift;
open IN,"<$infile" or die $!;


my $count=1;
my @part=();
while(<IN>){
	chomp;
	my $line=$_;
	push(@part,$line);
}
close(IN);
for(my $i=0;$i<@part;$i++){
    for(my $j=$i;$j<@part;$j++){
	open OUT,">$count.pbs" or die $!;
	if($count<350){
		print OUT "#BSUB -J $genome-Pair-$i-$j
#BSUB -o $count.out
#BSUB -n 2
#BSUB -q $queue
";
        }
        elsif($count>=350){
                print OUT "#BSUB -J $genome-Pair-$i-$j
#BSUB -o $count.out
#BSUB -n 2
#BSUB -q $queue
";
        }
	if($i==$j){
           print OUT " 
minimap2 -ax asm5 -t 4 $part[$i].mmi $part[$j] >./03-Pacbio-SelfAlignment/Part_SelfAlignment_$i-$j.sam
perl $script/sam2blasr.pl ./03-Pacbio-SelfAlignment/Part_SelfAlignment_$i-$j.sam ./03-Pacbio-SelfAlignment/Part_SelfAlignment_$i-$j.txt
rm -f ./03-Pacbio-SelfAlignment/Part_SelfAlignment_$i-$j.sam
";
	}
	else{
	   print OUT "
minimap2 -ax asm5 -t 4 $part[$i].mmi $part[$j] >./03-Pacbio-SelfAlignment/Part_SelfAlignment_$i-$j.sam
perl $script/sam2blasr.pl ./03-Pacbio-SelfAlignment/Part_SelfAlignment_$i-$j.sam ./03-Pacbio-SelfAlignment/Part_SelfAlignment_$i-$j.txt
rm -f ./03-Pacbio-SelfAlignment/Part_SelfAlignment_$i-$j.sam
";
	}
        close(OUT);
        system("qsub -cwd -V -l vf=10g,p=2 -o $count.out -q $queue -e log/ -N $genome-Pair-$i-$j < $count.pbs");
        $count++;
    }
}