This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Replace header line of fastq file with third line

Hi,

I want to change my fastq file

@1/1

GTCGCTGACGAGCGAGTACT

+SRR00000.1 1 length=150

FFFFFFFFFFFFFFFFFFFFFFFF

into something like

@SRR00000.1 1 length=150

GTCGCTGACGAGCGAGTACT

+SRR00000.1 1 length=150

FFFFFFFFFFFFFFFFFFFFFFFF

Replacing header line with third line. Iam trying something like this awk 'NR%4==1 {$0=$3} {print}' but ofcourse $3 would be empty in the start so it gives empty header line.

Any help with awk command or anything else would be really appreciated?

Thanks

bash fastq
• 1,693 views • link
updated by cpad0112 cpad0112 • written by madiha_shabbir46 • madiha_shabbir46 
$ zcat test.fq.gz| sed '0~4 s/$/\n/'  | awk  -F "\n" -v OFS="\n" -v RS="\n\n" '{$1=$3; sub("+","@",$1)}1'
• 0 views • link
written by cpad0112 cpad0112

1 answer 

 gunzip -c in.fq.gz | paste - - - - | awk -F '\t'  '{OFS="\n";$3=sprintf("+%s",substr($1,2));print;}'
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thankyou, can you please explain this? specially | paste - - - - | part? Plus we also need to replace + with @ here

• 0 views • link
written by madiha_shabbir46 • madiha_shabbir46

Thank you so much. With little update in line numbers, it's perfect.

gunzip -c in.fq.gz | paste - - - - | awk -F '\t' '{OFS="\n";$1=sprintf("@%s",substr($3,2));print;}'

• 0 views • link
written by madiha_shabbir46 • madiha_shabbir46

Log in to answer this question.

More posts like this