| 查看: 2041 | 回复: 4 | |||
[交流]
血泪求助：Trinity的normalize 已有2人参与
|
|
鄙人最近用Trinity软件normalize一个trim后的数据，奈何一直报错如下： Error, pairs.K25.stats is empty. Be sure to check your fastq reads and ensure that the read names are identical except for the /1 or /2 designation. at /diag/home/aubzyy/trinityrnaseq_r20140717/util/insilico_read_normalization.pl line 865. 这是咋回事啊，俺是新手，跪求大侠出手相助！ 我的数据是Pair-end。。。   |
回复此楼» 猜你喜欢
Wnt-C59：PORCN靶向化学探针，Wnt信号通路阻断机制与科研应用全解析
已经有0人回复
-
Fexaramine：肠道限制性FXR化学探针，结构特征与信号通路研究全解析
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有145人回复
-
NVS-STG2：靶向STING跨膜域的分子胶类激动剂，结构特征与寡聚化机制研究全解析
已经有0人回复
-
DOTA-GPC3肽（CAS 3066439-94-6）：双功能分子探针前体的结构特征与标记应用解析
已经有0人回复
-
科研用 TLR7/8 双激动剂怎么选？现货与质控要点
已经有0人回复
-
SEL120-34A HCl：CDK8/CDK19高选择性激酶抑制剂
已经有0人回复
-
氢氧化钾活化操作
已经有18人回复
-
散金
已经有68人回复
-
生物质超级电容器
已经有0人回复
2楼2015-01-15 10:29:15
wizardfan
至尊木虫 (著名写手)
- MolEPI: 4
- 应助: 599 (博士)
- 贵宾: 1.818
- 金币: 24648.2
- 散金: 197
- 红花: 48
- 沙发: 2
- 帖子: 2254
- 在线: 400.7小时
- 虫号: 1879241
- 注册: 2012-07-05
- 性别: GG
- 专业: 生物信息学
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
it says that the name of sequences in your fastq can not match In fastq format, Line 1 begins with a '@' character and is followed by a sequence identifier and an optional description (like a FASTA title line). e.g. @HWI-EAS209_0006_FC706VJ:5:58:5894:21141#ATCACG/1 another entry would be something like @HWI-EAS209_0006_FC706VJ:5:58:5894:21141#ATCACG/2 |
3楼2015-01-16 08:09:27
|
Thank you. But my sequence starts like this: @SRR1559285.21 HWI-ST531R:238:C346VACXX:5:1101:1193:2038 length=101 CCTTCATTCATGTGCTCACGCCAGTATGCAGCATTTTTAACAGCACAGTCCCAGT +SRR1559285.21 HWI-ST531R:238:C346VACXX:5:1101:1193:2038 length=101 @CCFFFFFGHHGHIIJJJ4FHIJJGIIJJIIJIJJJJJHIJJJIIJJJHHJJJJI @SRR1559285.22 HWI-ST531R:238:C346VACXX:5:1101:1126:2077 length=101 CCTATTCATTTACAGCAGAAGTGCTATTTCTTATTCATGATGTACAAGTCGTTGTTCCGTGAAAATGAAACAGAGACGGG CGTCGTAGGATCCTTAGACGG +SRR1559285.22 HWI-ST531R:238:C346VACXX:5:1101:1126:2077 length=101 CCCFFFFFHHGHHFIJJJHIIIHIJJJIJIJJIHJGGIFHIIGGIIJJFGHGEGHGIGGDGGGJJIIJGHIIIHJJHHHF FD?BDBBBDCCCDDDDC:@@B @SRR1559285.23 HWI-ST531R:238:C346VACXX:5:1101:1228:2083 length=101 CAGACGAATTATAGAAAGGCTGTTTAATTAGAAATTTACTTAATCCAAGTTTTGGCTGAAATGGATTGAAGTTTGTCCTG GAATGTTTCGAACACATCTGT +SRR1559285.23 HWI-ST531R:238:C346VACXX:5:1101:1228:2083 length=101 CCCFFFFFHFHGFIBHIJIJJJJGIGAHIEIHIHHIHIFIIIIIJJJEHHGIJIIHIIGJGIII@FHJGIGCAHEEHHFF FFFCEEEECDDDBDDDDCD:@ @SRR1559285.24 HWI-ST531R:238:C346VACXX:5:1101:1246:2099 length=101 TGAAACTGAGGAATGTATGCCTCTCACGGTGACCCTGTCTCACCGTCTGAACCAAAAAATTGCTGAGCTTCGTAGAGATG GTACATTCTGGTGGGCAAGAC +SRR1559285.24 HWI-ST531R:238:C346VACXX:5:1101:1246:2099 length=101 @@@DF;DDBFFFFGGGHIJJAHGHIGEHGACGIJFIGIIGGIJJJHGIJIIJJEIIIIGGCACEEEFDEFFFDEC;@CCC D@DCDFDDEEDCDDDB?@B<9 @SRR1559285.25 HWI-ST531R:238:C346VACXX:5:1101:1179:2100 length=101 Is this okay? I am a first year graduate, thank you for your patience. |
4楼2015-01-19 13:40:24
wizardfan
至尊木虫 (著名写手)
- MolEPI: 4
- 应助: 599 (博士)
- 贵宾: 1.818
- 金币: 24648.2
- 散金: 197
- 红花: 48
- 沙发: 2
- 帖子: 2254
- 在线: 400.7小时
- 虫号: 1879241
- 注册: 2012-07-05
- 性别: GG
- 专业: 生物信息学
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
how sure you are that this is a pair-end data? Normally two files for pair-end data, named like ***_1.fastaq and ***_2.fastq. In ***_1.fastq, the headers of sequences are like ********/1, so as in the second file. There is a public server for you to use Trinity https://galaxy.cbio.mskcc.org/ |
5楼2015-01-20 10:18:19
5