extracting seq from viral genome

[honey]

I have coordinates of a viral genome and would like to extract seq from viral genome. Is there a tool or easy way to upload the txt file of ordinates and it will give me seq (equal to what do in biomart or UCSC genome page for human genes)

Thanks

honey

[dariober]

I have coordinates of a viral genome and would like to extract seq from viral genome. Is there a tool or easy way to upload the txt file of ordinates and it will give me seq (equal to what do in biomart or UCSC genome page for human genes)

Thanks

honey

If you have the reference genome in FASTA format and the coordinates in BED format (chrom <tab> start <tab> end) you can you use bedtools getfasta:

bedtools getfasta

Tool:    bedtools getfasta (aka fastaFromBed)
Version: v2.17.0
Summary: Extract DNA sequences into a fasta file based on feature coordinates.

Usage:   bedtools getfasta [OPTIONS] -fi <fasta> -bed <bed/gff/vcf> -fo <fasta> 

Options: 
	-fi	Input FASTA file
	-bed	BED/GFF/VCF file of ranges to extract from -fi
	-fo	Output file (can be FASTA or TAB-delimited)
	-name	Use the name field for the FASTA header
	-split	given BED12 fmt., extract and concatenate the sequencesfrom the BED "blocks" (e.g., exons)
	-tab	Write output in TAB delimited format.
		- Default is FASTA format.

	-s	Force strandedness. If the feature occupies the antisense,
		strand, the sequence will be reverse complemented.
		- By default, strand information is ignored.

Does it help?
Dario

[swbarnes2]

samtools faidx will do that too, but one at a time. So embed it in some kind of shell script to get a bunch.

[honey]

extracting viral seq

Thanks that was helpful !!!