Split RNA-seq reads
Hi all,
I have 75bp RNA-seq reads. I would like to split it to 38*37 bp. Is somebody know how is the easiest way to do that? The reads files are in *.gz format and I don't want to unzip it. Thanks Morteza
• 1,638 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
mkh •
•
written
by
mkh •
1 answer
gunzip -c fastq.gz | awk '{if(NR%4==1 || NR%4==3) {print;} else {printf("%s\n",substr($0,1,38));}}' | gzip > 38.fq.gz
gunzip -c fastq.gz | awk '{if(NR%4==1 || NR%4==3) {print;} else {printf("%s\n",substr($0,39));}}' | gzip > 37.fq.gz
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
Split DNA sequences from 150pb to 75 pb
written by Ismael •
-
GRCH38 counts significantly lower than for version 37?
written by ab123 •
-
A pipeline/software to analyse Drip-seq
written by mkh •
-
How to visualize Hi-C bed file format?
written by mkh •
-
chipseq identified sequences
written by mkh •
-
FIMO: scan motif on a DNA sequence
written by mkh •
-
Difference Between Gene Expression And Transcript Expression
written by mkh •
-
metagenomic profiling and population analysis
written by mkh •
-
How to merge the reads into single file
written by Bioblazer
-
Split Fastq Files Into Chunks Of 1M Reads
written by Bioscientist
