Python script to download covid genome !
...w.gisaid.org/ has more (1200?), but they require you to sign up import requests import yaml seqs = yaml.load(requ...ank-sequences'] print("got %d sequences" % len(seqs)) from Bio import Entrez allseq = {} for x in...1130 days ago
Blast result parser with Perl and Bioperl
...ue, bit score, query frame, query start, query end, hit start, hit end, positives, and identical # The results are tab-deliminated and ready for import into a spreadsheet program fo...2908 days ago
Plot the density of genes in R
...if("ggplot2" %in% rownames(installed.packages())){ library(ggplot2) } else { install.packages("ggplot2") library(ggplot2) } # import a text file with gene positio...2278 days ago
Plot custom gene density with R
library(karyoploteR) pp2243 days ago
2069 days ago
Backward Elimination with Python
#Backward Elimination with p-values only: import statsmodels.formula.api as sm def backwardElimination(x...pt, SL) #Backward Elimination with p-values and Adjusted R Squared: import statsmodels.formula.api as sm...1949 days ago
Installing docker for Bioinformatics on Ubuntu !
...er export Export a container's filesystem as a tar archive history Show the history of an image images List images import Import the contents from a tarball t...1507 days ago
Commands to install conda in Ubuntu !
...urr ently reside may have restrictions on the import, possession, use, and/or re-e...aws, regulations and policies concerning the import, possession, or use, and re-e...alling: html5lib-1.0.1-py37_0 ... installing: importlib_metadata-0.8-py37_0 ... i...1505 days ago
Python script to check sequence length in multifasta file
#!/usr/bin/python from Bio import SeqIO import sys cmdargs = str(sys.argv) for seq_record in SeqIO.parse(str(sys.argv[1]), "fasta"): output_line = '%s\t%i' % \ (seq_record.id, len(seq_record)) print(output_line)1209 days ago
Print in terminal with python !
#!/usr/bin/env python import time import curses def pbar(window): height, width = window.getmaxyx() for i in range(10): window.addstr(height -1, 0, "[" + ("=" * i) + ">" + (" " * (10 - i )) + "]") window.refresh() time.sleep(0.5) curses.wrapper(pbar)972 days ago