Mercurial > repos > lxue > ageseq
comparison AGEseq_web.pl @ 7:449c8cf8fa3f draft default tip
Handling of large data set.
author | lxue |
---|---|
date | Tue, 09 Jun 2015 12:44:58 -0400 |
parents | a9c5e846dd76 |
children |
comparison
equal
deleted
inserted
replaced
6:8d5914e21ec8 | 7:449c8cf8fa3f |
---|---|
525 | 525 |
526 my $total_num = $hash_out_num{$fas_file_in}{'total'} ; | 526 my $total_num = $hash_out_num{$fas_file_in}{'total'} ; |
527 my $total_non_redun = $hash_out_num{$fas_file_in}{'sub'} ; | 527 my $total_non_redun = $hash_out_num{$fas_file_in}{'sub'} ; |
528 | 528 |
529 for my $ref_name (sort keys %{$hash_out{$fas_file_in }}){ | 529 for my $ref_name (sort keys %{$hash_out{$fas_file_in }}){ |
530 if(not exists $hash_out{$fas_file_in}{$ref_name}{"data"}){ | |
531 # print "No data for $fas_file_in $ref_name \n"; | |
532 next; | |
533 } | |
534 | |
530 my @data = @{$hash_out{$fas_file_in}{$ref_name}{"data"}}; | 535 my @data = @{$hash_out{$fas_file_in}{$ref_name}{"data"}}; |
531 | 536 |
532 for (@data){ | 537 for (@data){ |
533 $print_tracking++; | 538 $print_tracking++; |
534 if($print_tracking==1){ | 539 if($print_tracking==1){ |
556 print REPORT "Sum:INPUT\t","Target\t","AlignedTarget\t","AlignedRead\t","Total Hits\t","Sub Hits\t","Indel or WT Hits\t","Indel or WT rate %\t","Pattern\n"; | 561 print REPORT "Sum:INPUT\t","Target\t","AlignedTarget\t","AlignedRead\t","Total Hits\t","Sub Hits\t","Indel or WT Hits\t","Indel or WT rate %\t","Pattern\n"; |
557 | 562 |
558 | 563 |
559 | 564 |
560 for my $ref_name (sort keys %{$hash_out{$fas_file_in }}){ | 565 for my $ref_name (sort keys %{$hash_out{$fas_file_in }}){ |
566 | |
567 if(not exists $hash_out{$fas_file_in}{$ref_name}{"data"}){ | |
568 # print "No data for $fas_file_in $ref_name \n"; | |
569 next; | |
570 } | |
561 my @data = @{$hash_out{$fas_file_in}{$ref_name}{"data"}}; | 571 my @data = @{$hash_out{$fas_file_in}{$ref_name}{"data"}}; |
562 | 572 |
563 my $wt_pair = ''."\t".''; | 573 my $wt_pair = ''."\t".''; |
564 my $indel_pair = ''."\t".''; | 574 my $indel_pair = ''."\t".''; |
565 my $indel_out = ''; | 575 my $indel_out = ''; |