Tweet
If I have sequence start and end coordinates, how can I get the sequence from local hg19.fa file. Either using samtools or perl will be good for me.
hg19 has been indexed using samtools.
Thanks!
hg19 has been indexed using samtools.
Thanks!
-
-
You can use the same command you used to index the fasta. That is, construct the index:Originally posted by lewewoo View Post
Then extract your subsequence:Code:samtools faidx hg19.fa
You'll have to adjust the coordinates and chromosome for what you want. There is also a Perl package for working with SAM/BAM files that may be helpful if you have a lot of operations to perform.Code:samtools faidx hg19.fa 1:20-200
-
thanks
will try it
Comment
-
Or you can use seqret in Emboss. You do not need to index your reference.Comment
-
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...04-22-2024, 07:01 AM -
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...04-04-2024, 04:25 PM
Comment