code to remove reads whose QUAL != SEQ length

fix.seq.sh

#!/bin/bash

cat test.fastq | awk 'BEGIN{OFS="\n"} { 
a[NR % 4] = $0; 
if(NR % 4 == 0 && length(a[2]) == length(a[0])){
print a[1],a[2],a[3],a[0] 
}
}' > sim.fastq