I’m trying to run Kneaddata with some pair-end NovaSeq data. I don’t know why I keep getting empty files as outputs. Also, the run will never stop. I’m wondering if there is any issues in my loop. Here is the loop I’m using and attached is the log file.
ls | grep "R1_new.fastq"
for SAMPLE_READ1 in $files_knead_R1;
echo $SAMPLE_READ1 | cut -d R -f1SAMPLE_OUT.fastq
echo $SAMPLE_READ1 | cut -d R -f1R2_new.fastq
kneaddata -t 24 --input $SAMPLE_READ1 --input $SAMPLE_READ2 -db /blue/lorca/m.torrezlamberti/NS2074/RAW2/human_genome --output kneaddata_clean/$SAMPLE_OUTlog file.txt (4.1 KB)
On the other hand, I also tried the following codes
sed '1~4 s///1/g’ < original_R1.fastq > new_R1.fastq
sed '1~4 s///2/g’ < original_R2.fastq > new_R2.fastq
because I thought the issue could be that the R1 and R2 were not able to link. However that did not fix the issue above.
I’ll appreciate any help