Skip to content

Instantly share code, notes, and snippets.

View ravila4's full-sized avatar

Ricardo Avila ravila4

View GitHub Profile
#!/usr/bin/env python
import pandas as pd
import click
from Bio.Seq import Seq
from Bio.SeqRecord import SeqRecord
from Bio import SeqIO
@click.command()
#!/bin/bash
TYPE=${TYPE:-prot}
[[ ! -z ${1} ]] && INFILE=${1} || exit 1
shift
makeblastdb -in ${INFILE} -dbtype ${TYPE} -parse_seqids ${@} -blastdb_version 5
blastp -db fasta.fa -query database.fa \
-outfmt "6 std stitle qcovs" -num_threads 10 -out out.blast
@ravila4
ravila4 / parallel.py
Last active September 11, 2019 18:40
Functions for parallelizing things
# Functions for parallelizing things
def init_spark(nproc=-1, appname="sparksession"):
"""Function to start a Spark executor."""
from pyspark.sql import SparkSession
if nproc == -1:
# Use all CUPs
spark = SparkSession.builder.master(
"local[*]").appName(appname).getOrCreate()
else:
@ravila4
ravila4 / .tmux.conf
Created August 5, 2019 18:26
Tmux configuration
# $Id: vim-keys.conf,v 1.2 2010-09-18 09:36:15 nicm Exp $
#
# vim-keys.conf, v1.2 2010/09/12
#
# By Daniel Thau. Public domain.
#
# This configuration file binds many vi- and vim-like bindings to the
# appropriate tmux key bindings. Note that for many key bindings there is no
# tmux analogue. This is intended for tmux 1.3, which handles pane selection
# differently from the previous versions
@ravila4
ravila4 / utility.py
Created March 8, 2019 04:05
Utility script for generating TPATF and TPAPF fingerprints from SMILES strings
# Utility script for feature generation
# Md Mahmudulla Hassan
# The University of Texas at El Paso
# Last Modified: 12/19/2018
import os
from rdkit import Chem
from rdkit.Chem import AllChem
import tempfile
import shutil
@ravila4
ravila4 / parse_drugbank_xml.py
Created March 8, 2019 04:03
Python script for parsing an xml database dump from DrugBank for extracting Log P values
import xmltodict
import pandas as pd
with open("full_database.xml") as db:
doc = xmltodict.parse(db.read())
values = []
for item in doc['drugbank']['drug']:
logp = None
try:
@ravila4
ravila4 / Pharos_MySQL_connect.ipynb
Last active May 11, 2022 23:49
Connecting to the Pharos MySQL database with Python
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@ravila4
ravila4 / BioservicesUniProt.ipynb
Created July 7, 2018 17:37
Getting UniProt Annotations
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@ravila4
ravila4 / nglview.ipynb
Created April 29, 2018 22:25
NGLView - Basics
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.