Wednesday 30 July 2014

Convert all fastq files in folder to fasta

for fq in *; do fn=$(echo $fq|sed s'/.fq//'); echo $fn;awk 'NR % 4 == 1 {print ">" $0 } NR % 4 == 2 {print $0}' $fq > ../MergedFasta/$fn.fa; done;

No comments:

Post a Comment

Note: only a member of this blog may post a comment.