Create genome scaffolding with Perl
...$b1 eq " ") || ($b2 eq " ")){ ## equal bases, or absent bases, so consensus is easy return($b1); } # if different, convert to upper case to simplify lookup my $bc...2302 days ago
Install Parrot Virtual Machine !
...ystem to figure out how to build Parrot. The process is completely automated, unless you passed in the `--ask' flag on the command line, in which case I'll prompt you for a few pi...1531 days ago
Generates a genome coverage plot with R
library(CoverageView) ##draw a coverage plot for a test case BAM file #get a BAM test file treatBAMfile2049 days ago
Perl script to delete the adjacent repeats !
...ring; #For more helps ... check out this #This bit of code will do the job: $string =~ s/(.)\1/$1/g #Alternatively, if you want to ignore case when comparing, you could use...1544 days ago
Parse the NCBI taxonomy database with Perl !
...nt $_->scientific_name; print "\n";} my $decision =recheck_class (1974465, "metazoa"); print $decision; #It return 1 of matches otherwise 0 #case-insesitive #Spelling mistake...1166 days ago
Oneliner to convert lower-case to sequence masked with Ns
perl -pe '/^[^>]/ and $_=~ s/[a-z]/N/g' genomic.fna > genomic.N-masked.fna awk '{if(/^[^>]/)gsub(/[a-z]/,"N");print $0}' genomic.fna > genomic.N-masked.fna944 days ago
944 days ago
943 days ago
Get the Linux system information !
...e Space Utilization" \ 2>&1 1>&3) exit_status=$? exec 3>&- case $exit_status in $DIALOG_...ar echo "Program aborted." >&2 exit 1 ;; esac case $selection in 1 )...942 days ago
Downloading mmseqs databases !
...download and setup the Swiss-Prot at the output path outpath/swissprot: mmseqs databases UniProtKB/Swiss-Prot outpath/swissprot tmp #In this case, since Swiss-Prot has a value...936 days ago