根据 基因名、bed 文件的基因位置,提取 DNA 序列 bedtools
1、根据 Gene Symbol 查找在序列上的位置
从UCSC下载匹配的 文件
链接: http://genome.ucsc.edu/cgi-bin/hgTables.
1)、主要注意版本的信息,按默认就行
2)、点击 get output,至下一页,就可以选择自己勾选自己想要的信息了
3)、下载后的文件,三列分别是 Transcription位置,Exon位置,geneSymbol
2、根据 基因位置 提取参考上的序列
bedtools getfasta -fi hg19.fa -bed target.bed > target.fa
bed 文件格式:
chr10 131265453 131565783