Entering edit mode
10.8 years ago
bioinforupesh2009.au
▴
140
I would like to have sequence with corresponding quality values according to length distribution of sequences of fastq in as such fastq output. suppose i want to extract that reads which falls sequence length of 16, 17, and 18. please keep in mind that my fastq reads in color space.. exp:
@5521_14_238_F3
T0133103120031002212223
+
?><5=;<<<12>=<;1;;=5);
@5521_16_404_F3
T03310320002130202331112
+
78;:;;><>9=9;<<2=><<1;5
any tricks by one line bash commands or by perl (without using bioperl modules)???
Thanks ashutosh its done !!!