1566 days ago
Install Ragout genome assembler
$ conda install -c bioconda ragout Collecting package metadata (repodata.json): done Solving environment: done ## Package Plan ## environment location: /hom...1521 days ago
1499 days ago
get GC across the entire CDS !
#look at GC across the entire CDS. gffread -x - -g | \ seqtk comp - | \ awk -v OFS="\t" '{ print $1, "0", $2, ($4 + $5) / $2 }'1450 days ago
Perl One-Liner to print only non-uppercase letters
#Go through file and only print words that do not have any uppercase letters. perl -ne 'print unless m/[A-Z]/' dna.fa > dnaOnlyLowercase.fa #To lowercase everything perl -pne 'tr/[A-Z]/[a-z]/' dnaUpperCase.fa >dnawithoutuppercase.fa;1441 days ago
Script to extract the cluster detail !
$ lsb_release -a No LSB modules are available. Distributor ID: Ubuntu Description: Ubuntu 18.04.1 LTS Release: 18.04 Codename: bionic $ cat /proc/cpuinfo | grep -i 'model name' | head -n 1 model name : Intel(R) Xeon(R) CPU E5-2620 v4 @ 2.10GHz1435 days ago
Bash script to get exon fragments from genome files !
#Exons are already defined in the GTF file, so we simply need to print lines that are marked exonic. gunzip -c genome_file.gtf.gz | awk 'BEGIN{OFS="\t";} $3=="exon" {print $1,$4-1,$5}' | bedtools sort | bedtools merge -i - | gzip > my_exon.bed.gz1423 days ago
Bash script to extract intronic fragments !
...simply need the gene and exonic coordinates; #by subtracting the exonic regions from the genic region, we have the intronic region. gunzip -c genome_file.gtf.gz | awk 'BEGIN{O...1423 days ago
Bash script to get intergenic region from genome files !
...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 genome_file.gtf.gz | awk 'BEGIN{O...1423 days ago
Print the screenshot of read coverage remotely !
#Chrom_2:9700000-9760000 ASCIIGenome -ni --fasta A_vaga.NDPD.fasta -r Chrom_2:9700000-9760000 -x "save >> test.pdf" A_vaga.NDPD.fasta.out.sam.bam.sorted.bam1340 days ago