Insert the sequence at desire location in multi-fasta file with Perl
#!/usr/bin/perl use warnings; use strict; use Bio::SeqIO; use Bio::Seq; use File::Copy; #ARGV[0] should be in following format --- Keep the coordinate sorte...2313 days ago
Perl script to read multi fasta sequence one by one
#!/usr/bin/env perl use strict; use warnings; #USAGE #perl rohanRun.pl seq.fa my $outfile='tmp.fa'; my $fastaSeq_ref = readfasta ("$ARGV[0]"); my %fastaSeq = %$fastaSeq_r...2221 days ago
2127 days ago
2026 days ago
Installing Platypus on Ubuntu !
...c_LIBS@/-lz -lm -lbz2 -llzma/;s#@[^-][^@]*@##g' htslib.pc.in > htslib.pc.tmp sed -e 's#@-includedir@#/usr...ude#g;s#@-libdir@#/usr/local/lib#g;s#@-PACKAGE_VERSION@#1.9#g' htslib.pc.tmp > /usr/local/lib/pkgconfig/ht...2019 days ago
Bash commandline to install Anaconda !
#The line begins with $ are the commands $ mkdir tmp $ cd tmp/ $ curl -O https://repo.anaconda.com/archive/Anaconda3-2019.03-Linux-x86_64.sh % Total % Received % Xferd A...1555 days ago
Bash command to install Miniconda !
(base) [wsu29@bladeamd-2 tmp]$ wget https://repo.anaconda.com/miniconda/Miniconda3-latest-Linux-x86_64...-Linux-x86_64.sh' saved [71785000/71785000] (base) [wsu29@bladeamd-2 tmp]$ sh Miniconda3-latest-Linux-...1555 days ago
Bash command to install bwa, samtools, picard !
(base) [wsu29@bladeamd-2 tmp]$ conda install bwa Collecting package metadata...transaction: done (my_GATK) [wsu29@bladeamd-2 tmp]$ conda install samtools Col...transaction: done (my_GATK) [wsu29@bladeamd-2 tmp]$ conda install picard Colle...1555 days ago
Bash command to install GATK, Bedtools and SnpEff !
(my_GATK) [wsu29@bladeamd-2 tmp]$ conda install gatk Collecting package metada...ransaction: done (my_GATK) [wsu29@bladeamd-2 tmp]$ conda install bedtools Col...transaction: done (my_GATK) [wsu29@bladeamd-2 tmp]$ conda install SnpEff Colle...1555 days ago
Bash script to get intergenic region from genome files !
...ic region, we will require the size of the chromosomes. wget http://xxx.chrom.sizes cat xxx.chrom.sizes | sed 's/^chr//' | sed 's/Cp/Pt/' > tmp mv tmp xxx.chrom.sizes gunzip -c...1365 days ago