hisat2 align

begining to analysis raw data, you should be clean the fastq data. Here we use the hisat2 software to do this project.

#usr/bin/perl
use strict;
open IN,"sqf.sh";

my %fg;
my @gg;
while (<IN>) {
print $_;
}

#usr/bin/perl
use strict;
open I,"Homo_sapiens.GRCh38.81.gtf";
open IO,">Homo_sapiens.gtf";

while (<I>) {
    chomp;
    my @ss = split/\t/;
    next if($ss[0] =~ m/CHR/);

    my ($a)= (split/ /,$ss[-1])[1];
    my $al = join"\t",@ss[0..(@ss-2)];
     $a =~ s/\"//g;
     $a =~ s/\;//g;
    print IO $al."\t".$a."\n";
    

    # body...
}

<img src="C:\Users\Lenovo\Desktop\iI\微信图片_20200814122753.jpg" alt="微信图片_20200814122753" style="zoom: 25%;" />

<img src="C:\Users\Lenovo\Desktop\iI\未标题-1.png" style="zoom:10%;" />