BASH script for SelfBLAST a genome
...-dbtype nucl -out $MYDB fi if [ $1 = "extract" ] then echo "Extracting the sequence $2 for you from $FASTAFILE -- MAKE SURE U HAVE ADDED CORRECT NAME" samtools faidx Merge...2640 days ago
2139 days ago
Installing Platypus on Ubuntu !
...ting up python3-dev (3.6.5-3ubuntu1) ... Processing triggers for libc-bin (2.27-3ubuntu1) ... ➜ htslib-1.9 git:(master) ✗ make install zsh: correct 'install' to 'INSTALL' [nyae]...2007 days ago
1541 days ago
Correct bash command to zip your vcf file
bgzip -c file.vcf > file.vcf.gz tabix -p vcf file.vcf.gz #More at http://www.htslib.org/doc/tabix.html1511 days ago
picard tools command to get some insert statistics
#picard tools to get some insert statistics to see whether our reads seem to be in the correct place #module load picard/2.0.1 java -Xmx16g -XX:PermSize=8g -jar $PICARD_HOME/picar...1325 days ago
Set up WGD environment using conda !
...kit-learn-0.24.0 scipy-1.6.0 seaborn-0.11.1 six-1.15.0 sklearn-0.0 threadpoolctl-2.1.0 tornado-6.1 wgd-1.2 (wgd) ➜ wgd git:(master) wgd zsh: correct 'wgd' to 'gd' [nyae]? n Usag...1191 days ago
Parse the NCBI taxonomy database with Perl !
...s', 'metazoa'); #exit; #my $nam = $dbh->get_taxon(-name => $sps); $nam = $dbh->get_taxon(-taxonid => 9606); if(!$nam) { print "Try again with correct scientific name\n"; exit;} p...1160 days ago
Tadpole is 250x faster than SPADes assembler !
...out2=ecc_libA_r2.fq,null mode=correct Extending contigs with rea...s with overall probability of correctness below this. minprobmain=...perform error correction. corr...-counting, and again prior to correction. Output...964 days ago
Raku script to find overlaps between two bed files !
#!/usr/bin/env raku # Check if the correct number of arguments are provided if @*ARGS.elems != 2 { say "Usage: ./compare_bed_files.raku file1.bed file2.bed"; exit 1; }...82 days ago