| 查看: 629 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
用perl程序将fq文件的序列行反向输出已有4人参与
|
|||
|
本人实属perl小小白一枚 现遇到生信方面的一个问题 需要将fastq文件中含序列的行做反向处理，其他行不动，最后输出 这里有师兄帮我写的代码 但是在Linux下运行总是提示错误 大牛能否帮忙一看，小小白不胜感激！ 首次接触木虫，金币为0，攒了金币后会送给帮助我的牛牛们！ 大谢~ 举个栗子 @HWI-ST833:306:C7AYCACXX:1:1101:1491:2080#0/1 CGCTGAGAATGAGCCCCGTGTGGTTGGTGCGCGGACACGCACTGCCTGCGTAACTAGAGGGAGCTGACGGATGACGCCCCCGCGCCACGCCGCT + CCCFFFFFHHHGHJJJJJFHGIJHIJJDGHHIJJIHIJHHGFFEEEEEED?BDDEDDDDDDBDDDDDDDDBBDDDBDDDDDDDDDDDDDDDDDD 这是fq文件的某条数据，需要取出第二行（序列行），反向（TCGCCGCA……AAGAGTCGC），其他行不动，输出酱紫的： @HWI-ST833:306:C7AYCACXX:1:1101:1491:2080#0/1 TCGCCGCA……AAGAGTCGC + CCCFFFFFHHHGHJJJJJFHGIJHIJJDGHHIJJIHIJHHGFFEEEEEED?BDDEDDDDDDBDDDDDDDDBBDDDBDDDDDDDDDDDDDDDDDD 代码如下： #!/usr/bin/perl use warnings; open (FASTQ,'XXX.fastq'); while (<FASTQ> ) { my $tmp1=$_; chomp($tmp2); my $tmp2=<FASTQ>; chomp($tmp2); my @seq=split//,$tmp2; @seq=reverse @seq; my $s=join"",@seq; my $tmp3=<FASTQ>; chomp($tmp3); my $tmp4=<FASTQ>; chomp($tmp4); my @quality=split//,$tmp4; @quality=reverse @quality; my $q=join"",@quality; print "$tmp1\n$s\n$tmp3\n$q\n"; } close FASTQ; Linux提示的错误： Name "main::tmp2" used only once: possible typo at reverse.pl line 6. readline() on closed filehandle FASTQ at reverse.pl line 4. 感谢牛牛的帮助~~ |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有5人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
