sh 取终端基因C. elegans
Posted
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了sh 取终端基因C. elegans相关的知识，希望对你有一定的参考价值。
end_dist=30000
chrom_I=15072434
chrom_II=15279421
chrom_III=13783801
chrom_IV=17493829
chrom_V=20924180
chrom_X=17718942
tabix gene.gff.gz I:1-$end_dist > terminal_genes.txt
tabix gene.gff.gz II:1-$end_dist >> terminal_genes.txt
tabix gene.gff.gz III:1-$end_dist >> terminal_genes.txt
tabix gene.gff.gz IV:1-$end_dist >> terminal_genes.txt
tabix gene.gff.gz V:1-$end_dist >> terminal_genes.txt
tabix gene.gff.gz X:1-$end_dist >> terminal_genes.txt
tabix gene.gff.gz I:`expr $chrom_I - $end_dist`-$chrom_I >> terminal_genes.txt
tabix gene.gff.gz I:`expr $chrom_II - $end_dist`-$chrom_II >> terminal_genes.txt
tabix gene.gff.gz I:`expr $chrom_III - $end_dist`-$chrom_III >> terminal_genes.txt
tabix gene.gff.gz I:`expr $chrom_IV - $end_dist`-$chrom_IV >> terminal_genes.txt
tabix gene.gff.gz I:`expr $chrom_V - $end_dist`-$chrom_V >> terminal_genes.txt
tabix gene.gff.gz I:`expr $chrom_X - $end_dist`-$chrom_X >> terminal_genes.txt
cat terminal_genes.txt | grep 'Gene' | cut -f 9 | egrep -o 'Alias=([^,]+)'
以上是关于sh 取终端基因C. elegans的主要内容，如果未能解决你的问题，请参考以下文章
sh 得到含糊不清的等位基因
sh 重叠群和基因重命名
sh 获得1000G亚群的等位基因频率
sh 使用登录号从ncbi下载基因组
sh 从gt_hmm生成断点基因型
sh 将基因组文件拆分为单独的读数。 #fasta