Commit 6c2541b2 authored by Kosmas Hench's avatar Kosmas Hench

update 2.2.9 selscan scripts

parent 2d36059b
......@@ -27,3 +27,15 @@ for k in nig pue uni nigbel puebel unibel nighon puehon unihon nigboc pueboc uni
END{for(i in A){print H[i],(G[i]-1)*slide,(G[i]-1)*slide+window,B[i],A[i]/B[i],G[i],k+1,B[i]/window,C[i]/B[i];k++}}' >> selscan_merge_10kb/$k.selscan.ihh12.10k.1k.smooth.out
done
done
cd selscan_merge_10kb
for k in pue-nig pue-uni nig-uni puebel-nigbel puebel-unibel nigbel-unibel puehon-nighon puehon-unihon nighon-unihon pueboc-nigboc pueboc-uniboc nigboc-uniboc; do
echo $k
for j in {01..24};do
LG="LG"$j;
echo $LG;
mkdir -p $LG;
head -n 1 $k.selscan.ihh12.10k.1k.smooth.out > $LG/iHH12.$k.pue-nig.selscan-$LG-10kb-1kb.txt
done
done
\ No newline at end of file
#PBS -l elapstim_req=00:20:00
#PBS -l elapstim_req=01:20:00
#PBS -l memsz_job=80gb
#PBS -b 1
#PBS -l cpunum_job=1
......@@ -25,3 +25,15 @@ for k in pue-nig pue-uni nig-uni puebel-nigbel puebel-unibel nigbel-unibel pueho
END{for(i in A){print H[i],(G[i]-1)*slide,(G[i]-1)*slide+window,B[i],A[i]/B[i],G[i],k+1,B[i]/window,C[i]/B[i];k++}}' >> selscan_merge_10kb/$k.selscan.xpEHH.10k.1k.smooth.out
done
done
cd selscan_merge_10kb
for k in pue-nig pue-uni nig-uni puebel-nigbel puebel-unibel nigbel-unibel puehon-nighon puehon-unihon nighon-unihon pueboc-nigboc pueboc-uniboc nigboc-uniboc; do
echo $k
for j in {01..24};do
LG="LG"$j;
echo $LG;
mkdir -p $LG;
head -n 1 $k.selscan.xpEHH.10k.1k.smooth.out > $LG/xpEHH.$k.pue-nig-$LG-10kb-1kb.txt
done
done
\ No newline at end of file
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment