INTELLIGENT WORK FORUMS
FOR COMPUTER PROFESSIONALS

Log In

Come Join Us!

Are you a
Computer / IT professional?
Join Tek-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!

*Tek-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Jobs

AWK for file .fastq

AWK for file .fastq

(OP)
Hi folks!
this is a .fastq file example:

@IRIS:7:1:17:800#0/1
GGAAACACTACTTAGGCTTATAAGATCNGGTTGCGG
+IRIS:7:1:17:800#0/1
ababbaaabaaaaa`]`ba`]`aaaaYD\\_a``XT
@IRIS:7:1:17:150#0/1
TGATGTACTATGCATATGAACTTGTATGCAAAGTGG
+IRIS:7:1:17:150#0/1
abaabaa`aaaaaaa^ba_]]aaa^aaaaa_^][aa
@IRIS:7:1:18:16#0/1
AGGTTCGTGTTGAGTGTTGCCTCTTTTTCTGTTAAT
+IRIS:7:1:18:16#0/1
abb`bbbabbbaWa`baba_Z``babaa[_aa`\X_
@IRIS:7:1:19:1343#0/1
GGCATCTCCAGAGGAGGCTGTACCTGTGGAATAGCA
+IRIS:7:1:19:1343#0/1
abbbbbbbbababaaba_a\FXQXNVKO^F[RWYSQ

and so on for thousand of record
Each record is made of four line i.e

@IRIS:7:1:17:800#0/1
GGAAACACTACTTAGGCTTATAAGATCNGGTTGCGG
+IRIS:7:1:17:800#0/1
ababbaaabaaaaa`]`ba`]`aaaaYD\\_a``XT

Now,
I was wondering if anyone knows how to get through AWK or SED, RANDOM records from a fastq file.
thanks!

RE: AWK for file .fastq

Hi

CODE --> command-line

gawk 'ENDFILE{srand();print r=int(NR/4*rand())*4}FNR!=NR&&FNR>r&&FNR<=r+4' .fastq .fastq 
( Note the input file name specified twice. )

Tested with GNU awk, requires version 4.0.0. ( If necessary, is possible to replace the ENDFILE pattern with a workaround. )

Feherke.
feherke.ga

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Tek-Tips Forums free from inappropriate posts.
The Tek-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Tek-Tips forums is a member-only feature.

Click Here to join Tek-Tips and talk with other members!

Resources

Close Box

Join Tek-Tips® Today!

Join your peers on the Internet's largest technical computer professional community.
It's easy to join and it's free.

Here's Why Members Love Tek-Tips Forums:

Register now while it's still free!

Already a member? Close this window and log in.

Join Us             Close