This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How I can find sequences

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

sequence

2 answers

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

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

Log in to answer this question.