in reply to extract sequence given positions from fasta
Yes, you will want to read your sequence ID/positional data into an array, ideally an array of hashes were the hash keys would be something like "id", "start, "end".
Then you should probably use Bio::DB::Fasta to retrieve fragments of sequences from your FASTA file, but you could also use Bio::SeqIO and extract regions from your sequences using "substr" on the DNA string.
|
|---|