Entering edit mode
7.9 years ago
a.rex
▴
350
I have a bed file that contains the start and end co-ordinates of intronic sequences. How can I use this bed file to extract out the dna sequences of these same introns from a genome fasta file?
The bedfile of intron looks like this:
scaffold1 206 262 transcriptA . + intron
scaffold1 749 1391 transcriptA . + intron
scaffold1 1549 1595 transcriptA . + intron
scaffold1 2090 2657 transcriptA . + intron
scaffold1 3013 3065 transcriptA . + intron
scaffold1 3315 3371 transcriptA . + intron
scaffold1 3930 4522 transcriptA . + intron
scaffold1 4648 4699 transcriptA . + intron