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
• 1,340 views
•
link
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
• 1 views
•
link
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
• 1 views
•
link
Log in to answer this question.