免费注册 查看新帖 |

Chinaunix

  平台 论坛 博客 文库
最近访问板块 发新帖
查看: 1979 | 回复: 6
打印 上一主题 下一主题

[文本处理] 测序数据分析 [复制链接]

论坛徽章:
0
跳转到指定楼层
1 [收藏(0)] [报告]
发表于 2014-03-18 11:35 |只看该作者 |倒序浏览
各位大神,我想请教一个问题,
我有一个文本,格式如下,我想去掉第一个红色字符串左边的所有字母,第二个红色字符串右边的所有字母,保留中间的序列,红色的字符串是唯一的,没有重复,谢谢了
>K01_313_324
TGATTCGAGCTCGGTACCCGGGGAACTGAGGACACAGCTGTGTATTATCC
GCCAGAGAGAGATACGGGGAAGCTGCTTTTGACTACTGGGGTAAAGGCACACAAGTGACCGTT
CTTTTTTTCAACTCTTTTAGTAGCCTCTATTACCTTGTGAGGTGCAGTAGTA
>K02_313_324
.......
>K03_313_324

论坛徽章:
23
15-16赛季CBA联赛之吉林
日期:2017-12-21 16:39:27白羊座
日期:2014-10-27 11:14:37申猴
日期:2014-10-23 08:36:23金牛座
日期:2014-09-30 08:26:49午马
日期:2014-09-29 09:40:16射手座
日期:2014-11-25 08:56:112015年辞旧岁徽章
日期:2015-03-03 16:54:152015年迎新春徽章
日期:2015-03-04 09:49:0315-16赛季CBA联赛之山东
日期:2017-12-21 16:39:1915-16赛季CBA联赛之广东
日期:2016-01-19 13:33:372015亚冠之山东鲁能
日期:2015-10-13 09:39:062015亚冠之西悉尼流浪者
日期:2015-09-21 08:27:57
2 [报告]
发表于 2014-03-18 11:41 |只看该作者
回复 1# pdzcowboy


这个意思?
  1. sed 's/.*\(ACTGAGGACACAGCTGTGTATTA\)/\1/;s/\(ATTACCTTGTGAGGTGCAGTAGT\).*/\1/' file
复制代码

论坛徽章:
145
技术图书徽章
日期:2013-10-01 15:32:13戌狗
日期:2013-10-25 13:31:35金牛座
日期:2013-11-04 16:22:07子鼠
日期:2013-11-18 18:48:57白羊座
日期:2013-11-29 10:09:11狮子座
日期:2013-12-12 09:57:42白羊座
日期:2013-12-24 16:24:46辰龙
日期:2014-01-08 15:26:12技术图书徽章
日期:2014-01-17 13:24:40巳蛇
日期:2014-02-18 14:32:59未羊
日期:2014-02-20 14:12:13白羊座
日期:2014-02-26 12:06:59
3 [报告]
发表于 2014-03-18 12:09 |只看该作者
回复 1# pdzcowboy

# awk '{if(t=index($1,"ACTGAGGACACAGCTGTGTATTA"))$1=substr($1,t);if(t=index($1,"ATTACCTTGTGAGGTGCAGTAGT"))$1=substr($1,1,t+22)}1' FILE
>K01_313_324
ACTGAGGACACAGCTGTGTATTATCC
GCCAGAGAGAGATACGGGGAAGCTGCTTTTGACTACTGGGGTAAAGGCACACAAGTGACCGTT
CTTTTTTTCAACTCTTTTAGTAGCCTCTATTACCTTGTGAGGTGCAGTAGT
>K02_313_324
.......
>K03_313_324

   

论坛徽章:
0
4 [报告]
发表于 2014-03-18 13:17 |只看该作者
这红色字符串的上面可能还有好几行,这都得删掉回复 2# ly5066113


   

论坛徽章:
0
5 [报告]
发表于 2014-03-18 13:18 |只看该作者
这红色字符串的上面可能还有好几行,这都得删掉回复 3# jason680


   

论坛徽章:
32
处女座
日期:2013-11-20 23:41:20双子座
日期:2014-06-11 17:20:43戌狗
日期:2014-06-16 11:05:00处女座
日期:2014-07-22 17:30:47狮子座
日期:2014-07-28 15:38:17金牛座
日期:2014-08-05 16:34:01亥猪
日期:2014-08-18 13:34:25白羊座
日期:2014-09-02 15:03:55金牛座
日期:2014-11-10 10:23:58处女座
日期:2014-12-02 09:17:52程序设计版块每日发帖之星
日期:2015-06-16 22:20:002015亚冠之塔什干火车头
日期:2015-06-20 23:28:22
6 [报告]
发表于 2014-03-18 13:26 |只看该作者
这样?
  1. sed -r ':1;N;$!b1;s/.*(ACTGAGGACACAGCTGTGTATTA.*ATTACCTTGTGAGGTGCAGTAGT).*/\1/'
复制代码

论坛徽章:
145
技术图书徽章
日期:2013-10-01 15:32:13戌狗
日期:2013-10-25 13:31:35金牛座
日期:2013-11-04 16:22:07子鼠
日期:2013-11-18 18:48:57白羊座
日期:2013-11-29 10:09:11狮子座
日期:2013-12-12 09:57:42白羊座
日期:2013-12-24 16:24:46辰龙
日期:2014-01-08 15:26:12技术图书徽章
日期:2014-01-17 13:24:40巳蛇
日期:2014-02-18 14:32:59未羊
日期:2014-02-20 14:12:13白羊座
日期:2014-02-26 12:06:59
7 [报告]
发表于 2014-03-18 13:49 |只看该作者
回复 5# pdzcowboy

provide your example answer
   
您需要登录后才可以回帖 登录 | 注册

本版积分规则 发表回复

  

北京盛拓优讯信息技术有限公司. 版权所有 京ICP备16024965号-6 北京市公安局海淀分局网监中心备案编号:11010802020122 niuxiaotong@pcpop.com 17352615567
未成年举报专区
中国互联网协会会员  联系我们:huangweiwei@itpub.net
感谢所有关心和支持过ChinaUnix的朋友们 转载本站内容请注明原作者名及出处

清除 Cookies - ChinaUnix - Archiver - WAP - TOP