Mercurial > repos > mingchen0919 > statonlab_fasta_splitter
comparison split.pl @ 2:38435c74e795 draft default tip
working version
| author | mingchen0919 |
|---|---|
| date | Mon, 09 Apr 2018 14:20:42 -0400 |
| parents | 3918528c19b4 |
| children |
comparison
equal
deleted
inserted
replaced
| 1:ecf34bd52987 | 2:38435c74e795 |
|---|---|
| 22 open IN, $homfile || die "File not found - 1\n"; | 22 open IN, $homfile || die "File not found - 1\n"; |
| 23 $lineIn = <IN>; | 23 $lineIn = <IN>; |
| 24 | 24 |
| 25 for($i = 1; $i <= $numOfFiles; $i++){ | 25 for($i = 1; $i <= $numOfFiles; $i++){ |
| 26 print "$i\n"; | 26 print "$i\n"; |
| 27 open FILE, ">".$homfile.".".$i || die "Can't open file"; | 27 open FILE, ">".$homfile.".".$i.".fasta" || die "Can't open file"; |
| 28 print FILE $lineIn; | 28 print FILE $lineIn; |
| 29 $seqs = 1; | 29 $seqs = 1; |
| 30 $lineIn = <IN>; | 30 $lineIn = <IN>; |
| 31 while(defined $lineIn && $seqs < $numPerFile){ | 31 while(defined $lineIn && $seqs < $numPerFile){ |
| 32 print FILE $lineIn; | 32 print FILE $lineIn; |
| 38 $lineIn = <IN>; | 38 $lineIn = <IN>; |
| 39 } | 39 } |
| 40 close FILE; | 40 close FILE; |
| 41 } | 41 } |
| 42 $i = $i -1; | 42 $i = $i -1; |
| 43 open FILE, ">>".$homfile.".".$i; | 43 open FILE, ">>".$homfile.".".$i.".fasta"; |
| 44 while ($lineIn = <IN>){ | 44 while ($lineIn = <IN>){ |
| 45 print FILE $lineIn; | 45 print FILE $lineIn; |
| 46 } | 46 } |
| 47 close FILE; | 47 close FILE; |
| 48 | 48 |