Create genome scaffolding with Perl
...RSION = "1.00"; =head1 NAME psl_scaffolder.pl - use self-mapped PSL file to scaffold a genome =head1 SYNOPSIS...query sequences "prefix" => "psl_scaffold_", # prefix for contig names...for (my $i = 0; $i {"prefix"}, $nextScaffoldID++); if(!exists($repl...2298 days ago
Perl script to break the contigs by 'N'
...g sequence my $seq = $seqobj->seq; my $lenseq=length $seq; $seq = uc $seq; #now all bases are in uppercase #Searches for NNNNN regions (scaffold breaks) #Hash nregions stor...2053 days ago