How I can find sequences

Question

[hosin]

I have some coordinates ( like: chr24:40902129-41709809, chr25:11035146-11097889) from Ovis aries:Oar_v4.0 file. How I can find sequences for them. Thanks

Answer 1

[h.mon]

There are a number of possibilities. For example, download and index the whole genome, and use samtools to extract regions of interest:

samtools faidx genome.fa
samtools faidx genome.fa chr24:40902129-41709809 > region1.fa

Answer 2

[Tm]

If you have bed/gff file of coordinates, then you can use bedtools getfasta and if you have just couple of coordinates to extract the you can try seqkit faidx which is similar to samtools faidx