This is a test version of Biostars. For the public version, visit https://www.biostars.org.
extract out intronic sequences from genome fasta file

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
genome

1 answer

Log in to answer this question.