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

登录

该问题已被锁定！

1: 关注

773: 浏览

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

生物信息学

我已经筛选出了差异表达基因，但是存在一行有很多个基因的现象，比如某一行有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 人关注

相关问题

在很多分析表达谱芯片的教程中，使用limma 包寻找差异表达基因，都是先找到差异表达探针，再把探针注释为基因id，为何不先进行注释再计算差异表达呢？这两种方法有什么区别呢？: 807 浏览 2 关注 1 回答 0 评论

请问UCSC上下载的参考基因组中这个文件是什么含义: 565 浏览 2 关注 1 回答 0 评论

GEO数据芯片数据基因名转换: 545 浏览 2 关注 1 回答 0 评论

信号通路特征基因集的查询？用于作ssGSEA分析。: 744 浏览 2 关注 1 回答 0 评论

RNA-seq不同样本多个生物学重复不同处理条件下的如何找差异基因: 772 浏览 2 关注 2 回答 0 评论

SRR1924430的sra数据与hdf5提取的fastq不一致: 520 浏览 1 关注 0 回答 0 评论

Bulk-RNAseq多组数据差异表达基因的筛选: 881 浏览 2 关注 2 回答 0 评论

关于tophat2 + cufflinks 得到的差异基因和featureCounts差异基因数目相差太大原因: 858 浏览 2 关注 2 回答 0 评论

请问一下，我的噬菌体基因组fasta文件还是config打头的，是不是需要进一步拼接成scaffold？还是挑选最大的config进行后续分析？: 551 浏览 2 关注 1 回答 0 评论

细菌基因组分析: 742 浏览 2 关注 1 回答 0 评论

推荐内容

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

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

问如何筛选小鼠lncRNA?: 740 浏览 2 关注 1 回答 0 评论

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

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

问生信问题（我也不知道该怎么问）: 1435 浏览 2 关注 1 回答 0 评论

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

问关于cox回归分析问题: 355 浏览 1 关注 0 回答 0 评论

问生信到底如何入门？？？: 1531 浏览 6 关注 6 回答 0 评论

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

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024