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
-
CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).
Despite this, “CRISPR helped turn genome editing from a specialized technique into...07-31-2026, 11:01 AM
Comment