Commit 23d623fc authored by Chaos's avatar Chaos
Browse files

fix a bug in num1 calculation

parent ab6a9d0d

hic_local.sh

+5 −5
Original line number Diff line number Diff line
@@ -38,10 +38,10 @@ do 


	if [ ! -e $res_file ];then

		num1=$(\

			paste

			<(gunzip -c $read1_file)

			<(gunzip -c $read2_file) |

			$awk_cmd '!((NR+2)%4)' |

			paste \

			<(gunzip -c $read1_file) \

			<(gunzip -c $read2_file) | \

			$awk_cmd '!((NR+2)%4)' | \

			grep -cE $ligation \

		)

CRA Git | Maintained and supported by SUSTech CRA and CCSE