Generate simulated polyploid genome !
#Generate 3% divergence msbar -point 4 -count 16558 toy.fasta > toyheterozygous3percent.fasta #Cat both files cat toy.fasta toymutated3percent.fasta > toyheterozygo...2107 days ago
Installing pb-assembly on Linux !
[jnarayan@hmem00 ~]$ module avail -------------------------------------------------------- /usr/share/Modules/modulefiles ---------------------------------------------...2089 days ago
2081 days ago
Perl script to delete the adjacent repeats !
/usr/bin/perl #Mostly the interview question for bioinfomatician ! #Write a code to delete the adjacent repeated character .... $string='ATTTTTTGGC'; # This should...1605 days ago
Exit Bash Script When Any Command Fails
This can actually be done with a single line using the set builtin command with the -e option. #1-----------------------------------------------------------------------...1575 days ago
Onliner to split the multifasta to singlefasta files !
#Split the multifasta to singlefasta # Multi fasta #Single fasta awk '$0 ~ "^>" { match($1, /^>([^:]+)/, id); filename=id[1]} {print >> filename".fa"}' sequence.fasta1455 days ago
Bash script to handle Multifasta files
#Convert all lowercase residues to uppercase in a FASTA sequence file $ awk 'BEGIN{FS=" "}{if(!/>/){print toupper($0)}else{print $1}}' input.fasta > output.fasta #Re...1412 days ago
1005 days ago
Command line to download blast database / protein
#download all available nr - protein database as a single file #Database location - NCBI where all databases are available ftp://ftp.ncbi.nlm.nih.gov/blast/db/ https...998 days ago
Onliner to convert multi line fasta to single line fasta !
#Oneliner to convert awk '/^>/ {printf("\n%s\n",$0);next; } { printf("%s",$0);} END {printf("\n");}' < file.fa > fileres.fa #Then delete the first empty line tail -n +2 fileres.fa > fileout.fa983 days ago