[Bioperl-l] divide and blast blastunsplit blast subsequence
rrc22 at cam.ac.uk
Tue Oct 2 09:38:45 EDT 2007
> but isn't there an easier way like i can specify blast subsequence 200-900
> of fasta file and it will return the blastx hits in coords in terms of the
> whole 200kb?
Once you have split up your sequence (as Chris suggested), and run your
BLAST, then you can add the hits to each subsequence as features. The
subsequences can then be re-assembled using the cat method from
Bio::SeqUtils, which will adjust the coordinates of the features
Dr. Roy Chaudhuri
Department of Veterinary Medicine
University of Cambridge, U.K.
More information about the Bioperl-l