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 = "[email protected]" |
| # 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' |