24小时热门版块排行榜    

查看: 2017  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 正版小辉辉 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿重庆大学085700资源与环境,总分308求调剂 +7 墨墨漠 2026-03-23 8/400 2026-03-23 20:36 by Creta
[考研] 一志愿陕师大生物学071000,298分,求调剂 +3 SYA! 2026-03-23 3/150 2026-03-23 19:09 by macy2011
[考研] 303求调剂 +4 元夕元 2026-03-20 4/200 2026-03-23 19:00 by macy2011
[考研] 求调剂一志愿武汉理工大学材料工程(085601) +3 WW.' 2026-03-23 5/250 2026-03-23 17:18 by 枫翼ljj
[考研] 307求调剂 +11 冷笙123 2026-03-17 11/550 2026-03-22 20:16 by edmund7
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-20 3/150 2026-03-22 16:00 by ColorlessPI
[考研] 化学调剂 +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[考研] 材料工程专硕 348分求调剂 +3 冬辞. 2026-03-17 5/250 2026-03-21 18:47 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
信息提示
请填处理意见