| 查看: 1885 | 回复: 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。。。   |
回复此楼» 猜你喜欢
蛋白质检测：精准分析，解锁生物分子的密码
已经有0人回复
-
不合理蛙科研实验之小鼠实验：严谨设计，解析生命机制的重要载体
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有233人回复
-
不合理蛙科研实验之重金属检测：精准筛查，守护健康与环境的防线
已经有0人回复
-
不合理蛙科研实验之“蛙测重金属我背锅三千”
已经有0人回复
-
不合理蛙科研实验之“鼠逃三次我发三篇SCI”
已经有0人回复
-
纳米粒度分析
已经有3人回复
-
林科院林化所招收材料与化工专业硕士，坐标南京
已经有0人回复
-
留学--博士招生
已经有15人回复
-
化学工程，本211，求调剂，ab区不限
已经有4人回复
2楼2015-01-15 10:29:15
wizardfan
至尊木虫 (著名写手)
- MolEPI: 4
- 应助: 599 (博士)
- 贵宾: 1.818
- 金币: 24632.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
- 金币: 24632.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
