in reply to converting fastaq to fasta
Try perl -nle"next if ($.-1)%4 > 1; s[^\@][>]; print" filein.FASTQ > fileout.FASTA like this:
c:\test>type filein.FASTQ @SEQ_ID1 GATTTGGGGTTCAAAGCAGTATCGATCAAATAGTAAATCCATTTGTTCAACTCACAGTTT +SEQ_ID1 !''*((((***+))%%%++)(%%%%).1***-+*''))**55CCF>>>>>>CCCCCCC65 @SEQ_ID2 GATTTGGGGTTCAAAGCAGTATCGATCAAATAGTAAATCCATTTGTTCAACTCACAGTTT +SEQ_ID2 !''*((((***+))%%%++)(%%%%).1***-+*''))**55CCF>>>>>>CCCCCCC65 @SEQ_ID3 GATTTGGGGTTCAAAGCAGTATCGATCAAATAGTAAATCCATTTGTTCAACTCACAGTTT +SEQ_ID3 !''*((((***+))%%%++)(%%%%).1***-+*''))**55CCF>>>>>>CCCCCCC65 @SEQ_ID4 GATTTGGGGTTCAAAGCAGTATCGATCAAATAGTAAATCCATTTGTTCAACTCACAGTTT +SEQ_ID4 !''*((((***+))%%%++)(%%%%).1***-+*''))**55CCF>>>>>>CCCCCCC65 c:\test>perl -nle"next if ($.-1)%4 > 1; s[^\@][>]; print" filein.FASTQ + > fileout.FASTA c:\test>type fileout.FASTA >SEQ_ID1 GATTTGGGGTTCAAAGCAGTATCGATCAAATAGTAAATCCATTTGTTCAACTCACAGTTT >SEQ_ID2 GATTTGGGGTTCAAAGCAGTATCGATCAAATAGTAAATCCATTTGTTCAACTCACAGTTT >SEQ_ID3 GATTTGGGGTTCAAAGCAGTATCGATCAAATAGTAAATCCATTTGTTCAACTCACAGTTT >SEQ_ID4 GATTTGGGGTTCAAAGCAGTATCGATCAAATAGTAAATCCATTTGTTCAACTCACAGTTT
Note:Switch "s to 's if you're on *nix.
|
|---|
| Replies are listed 'Best First'. | |
|---|---|
|
Re^2: converting fastaq to fasta
by Anonymous Monk on Feb 08, 2012 at 07:12 UTC | |
by BrowserUk (Patriarch) on Feb 08, 2012 at 07:25 UTC |