想要提取指定基因的行信息 - 做专业的生信问答网站

登录

该问题已被锁定！

1: 关注

1834: 浏览

想要提取指定基因的行信息

生物信息学

我已经筛选出了差异表达基因，但是存在一行有很多个基因的现象，比如某一行有RPL13AP5///RPL13AP6///SNORD32A///SNORD33///SNORD34///SNORD35A///RPL13A 所以我想用shell脚本把这样的每个基因弄成一行，试了好多种办法，还是没有达到想要的效果，请问各位前辈有没有什么办法可以提供一下？

好问题 0 评论收藏举报

1 回答

minipeach 前台管理员用户来自于: 重庆市
2018-09-03 14:05

name=`cat /home/Taoct_16/anaysis/ICM/deal/name_deal.txt` for dd in $name do grep $dd 5406.txt >> result1.txt done 这是我提取信息的脚本，其中name_deal.txt是我要提取的所有的基因的名字，5406.txt是我要从中提取的文件

赞同 0 1评论

关于作者

: minipeach 前台管理员
这家伙很懒，还没有设置简介

10: 回答

0: 文章

16: 问题

问题动态

发布时间: 2018-09-03 14:03

更新时间: 2018-09-03 14:06

关注人数: 1 人关注

相关问题

bowtie2 参考基因组注释比对: 3037 浏览 2 关注 2 回答 0 评论

创建NT子库以及NT库提取特定物种分类的序列: 2553 浏览 3 关注 2 回答 0 评论

enhancer结合ATAC-seq进行下游靶基因的寻找策略: 2087 浏览 2 关注 1 回答 0 评论

进行转录组数据分析时，进行cuffdiff后的输出文件gene_exp.diff中，一个基因出现了两个不同的表达量数据，应该如何处理？: 1700 浏览 2 关注 1 回答 0 评论

转录本坐标转换成基因组坐标: 3083 浏览 4 关注 2 回答 0 评论

为什么差异基因分析会出现被KO的基因log2FC的数值大于1的情况: 2679 浏览 2 关注 1 回答 0 评论

从bam文件还能提取出来原来比对前的fastaq文件吗？: 2665 浏览 2 关注 2 回答 0 评论

为啥NCBI和ensemble相同版本的基因组注释有这么大的差别？: 2469 浏览 2 关注 2 回答 0 评论

一般报告基因组某个碱基突变是指正链还是负链: 1609 浏览 2 关注 1 回答 0 评论

关于基因对某种疾病的影响的问题: 1459 浏览 2 关注 1 回答 0 评论

推荐内容

问 sambamba使用: 2120 浏览 2 关注 1 回答 0 评论

问生存分析KM-plot交叉问题: 2630 浏览 2 关注 2 回答 0 评论

问翻译序列: 1489 浏览 2 关注 1 回答 0 评论

问 infercnv运行报错: 5126 浏览 2 关注 1 回答 0 评论

问关于蛋白质文件的疑问: 2152 浏览 2 关注 1 回答 0 评论

问关于基因间的相关性分析: 2442 浏览 2 关注 2 回答 0 评论

问 riboseq中frame指的是什么，是如何确定的，通过riboseq如何知道一个基因的开放阅读框: 2060 浏览 2 关注 1 回答 0 评论

问关于GO分析的问题: 1644 浏览 2 关注 1 回答 0 评论

问关于基因对某种疾病的影响的问题: 1459 浏览 2 关注 1 回答 0 评论

问根据GFF和fasta等文件提取某一基因的ATG位点信息，比如具体的位置？: 1788 浏览 2 关注 1 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025