SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   RNA Sequencing (http://seqanswers.com/forums/forumdisplay.php?f=26)
-   -   How to parse FASTQ file? (http://seqanswers.com/forums/showthread.php?t=27063)

Shishir 02-01-2013 04:37 AM

How to parse FASTQ file?
 
I am new in NGS. I need to format some fastq files. For example I want to replace characters /1 and /2 with -1 and - 2 in the files having data like @HWUSI-EAS1671_0001:5:1:1022:10290/1.
I tried to do it in perl with pattern matching but the files are too big to parse. Could anyone please give me suggestion on share any code?

GenoMax 02-01-2013 05:03 AM

Have you considered using "awk"?

Shishir 02-01-2013 05:34 AM

Thanks!
I searched for "awk" and realize "sed" can do it easily :)

kopi-o 02-01-2013 05:42 AM

The files are surely not too big to parse with perl (for your purposes) if you do it right - you only need to have four lines at a time (or maybe just one) in memory. Also it might be a good opportunity to learn some BioPerl :-)


All times are GMT -8. The time now is 09:49 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.