fixed a bug that could lead to incorrect read numbers

ShuXingYu94 · Jan 10, 2023 · 320812f · 320812f
1 parent 32d5d4c
commit 320812f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/statistics.sh b/statistics.sh
@@ -77,7 +77,7 @@ do
   # Reads Count (Raw Data) - with seqkit stats
   cd ${basecall_dir}
   fn1=`ls ${basecall_dir} | grep *$file*R1*.gz`
-  fn2=`ls ${basecall_dir} | grep *$file*R1*.gz`
+  fn2=`ls ${basecall_dir} | grep *$file*R2*.gz`
   R1_reads=`seqkit stats $fn1 -T -j ${threads} | awk 'NR>1{print $4}'`
   R2_reads=`seqkit stats $fn2 -T -j ${threads} | awk 'NR>1{print $4}'`