| 查看: 675 | 回复: 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. 感谢牛牛的帮助~~ |
回复此楼» 猜你喜欢
303求调剂
已经有6人回复
-
07化学303求调剂
已经有5人回复
-
332求调剂
已经有6人回复
-
考研调剂
已经有7人回复
-
0856求调剂
已经有5人回复
-
总分293求调剂
已经有6人回复
-
275求调剂
已经有5人回复
-
085600 材料与化工 329分求调剂
已经有7人回复
-
材料与化工328分调剂
已经有6人回复
-
291 求调剂
已经有6人回复
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
2楼2015-10-03 20:40:12
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
3楼2015-10-03 20:46:08
4楼2015-10-13 20:04:49
glssg
木虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1973.5
- 散金: 343
- 红花: 2
- 帖子: 846
- 在线: 189.8小时
- 虫号: 1962301
- 注册: 2012-08-29
- 性别: GG
- 专业: 肿瘤学
5楼2015-10-14 08:04:00
HappyKuan
银虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 195.7
- 散金: 178
- 帖子: 86
- 在线: 75.9小时
- 虫号: 2219965
- 注册: 2013-01-03
- 专业: 生物信息学
6楼2015-10-14 14:02:50
