5 利用perl脚本提取基因上游的1500bp(指定位置)序列时,输出文件总是空文件。尝试过很多办法无法解决

输入文件1-基因组文件如下图:attachments-2018-08-sCRSQzOF5b8566d214cf2.jpg

输入文件2-位置信息文件如下图:注,没有空格并且txt文本是excl文件另存制表符分隔

attachments-2018-08-iqqwdX8m5b8566e191f93.jpg

输出文件以及所敲指令-如下图:指令未显示报错,但是输出文件是空文件

attachments-2018-08-PVbLz9SZ5b8566ef632d9.jpgattachments-2018-08-Gnk8GIMj5b8566fdcac87.jpg


所用代码:

attachments-2018-08-nkCGUi9P5b8549b94b202.jpgattachments-2018-08-QS1KCkHQ5b8549c3209aa.jpg

所用数据从jgi网站下载

请先 登录 后评论

1 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

通过调试,发现,输入的基因位置文件由于在windows当中用editplus编辑过,因此产生的换行符和linux当中不兼容,导致脚本不识别所以没有结果;

建议不要在windows当中操作修改文件再保存,如果要操作一定要统一编码格式utf-8:具体editplus设置方法如下:

https://www.omicsclass.com/article/395



请先 登录 后评论
  • 1 关注
  • 0 收藏,3430 浏览
  • 忽然慧慧都 提出于 2018-08-28 21:07

相似问题