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

登录

该问题已被锁定！

1: 关注

1301: 浏览

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

生物信息学

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

相关问题

基因组组装问题: 1802 浏览 2 关注 1 回答 0 评论

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

参考基因组添加外源基因序列进行比对: 2356 浏览 2 关注 5 回答 0 评论

基因组文件中chrUn具体是什么意思呢: 1068 浏览 0 关注 0 回答 1 评论

singularity查到不到指定输入文件位置: 1916 浏览 2 关注 1 回答 0 评论

log2(fold_change)可以用来做热图吗？遇到想要表达的目的基因log2(fold_change)值为0时，怎么办？: 1751 浏览 2 关注 1 回答 0 评论

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

从多个.fa文件中提取以“poptri.”开头的蛋白序列: 1108 浏览 2 关注 1 回答 0 评论

grep 能不能指定匹配第多少列？: 1623 浏览 3 关注 2 回答 1 评论

如何按一个列表对基因型文件进行过滤，剔除不需要的样本？: 1715 浏览 2 关注 1 回答 0 评论

推荐内容

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

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

问如何根据现有序列找到它们在三代基因组的编号: 1705 浏览 2 关注 1 回答 0 评论

问请问在不考虑测序质量的情况下，如何根据bam文件获得比对到指定position的碱基？: 2676 浏览 2 关注 2 回答 0 评论

问去除批次效应选定靶标细胞群后如何进行亚群分析？: 2093 浏览 2 关注 2 回答 0 评论

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

问考研生物信息学院校推荐: 2166 浏览 4 关注 3 回答 0 评论

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

问怎么准确提取某一基因的ATG起始密码子位点: 1780 浏览 2 关注 1 回答 0 评论

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

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025