sequence retrieval in R using fasta file
I have a Hg19.fasta file and I have a csv file which contains the genomic coordinates of all 23 chromosomes which the sequence is to be retrieved. How can I do it using R?
• 3,458 views
•
link
updated
by
Jean-Karim Heriche
•
written
by
rsabrina93 •
•
written
by
rsabrina93 •
1 answer
You could use the R package seqinr to read the FASTA file with read.fasta(), read the csv file with read.csv() then subset the sequence using the coordinates. There may also be something for you in the Bioconductor GenomicRanges package.
• 0 views
•
link
written
by
Jean-Karim Heriche
Log in to answer this question.
More posts like this
-
where to obtain hg19 reference transcriptome in gencode
written by Josh •
-
How to extract FASTA headers in R
written by WUSCHEL
-
Extracting genomic coordinates for a Variant calling pipeline
written by salman_96 •
-
Create an assembly with a ribosomic chromosome
written by ra2967 •
-
Visualize genome and blast hits
written by krzysieksobieszek •
-
How to splice introne sequences for each transcript.
written by Tim Padvitski •
-
unplaced, unlocalized, ... scaffolds from hg19 to GRCh38
written by tabone.olivier •
-
Identify restriction sites of a cDNA sequence using a dictionary
written by albamaria.muniesa •
-
Retrieving Genomic Sequence Based On Coordinate Information From Annotation File
written by es55spam •
-
Cuffmerge Is Looking For Contig Fasta Files That I Do Not Have!
written by jobinv
Not what you ask for, but
bedtools getfastawould be the way to do this if you wouldn't constrain to R.