24小时热门版块排行榜    

查看: 8432  |  回复: 21
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

seamas_gao1

银虫 (小有名气)

[求助] 转录组测序结果的“raw data”和“clean data”里为什么全是100bp的序列? 已有3人参与

我测了几个转录组样品,测序公司把“raw data”和“clean data”,都给我了,都是fq格式的文件,每个大概5个G。我打开发现里面全都是100 bp的序列,请问这正常吗?QC和拼接都是用这些100bp的序列吗?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seamas_gao1

银虫 (小有名气)

没有人知道吗?
2楼2014-03-31 11:10:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seamas_gao1

银虫 (小有名气)

引用回帖:
3楼: Originally posted by bio-li at 2014-03-31 12:09:23
拼接用的是clean data,raw data不能用的,你这测序采用的是什么平台

我正是准备用clean data来拼接,用的是Illumina,我只是不明白为什么clean data全都是100 bp的序列,这正常吗?
4楼2014-04-01 08:00:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seamas_gao1

银虫 (小有名气)

引用回帖:
5楼: Originally posted by 洛基殿下~ at 2014-04-01 10:58:42
正常的,如果是转录组、DNA重测序,illumina平台一般都是采用PE100测序,测序的过程中都是测DNA片断2端各100bp,所以下机出来的都是100bp的片断。

raw data是原始下机数据,clean data是经过过滤处理后留下的数 ...

能否麻烦再问一下,只是两端各测100bp的话,那怎么能把序列测通呢,比如序列片段是400bp的话,中间的200bp不就测不到了吗?是否是先把序列都打碎到不到200bp的长度呢,还麻烦解释一下,谢谢!
7楼2014-04-01 15:40:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seamas_gao1

银虫 (小有名气)

引用回帖:
6楼: Originally posted by XOooZzz at 2014-04-01 11:03:50
lz你好,我对生物信息学有点兴趣,想请问你用什么软件进行拼接?

CLC,花钱买的
8楼2014-04-01 15:40:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seamas_gao1

银虫 (小有名气)

引用回帖:
9楼: Originally posted by woody196 at 2014-04-02 08:05:16
如果你使用的是illumina测序的话都是片段为100bp的reads,raw reads就是测序完后的原始序列,而cleanreads是去除接头以及一些质量值较低的reads后产生的。

能否麻烦再问一下,只是两端各测100bp的话,那怎么能把序列测通呢,比如序列片段是400bp的话,中间的200bp不就测不到了吗?是否是先把序列都打碎到不到200bp的长度呢,还麻烦解释一下,谢谢!
10楼2014-04-03 10:06:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seamas_gao1

银虫 (小有名气)

引用回帖:
12楼: Originally posted by 洛基殿下~ at 2014-04-04 09:50:14
一半都是打断成500、300bp。以你的400bp为例,其实他是没有测通的,只是测了2边的各100bp,然后利用片断与片断之间这100bp序列的重叠关系(因为是DNA是随机打断的)利用生物信息软件拼接起来的。...

十分感谢
13楼2014-04-14 08:44:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seamas_gao1

银虫 (小有名气)

引用回帖:
11楼: Originally posted by addict_cy at 2014-04-03 11:25:21
考虑到一次测序的片断数目很多(比如5000万个片断),针对一对reads中间测不到的200bp,总会有其他的reads片断覆盖到,这样从全基因来看,基本上每个位置都会有足够的片断覆盖,从而实现对序列的resequencing....

请问测序时用到的“100倍覆盖”,它的意思是不是就是将完整的转录组随机打断100次,之后再进行拼接呢?
14楼2014-04-14 08:56:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seamas_gao1

银虫 (小有名气)

送红花一朵
引用回帖:
15楼: Originally posted by addict_cy at 2014-04-14 14:22:24
以人的基因组为例,总共3G bp, 如果总测序量是90G(90bp*1G reads),则基因组上每个位点平均被覆盖30次,叫做30X覆盖;注意这不意味着每个位点都有30次覆盖,实际上是一个分布;
转录组的话,由于每个 ...

十分感谢啊
17楼2014-04-15 07:47:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seamas_gao1

银虫 (小有名气)

引用回帖:
16楼: Originally posted by 洛基殿下~ at 2014-04-14 16:14:45
是的,转录组不讲测序深度,因为转录组的表达量不一样,谈不上平均测序深度。主要看饱和曲线,饱和曲线区域饱和,就表明大部分的转录组被测到,再加大测序量,被测到的转录组数量也不会明显增加。  重测序和基因组这 ...

谢谢啊~
18楼2014-04-15 07:48:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 seamas_gao1 的主题更新
信息提示
请填处理意见