[Bioperl-l] divide and blast blastunsplit blast subsequence

Roy Chaudhuri rrc22 at cam.ac.uk
Tue Oct 2 13:38:45 UTC 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 
appropriately.

Roy.
--
Dr. Roy Chaudhuri
Department of Veterinary Medicine
University of Cambridge, U.K.



More information about the Bioperl-l mailing list