- 论坛徽章:
- 8
|
本帖最后由 xiumu2280 于 2013-09-25 10:45 编辑
其实楼主就是想要个去除测序接头的程序···- #!/usr/bin/perl
- use strict;
- use warnings;
- while (<DATA>) {
- chomp;
- my $line = $_;
- my $line1 = <DATA>;
- my $line2 = <DATA>;
- my $line3 = <DATA>;
- if ($line1=~/^ATCA/) {
- $line1=substr ($line1,4);
- $line3=substr ($line3,4);
- }
- print "$line\n$line1$line2$line3";
- }
- __DATA__
- @HWI-ST531R:212:C22J6ACXX:8:2314:4159:37409
- AAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGAT
- +
- C@CFFFFFHHHHHIIHGIIJGIJIBGEIIGJCGHIJIJGIJJJ8>EDDFCED6@ACCCDDD@@BDDDDDDCDDCDCCCCADDDDECDDDCDDDCCDCDCAC
- @HWI-ST531R:212:C22J6ACXX:8:2102:3421:60710
- AAAAAAAAACCATCCAAATCTGGATGGCTTTTCATAATTCTGAGAAATTAGCTGCGCTGGCGCACCGCTTCAAATAAGCAAATTCCGGTCGCAACCGAAAC
- +
- @@CFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJIJIIJJIJJIJIJJJJIIIJJHHFFDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD
- @HWI-ST531R:212:C22J6ACXX:8:1303:19050:54218
- AAAAAAAAACCGCACCAGAATTCCCATCGATCCAAGGCCGTTGGTGTTCCTTTCCCCACAAACCGCCCCAACGTCTCAATCGAGACTCCCCCTCACAGCAG
- +
- ??<BDDDDFFFF08@FFBG>D8?FFFF9;;==5@@=(66:55?B(5;?BA;AA>AA<<;<?B9?5@08-<@BB<2<9(4+:28+25@:>8>B@@B<A>?##
- @HWI-ST531R:212:C22J6ACXX:8:2108:4215:95315 1:N:0:GAGTGG
- AAAAAAAAACCGGGCAATTGCCCGGTGCTGGAAGGGATCTCGTTATTCGGGGGAATAAGGTAGATGTGGATAGTCCAGCCAGTGCGTCAAATAGTGAGAGA
- +
- <<@DDDDDAH@DHG:FHFGIIGCGI@FHHIHIIFHE6?1;?@?==ACDCA@BB5<8>C@C33:@@A@>>:>:?::>4>9<<A::>9(55>A:>A4(+4@32
- @HWI-ST531R:212:C22J6ACXX:8:2208:5557:5216 1:N:0:GAGTGG
- AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAACTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTT
- +
- @@CFFFFFFDFFFEHGIGIHGHHGIIGEDEHFEHFBC@?DBA(;@CCDDC@CDCCCDDCDDD;DB7C>C@BB05?A>@CDDDDCCDDDDCCCCDDDDECDD
- @HWI-ST531R:212:C22J6ACXX:8:2214:1365:3705
- AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTT
- +
- @@@FDDFFF:DHFHGIGHEHEFGEAFHGIIHGCEDFFDE@A@>;ACDDD:>@C(>CD?@@CD@B@BDCA>;;0;B>CC@C<CCDCDAA>CACCCDDEC>CB
- @HWI-ST531R:212:C22J6ACXX:8:1107:2469:40337
- ATCAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTT
- +
- CCCFFFFFHHHHHJJJJJJIJIJJJJJJJJHGHGFFEFEEECEEEEDDDCDDDACDDDDDDDDDDDDDDDDDBDDDDDDDDDDDDDDCDDDEDEEDEEEED
- @HWI-ST531R:212:C22J6ACXX:8:2211:2015:43926
- AAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTT
- +
- CCCFFFFFFHHAFGGHDHIIIIIGHGGHGGHHAHEDFFFCAECCCEDDD>:CCCCDDDDDDD7@BDDCCDBB;B@CCDDCCDDDDDDCCCADDDEDEDEEE
复制代码 |
|