Create genome scaffolding with Perl
...l_scaffold_", # prefix for contig names "pid" => 90, # percent ID threshold "trimlimit" => 50, # max...od likelihood that this will work ## i.e. trimLength * (1-%id) < threshold my $trimTotal = ($preTri...2321 days ago
Perl script to remove fasta sequences in multifasta file with certain length threshold
#!/usr/bin/perl use strict; use warnings; my $minlen = shift or die "Error: `minlen` parameter not provided\n"; { local $/=">"; while() { chomp...2264 days ago
2140 days ago
Install Parrot Virtual Machine !
...-g --gc ms2|gms|ms|inf set GC type --gc-dynamic-threshold=percentage maximum memory wasted by GC --gc-min-threshold=KB --gc-nurs...1550 days ago
2016 days ago
Samtools commands for bioinformatician !
...l_vars.vcf ### note: new command is "call". type "bcftools" to see current version and commands ### calling snp and indels with no frequency threshold # omit -v parameter bcftool...1615 days ago
963 days ago