24小时热门版块排行榜    

查看: 2193  |  回复: 17

正版小辉辉

铁虫 (初入文坛)

[求助] shell求助---删除特定行中的字符串

我有如下一个文件:
>S6_1
GTAGGGAATCTTCCGCAATGGACGAAAGTCTAACGGAGCAACAGAAACCCGGGTAGTCC
>S6_2
GTAGGGAATCTTCCGCAATGGACGAAAGTCTAGATACCCGTGTAGTCC
>S6_3
GTAGGGAATCTTCCGCAATGGACGAAAGTCTGTTAGATACCCTAGTAGTCC

我想删除每一行(碱基)中的最后17个字符(AGATACCCTAGTAGTCC),但是不知道用语言如何实现?
貌似用sed和awk比较容易。
回复此楼
哈哈
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

前提是每行AGATACCCTAGTAGTCC唯一

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

不要说话
5楼2017-07-13 11:12:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

不要说话
2楼2017-07-12 22:53:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

正版小辉辉

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by ljw4010 at 2017-07-12 22:53:09
sed可以搞

我也明白,可是不会写代码

发自小木虫Android客户端
哈哈
3楼2017-07-12 22:56:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

sed -i 's/AGATACCCTAGTAGTCC//g'   filename

发自小木虫Android客户端
不要说话
4楼2017-07-13 11:11:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

正版小辉辉

铁虫 (初入文坛)

送红花一朵
引用回帖:
5楼: Originally posted by ljw4010 at 2017-07-13 11:12:27
前提是每行AGATACCCTAGTAGTCC唯一

但是如果不一样呢?这样是不是就不好办了?
哈哈
6楼2017-07-13 11:28:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

不一样就这么搞,分组正则,我没验过,你试下,sed -i 's/(.*)AGATACCCTAGTAGTCC$/\1/g'? ?filename

发自小木虫Android客户端
不要说话
7楼2017-07-13 11:39:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

sed -i 's/(.*)AGATACCCTAGTAGTCC$/\1/g' filename

发自小木虫Android客户端
不要说话
8楼2017-07-13 11:40:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

正版小辉辉

铁虫 (初入文坛)

引用回帖:
8楼: Originally posted by ljw4010 at 2017-07-13 11:40:50
sed -i 's/(.*)AGATACCCTAGTAGTCC$/\1/g' filename

qiime@qiime-190-virtual-box:~$ sed -i 's/(.*)ATTAGATACCCTAGTAGTCC$/\1/g' /home/qiime/S1_S14.fasta
sed: -e expression #1, char 36: invalid reference \1 on `s' command's RHS
哈哈
9楼2017-07-13 12:19:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

正版小辉辉

铁虫 (初入文坛)

引用回帖:
9楼: Originally posted by 正版小辉辉 at 2017-07-13 12:19:34
qiime@qiime-190-virtual-box:~$ sed -i 's/(.*)ATTAGATACCCTAGTAGTCC$/\1/g' /home/qiime/S1_S14.fasta
sed: -e expression #1, char 36: invalid reference \1 on `s' command's RHS...

这样看起来,有错误出现
哈哈
10楼2017-07-13 12:20:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 正版小辉辉 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +11 医学老男孩 2026-08-13 29/1450 2026-08-14 00:04 by 大快人心啦
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 28/1400 2026-08-13 17:51 by 医学老男孩
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见