Skip to content

Instantly share code, notes, and snippets.

View ChemicalJames's full-sized avatar

ChemicalJames

View GitHub Profile
@ChemicalJames
ChemicalJames / barcode.pl
Created December 5, 2012 05:46
Project Scripts
#!/bin/perl -w
#sort sequences into two files according to 5' barcode
use strict;
use warnings;
use Bio::SeqIO;
my $file = 'trimmed_seq.fa';
my $in = Bio::SeqIO->new(-format => 'Fasta',
#!/bin/perl -w
use strict;
use warnings;
use Bio::SeqIO;
my $seqio = Bio::SeqIO->new(-format=>'fasta',-file=>'ath_sRNA1.fa');
my $adapter='CTGTAGGCACCATCAAT';
open OUT,'>trimmed_seq.fa' or die "Cannot write to a file\n";
@ChemicalJames
ChemicalJames / Bio.DB.Fasta
Created November 1, 2012 04:52
GEN220 HW5 BioPerl
#!/usr/bin/perl -w
use strict;
use warnings;
use Bio::DB::Fasta;
use Bio::Seq;
use Bio::SeqIO;
my $file = 'sacharomyces_cerevisiae_S288C.fa';
my $db = Bio::DB::Fasta->new($file);
@ChemicalJames
ChemicalJames / Temperature.pm
Created October 24, 2012 00:47
Temperature Module
package Temperature;
use strict; use warnings;
sub f2c {
my $f = shift;
my $c = (5/9) * ($f - 32);
return $c;
}
sub c2f {
my $c = shift;
@ChemicalJames
ChemicalJames / print_line_num.pl
Created October 23, 2012 21:01
GEN220 HW4 - error debugging
#!/usr/bin/perl -w
use strict;
foreach (0..8)
{
my $line = <>;
print "$_ : $line";
}
#!bin/perl
use strict;
# number of 3'UTR
open(IN, "grep -c 'three_prime' ath_chr1.gff | ") || die $!;
while(my $count = <IN>) {
open(my $outfh => "UTR.txt") || die $!;
print $outfh "$count\n";
}
@ChemicalJames
ChemicalJames / protein_seqs
Created October 10, 2012 07:01
GEN220 Wk2 protein_seqs
#!/bin/perl
use strict;
my %seq;
$seq{'AAC35278'} = 'LLIAITYYNEDKVLTARTLHGVMQNPAWQKIVVCLVFDGIDPVLATIGV-VMKKDVDGKE';
$seq{'AnCSMA'} = 'AMCLVTCYSEGEEGIRTTLDSIALTPN-SHKSIVVICDGIIKVLRMMRD-TGSKRHNMAK';
$seq{'AfCHSF'} = 'ALCLVTCYSEGEEGIRTTLDSIAMTPN-SHKTIIVICDGIIKVLRMMRD-TGSKRHNMAK';
$seq{'AAF19527'} = 'AILLVTAYSEGELGIRTTLDSIATTPN-SHKTILVICDGIIKVLGMMKD-RGSKRHNMAK';
$seq{'P30573-1'} = 'TINLVTCYSEDEEGIRITLDSIATTPN-SHKLILVICDGIIKVLDMMSDAQGSKRHNMAK';
@ChemicalJames
ChemicalJames / mRNA_seq
Created October 10, 2012 06:49
GEN220 Wk2 mRNA_seq
#!/usr/bin/perl
use warnings;
use strict;
my $seq ="AGACAAGTCGGACGTTTCATCTGAGGGTTCTTCTGCCTCCGCACTTGGTGCACATCAGACAAGGCAATCA
TGGGGGACGCTCAGATGGCAGAGTTTGGAGCAGCAGCTTCTTACCTGCGAAAGTCAGATCGAGAGCGTCT
GGAAGCACAAACCCGTCCCTTTGATATGAAAAAGGAGTGTTTTGTGCCTGATCCAGATGAAGAGTATGTA
AAAGCTTCAATCGTCAGTCGTGAAGGTGACAAAGTCACTGTACAGACTGAGAAAAGAAAGACTGTAACTG
TAAAGGAAGCTGACATTCACCCCCAGAACCCTCCAAAGTTTGATAAAATTGAAGACATGGCAATGTTCAC
CTTCCTTCATGAGCCAGCCGTGCTGTTCAACCTCAAAGAGCGCTATGCAGCATGGATGATCTATACCTAC
TCAGGACTGTTTTGTGTCACTGTCAACCCCTACAAGTGGCTGCCGGTGTACAATCAGGAGGTGGTTGTAG