Skip to content

Instantly share code, notes, and snippets.

View wyim-pgl's full-sized avatar
😜

Won Cheol Yim wyim-pgl

😜
View GitHub Profile
@wyim-pgl
wyim-pgl / maker_genome_annotation.md
Last active July 26, 2023 07:19 — forked from darencard/maker_genome_annotation.md
In-depth description of running MAKER for genome annotation.

Genome Annotation using MAKER

MAKER is a great tool for annotating a reference genome using empirical and ab initio gene predictions. GMOD, the umbrella organization that includes MAKER, has some nice tutorials online for running MAKER. However, these were quite simplified examples and it took a bit of effort to wrap my head completely around everything. Here I will describe a de novo genome annotation for Boa constrictor in detail, so that there is a record and that it is easy to use this as a guide to annotate any genome.

Software & Data

Software prerequisites:

  1. RepeatModeler and RepeatMasker with all dependencies (I used NCBI BLAST) and RepBase (version used was 20150807).
  2. MAKER MPI version 2.31.8 (though any other version 2 releases should be okay).
  3. [Augustus](http://bio
@JEEN
JEEN / ascendolify.pl
Last active December 19, 2015 12:19
#!/usr/bin/perl
my $file1 = $ARGV[0];
my $file2 = $ARGV[1];
open my $R, "<", $file1 or die ("Could not open $file!");
open my $P, "<", $file2 or die ("Could not open $file!");
my %data = ();
while (my $Rline = <$R>) {
$Rline =~ s/[\r\n]//g;
@brantfaircloth
brantfaircloth / get_protein.py
Created April 3, 2011 23:50
Get protein sequences from Genbank given a genomic accession number and a gene name
import sys
import time
from Bio import Entrez
Entrez.email = "[email protected]"
if not Entrez.email:
print "you must add your email address"
sys.exit(2)
# create an empty list we will fill with the gene names
@brantfaircloth
brantfaircloth / gist:895282
Created March 30, 2011 20:55
Substitution models in mrbayes

##GTR lset applyto=() nst=6 # GTR lset applyto=() nst=6 rates=propinv # GTR + I lset applyto=() nst=6 rates=gamma # GTR + gamma lset applyto=() nst=6 rates=invgamma # GTR + I + gamma

##SYM

lset applyto=() nst=6                           # SYM

prset applyto=() statefreqpr=fixed(equal)