Biological Sequence handling with Perl !
...(__PACKAGE__ . '::Subsequence')->new($self,$self->length,1); } sub translate { my $self = shift; my...s,$frame); foreach (-3..-1,1..3) { my $translation = $self->translate($_); while ($translation...2181 days ago
Perl script to convert GFF 2 FASTA !
...tream3000.fasta" ); ###### Output type description ###### # cds - translated sequence (starting with ATG...) # cdna - transcribed sequence (devoid of introns, but containing untranslated exons) # protein - cds tran...2139 days ago
2019 days ago
Install and set up i-adhore for synteny and wgd analysis ! -- step by step --
#Need to download i-adhore-3.0.01.tar.gz from https://wdiceryfd4rjn74bjhhtz2k5di--bioinformatics-psb-ugent-be.translate.goog/webtools/i-adhore/licensing/ #Follow these steps (wg...1203 days ago
94 days ago
Perl and BioPerl script to extract protein sequences using GFF file !
#!/usr/bin/perl use strict; use warnings; use Bio::DB::Fasta; use Bio::SeqIO; # Paths to your GFF file and genome FASTA file my $gff_file = 'path/to/your/file.gf...94 days ago
Perl script for six frame translation !
...= (1, 2, 3, -1, -2, -3); foreach my $frame (@frames) { my $translated_seq = translate_frame($sequence, $frame);...anslation my $revcomp_seq = $sequence->revcom; $translated_seq = $revcomp_seq->translat...94 days ago
Python script for six frame translation of sequences !
...IO from Bio.Seq import Seq def translate_frame(sequence, frame): if frame > 0: translated_seq = sequence[frame-1:].tra...ence.reverse_complement() translated_seq = reverse_complement[abs...): protein_sequence = translate_frame(record.seq, frame)...94 days ago