Question: How I can find sequences
0
gravatar for hosein_salehi6
7 months ago by
hosein_salehi60 wrote:

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 • 229 views
ADD COMMENTlink modified 7 months ago by toralmanvar750 • written 7 months ago by hosein_salehi60
1
gravatar for h.mon
7 months ago by
h.mon23k
Brazil
h.mon23k wrote:

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
ADD COMMENTlink written 7 months ago by h.mon23k
1
gravatar for toralmanvar
7 months ago by
toralmanvar750
toralmanvar750 wrote:

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

ADD COMMENTlink modified 7 months ago • written 7 months ago by toralmanvar750
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2294 users visited in the last hour