24小时热门版块排行榜    

查看: 2544  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 19:21 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 6/300 2026-09-06 18:41 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 18:09 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:01 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才!!! +5 SCSIOyxguo 2026-09-03 8/400 2026-09-06 15:06 by 可爱小花朵
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见