Invoke thusly
$ python bogo.py
Then, enjoy!
You can edit the lines:
# Seriously, how should I format this? | |
domain: | |
kingdom: | |
subkingdom: | |
branch: | |
infrakingdom: | |
superphylum: | |
phylum: | |
subphylum: |
GTCTTCGATGTTGCCGTGGACCTGAGGGAGGGCTCTGCCACATACGGAAAGTGGTATGGGGTTGTGCTGTCTGCTGAAAATNAGAAGAAGTTCTTCATTC------------------------------- | |
-------------------------------GCTCTGCCACATACGGAAAGTGGTATGGGGTTGTGCTGTCTGCTGAAAATAAGAAGCAGTTCTTCATTCCGAAGAATTTTGCGCATGGGTTCCTTGTGTT | |
winner! 38 | |
GGAAACAGCACACGCCAACCCAAACCGAACGACGGAATGGCCTGTGTAGCTCCCCACATAGCGATGTAAGGTGCAAGAAACGAAGCAATCGATTTTACAA------------------------------------- | |
-------------------------------------TGGCCTGTGTAGCTCCCCACATAGCGATGTAAGGTGCAAGAAACGAAGCAATCGCTTTTACAAATTGTCCGAAGGTAAGGCTACTTGCCAGCCGGTCACC | |
winner! 43 | |
CCAAGAATACAAGATAGATGCGTGATCTTTTGGAATCATCAAATGAGCACGAGCTGTGAATTTTTTATCAAGACCAACGATAGATTCTGCAGAACTCCAC------------------------------------------ | |
------------------------------------------ATGAGCACGAGCTGTGAATTTTTTATCAAGACCAACGATAGATTCTGCAGAAATCCACTCTTTACTACGGCTATCATATACTGCATCACGTGCAATATTC | |
winner! 44 | |
TTAATCATGCCTGGAGCTCCTTGTTCCGCACGCGGCAGGATANGAAAACCAGTCACGTTATATTCCTTGTCGAGATCTATCACCACCTGTTGGGGTTAGG------------------------------------------- |
def function(argument1, argument2): | |
print argument1, | |
print argument2[::-1] | |
function('Austin','sivaD') |
bigblast.qsub: | |
#!/bin/sh | |
#$ -N Big_Blast | |
#$ -m e | |
#$ -e sge.err | |
#$ -o sge.out | |
#$ -pe threaded 8 | |
JOB_ID=bigblast |
#!/usr/bin/env ruby | |
# encoding: utf-8 | |
require 'rubygems' | |
require 'Bio' | |
require 'hpricot' | |
# always tell NCBI who you are | |
Bio::NCBI.default_email = "harekrishna@gmail.com" |
# List of nucleotide symbols used in FASTA files: | |
Nucleotides = { | |
'A': 'A', # Adenosine | |
'C': 'C', # Cytidine | |
'G': 'G', # Guanine | |
'T': 'T', # Thymidine | |
'U': 'U', # Uridine | |
'R': 'AG', # Purine | |
'Y': 'TC', # Pyrimidine |
# Playing code-golf | |
# Averages sequences in a FASTA file (wrapped or unwrapped) | |
a,b=0,[];File.new(ARGV.shift).each do |i|;if i[0]!=62;a+=i.length;else;b.push a;a=0;end;end;puts b.inject{|s,v|s+=v}/b.length |
##gff-version 3 | |
##sequence-region mitochondrial_genome 1 17221 | |
mitochondrial_genome GenBank mitochondrial_chromosome 1 17221 . + . ID=mitochondrial_genome;Name=mitochondrial_genome;Note=gi_189182363 | |
mitochondrial_genome GenBank misc_feature 15350 17221 . + . Name='control region' | |
mitochondrial_genome GenBank gene 2675 3635 . + . ID=ACD81883.1;Name=ND1 | |
mitochondrial_genome GenBank CDS 2675 3635 . + 0 ID=GI:189182364;Parent=ACD81883.1;Name='NADH dehydrogenase subunit 1';Note='TAA stop codon is completed by the addition of 3' A residues to the mRNA' | |
mitochondrial_genome GenBank gene 3843 4878 . + . ID=ACD81884.1;Name=ND2 | |
mitochondrial_genome GenBank CDS 3843 4878 . + 0 ID=GI:189182365;Parent=ACD81884.1;Name='NADH dehydrogenase subunit 2';Note='TAA stop codon is completed by the addition of 3' A residues to the mRNA' | |
mitochondrial_genome GenBank gene 5260 6804 . + . ID=ACD81885.1;Name=COX1 | |
mitochondrial_genome GenBank CDS 5260 6804 . + 0 ID=GI:189182366;Parent=ACD81885.1;Name='cytochrome c oxidase subunit I' |