美文网首页一起生信啦啦啦
gencode和UCSC的fasta相同

gencode和UCSC的fasta相同

作者: 苏牧传媒 | 来源:发表于2018-11-01 22:24 被阅读7次

1.构建一个tq.bed:

chr1 33344 33388

chr4    77777  77799

2.用bedtools分别对gencode和UCSC的fasta进行提取:

bedtools getfasta -fi GRCh37.chr.fa -bed tq.bed -split -name | fold -w 60 > tq.fa

gencode

bedtools getfasta -fi ../UCSC_hg19/hg19.fa -bed tq.bed -split -name | fold -w 60 > tq2.fa

UCSC

结果显示序列是相同的!

用samtools faidx xxx.fasta序列大小也是一样的!

相关文章

网友评论

    本文标题:gencode和UCSC的fasta相同

    本文链接:https://www.haomeiwen.com/subject/orpfxqtx.html