Perl script to find palindromic regions in DNA sequences
use strict; use warnings; my $pp = qr/(?: (\w) (?1) \g{-1} | \w? )/ix; my $filename = $ARGV[0]; open(my $fh, '2153 days ago
Perl script to find coding regions in DNA sequences
...hould be obtained with the two sequences: high values (positive) for coding protein sequences and low values (zero or even negative) for intronic regions. exon.fa >HUMHBB.2ex...2154 days ago
2074 days ago
Perl script to break the contigs by 'N'
...$seq; #now all bases are in uppercase #Searches for NNNNN regions (scaffold breaks) #Hash nregions stores a hash from the contig...inNlen) { my $end=$i; my $ini=$end-$n+1; push(@regions,"$ini $end"); } $n=0;...2062 days ago
Bash script to extract intronic fragments !
#To obtain introns, we simply need the gene and exonic coordinates; #by subtracting the exonic regions from the genic region, we have the intronic region. gunzip -c genome_file.g...1364 days ago