24小时热门版块排行榜    

查看: 2535  |  回复: 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的回帖

yuehedou

木虫 (小有名气)

引用回帖:
15楼: Originally posted by XOooZzz at 2014-03-04 09:23:04
呃,要实现这两种方法的话分别有这些问题:
方法1:
困难在于如果序列1和序列2不保证长度相等且碱基位置一一对应的话,首先要写一个两序列比对的算法。这个我做不来。

方法2:
总质量值怎么算?直接相加的话 ...

哦,我都没想到这么多……看来还真不容易,我可能要和公司给我的分析结果一样,只分析一个文件了……
谢谢你非常细心的关注和思考!
其实有时候我又有这样的想法:其实每两个序列的筛选都可以看做是极端的序列拼接操作(最终得到一条序列),如果有这样一个拼接软件,能用它实现两个fastq序列的拼接,而这种拼接时对于匹配区那些错配碱基的取舍是根据其质量值的大小来做决定的,这样一来,我的目的也就达到了。——不过当然我也还没见到过这种软件。
每天都为自己的无知而羞耻!
16楼2014-03-05 09:34:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

中丘子

捐助贵宾 (小有名气)


【答案】应助回帖


感谢参与,应助指数 +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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 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一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +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文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:40 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一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:39 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 两块石头 +6 阿美_Lml888 2026-09-03 8/400 2026-09-05 13:52 by 阿美_Lml888
[硕博家园] 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
信息提示
请填处理意见