24小时热门版块排行榜    

查看: 2027  |  回复: 18

yuehedou

木虫 (小有名气)

[求助] 有人做过或接触过microRNA的paired-end测序数据分析吗

我的两个植物样,送去做microRNA的illumina HiSeq测序,公司给我送回来的数据是paired-end类型的,每个样品都有两个同样大小的序列文件,分别存储microRNA的序列及其反向互补序列。
现在在前期过滤序列时遇到了问题了,同一个样品的两份数据因为其测序质量不一样,过滤掉的reads就可能不一样。我的问题是后续分析都只需一个样品的一份数据,这样在最终获得clean reads时,就需要有所取舍,现在不知道该怎么办了,求教各位大牛赐教?!
我现在有三种想法:1. 保留一份数据;2. 把两份数据加和起来;3. 保留两个文件中相同的序列的一份,外加不同的质量高的那条序列(这个难度比较大,还没想好怎么做)。
以前看过microRNA的文献资料描述的基本都是single end的,很少看到paired-end的,如果有人有这方面的文献资料或好的软件,也可以分享给我啊,先谢谢诸位了!
回复此楼
每天都为自己的无知而羞耻!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

XOooZzz

银虫 (小有名气)

能不能贴一小段上来看看是什么样子的?我最近在学编程,看能不能给你写个小程序
11楼2014-03-02 12:41:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

中丘子

捐助贵宾 (小有名气)


【答案】应助回帖


感谢参与,应助指数 +1
yuehedou: 金币+1, 有帮助, 感谢参与 2014-02-28 21:13:36
你们为什么不直接让测序公司帮你们分析一下呢,至少可以当面向他们请教下吧。他们应是有这个义务的。

[ 发自手机版 http://muchong.com/3g ]
2楼2014-02-28 16:53:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

中丘子

捐助贵宾 (小有名气)


【答案】应助回帖

你们为什么不直接让测序公司帮你们分析一下呢,至少可以当面向他们请教下吧。他们应是有这个义务的。

[ 发自手机版 http://muchong.com/3g ]
3楼2014-02-28 16:54:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaoyang636

木虫 (著名写手)

请问为啥microRNA 还要测PE呢?
4楼2014-02-28 21:02:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuehedou

木虫 (小有名气)

引用回帖:
2楼: Originally posted by 中丘子 at 2014-02-28 16:53:33
你们为什么不直接让测序公司帮你们分析一下呢,至少可以当面向他们请教下吧。他们应是有这个义务的。

这个他们当时大概分析了一下,但感觉他们给的数据和图片都很粗糙,于是想自己分析一下。
向公司请教这个倒是可以有,我试试。
每天都为自己的无知而羞耻!
5楼2014-02-28 21:12:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuehedou

木虫 (小有名气)

引用回帖:
4楼: Originally posted by gaoyang636 at 2014-02-28 21:02:55
请问为啥microRNA 还要测PE呢?

失误吧,老板当时不太懂,直接交给公司的,没想到拿到数据是paired-end的
每天都为自己的无知而羞耻!
6楼2014-02-28 21:15:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

XOooZzz

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
我想向lz请教一下,paired end的microRNA数据是怎么样的呢?测的是microRNA成熟体?那读长不是很短?
7楼2014-03-01 00:08:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuehedou

木虫 (小有名气)

引用回帖:
7楼: Originally posted by XOooZzz at 2014-03-01 00:08:24
我想向lz请教一下,paired end的microRNA数据是怎么样的呢?测的是microRNA成熟体?那读长不是很短?

是的,读长很短,两个文件数据分别去掉接头后,对应序列等长且反向互补。
每天都为自己的无知而羞耻!
8楼2014-03-01 11:46:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

XOooZzz

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★
yuehedou: 金币+1, 有帮助, 多谢关注! 2014-03-02 11:13:01
yuehedou: 金币+2, 有帮助, 加油 2014-03-03 15:05:53
西门吹雪170: 金币+1, 鼓励交流 2014-03-03 18:38:24
yuehedou: 金币+5, ★★★很有帮助, 非常感谢参与讨论! 2014-03-05 09:39:24
我没接触过,不懂,但同意方法3。
“保留两个文件中相同的序列的一份”这一步应该不难;
“外加不同的质量高的那条序列”这一步,不知道是怎么看到测序质量的信息的?
9楼2014-03-01 12:47:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuehedou

木虫 (小有名气)

引用回帖:
9楼: Originally posted by XOooZzz at 2014-03-01 12:47:33
我没接触过,不懂,但同意方法3。
“保留两个文件中相同的序列的一份”这一步应该不难;
“外加不同的质量高的那条序列”这一步,不知道是怎么看到测序质量的信息的?

fastq文件,第四行就是质量值啊!
想这样做,但没找到具有这种功能的软件……自己的编程能力恐怕还做不到这一点,纠结中……
每天都为自己的无知而羞耻!
10楼2014-03-02 11:12:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yuehedou 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 情人节自我反思:在爱情中有过遗憾吗? +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错,点赞! +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题,请指教! 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
请填处理意见