"Michael S. Robeson II" <[EMAIL PROTECTED]> wrote in message ...
> open(DNA_SEQ, $dna_seq)
Due to precedence, the parens are optional here.
> open(OUTFILE, ">indel_list_"."$dna_seq")
^
">indel_list_$dna_seq" or ">indel_list_".$dna_seq
> foreach (keys
I cleaned up the code a little. So, here it is for anyone interested:
#!usr/bin/perl
# By Michael S. Robeson II with the help from the folks at lernperl.org
and bioperl.org
# 10/16/2004
# Last updated: 10/17/2004
# This script was made for the purpose of searching for indels (gaps)
in aligned
# D
I just wanted to thank everyone for their help and suggestions. This is
the final full working code to count continuos gaps in a every sequence
in a multi-sequence FASTA file. It may not be elegant but it is fast
and works well. Sorry for the long post but I wanted to share this with
those that
