Repositories (1)
gramene-ensembl
4 snippets stored for this repo
print stderr "$name\t";
# cycle through the regexes and see where the bastard falls
my $matches = 0;
foreach my $analysis (keys %regexes){
foreach my $regex (keys %{$regexes{$analysis}}){
#f/ complete twat handling!
#y are there any actual exons or CDS - i.e. there should be NO 3o features w/o
#b junk
# TWAT YOU CAN'T PUT THIS CODE AFTER THE INCREMENTER - TAHT MUST BE LAST - just printing for diff
print STDOUT qq{\n-----------------\n| gene strand = $gene_strand},
#y to retive the translation directly
#print $trans[0]->translation().$n;
#/ twat the uc was applied in situ and not as an lvalue?!?
print '================================================================'.$n;
print 'Checking gene '.$gene.$n.$n;
print 'seq from gff file (has length: '.$GffPepLen.'): '.$n.$GffPeptideSeq.$n.$n;