Entering edit mode
7.8 years ago
Abdul Rafay Khan
★
1.2k
I try Edena but i can't get the same length of my reads, If you can help me I appreciated. ... [err] All reads within a file must be the same length.
edena -singleEnd FileName_1.fastq
Rapid file(s) examination... 42 93
[err] All reads within a file must be the same length
you can get the shortest read in your fasta/q file then trim the rest to be the same length ; but I do not understand why the software behaves in this way, what is the idea of all reads have the same length?
or you can first but a threshold to get all reads greater than length x (for example 100) them trim them all to be 100
this will give you the length of the least reads length
cat file.fastq | awk '{if(NR%2==0) print( length $0);}' | uniq | sort -n | head -n 10
What is the least reads length in above result?
sorry the command should be like
awk '{if(NR%2==0) print( length );}' | uniq| sort | head -n 10
and result will be something like50 51 52 55 ....