in reply to Seeking help with script
use warnings; use strict; my $i; my $scaf; while (<DATA>) { if (/^>(.*)$/) { $i = 0; $scaf = $1; print ">$scaf.$i\n"; } elsif (/^N/) { $i++; print ">$scaf.$i\n"; } else { print; } } __DATA__ >Long_DNA_Scaffold atgctagsctgtagctagctagtcgatgctagtagct ttagctagctgatgctagtcgatgctagtcatctagc NNNNNNNNNNNNNNNNNNN tagctagctagctgatcgtagctagtcgtagctagct ggtcgatgctgatcgtagctgatgctagtcgatgcta >Long_DNA_Scaffold2 atgctagsctgtagctagctagtcgatgctagtagct tagctagctgatcgtagctagctgatgctagctagtc NNNNNNNNNNNNNNNNNN tagctagctagctgatgctagctagtcgatcgtagct ggctagtgctagtgctagctgatgctagtcgtagctg
Prints:
>Long_DNA_Scaffold.0 atgctagsctgtagctagctagtcgatgctagtagct ttagctagctgatgctagtcgatgctagtcatctagc >Long_DNA_Scaffold.1 tagctagctagctgatcgtagctagtcgtagctagct ggtcgatgctgatcgtagctgatgctagtcgatgcta >Long_DNA_Scaffold2.0 atgctagsctgtagctagctagtcgatgctagtagct tagctagctgatcgtagctagctgatgctagctagtc >Long_DNA_Scaffold2.1 tagctagctagctgatgctagctagtcgatcgtagct ggctagtgctagtgctagctgatgctagtcgtagctg
Other comments:
|
|---|
| Replies are listed 'Best First'. | |
|---|---|
|
Re^2: Seeking help with script
by docalf (Initiate) on Jun 29, 2010 at 15:26 UTC |