-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathcombineTables_command.sh
executable file
·15 lines (6 loc) · 7.05 KB
/
combineTables_command.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
echo "Region 158_2 158_4 158_5 159_6 159_7 161_4 161_5 163_1 163_2 164_8 164_9 165_7 167_3 167_5 172_2 173_1 838_5 838_6 839_1 840_1 840_3 841_2 844_1 845_1 C57_14 C57_15 C57_27 C57_29 C57_31 C57_32 C57_35 C57_36" > CTR_edw23_0003_MeDIPS_Region_Coverage_Combo_5Kb.table
paste <(awk '{print $1"_"$2"_"$3" "1000000*($4/51147869) }' 158_2.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/29352775)}' 158_4.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/47656057)}' 158_5.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/43970107)}' 159_6.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/55865174)}' 159_7.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/46870534)}' 161_4.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/46537599)}' 161_5.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/45576173)}' 163_1.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/49687171)}' 163_2.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/52341120)}' 164_8.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/47459700)}' 164_9.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/54175204)}' 165_7.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/39018829)}' 167_3.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/37627119)}' 167_5.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/39549149)}' 172_2.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/43457002)}' 173_1.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/46502056)}' 838_5.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/44188002)}' 838_6.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/38974184)}' 839_1.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/41823568)}' 840_1.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/63925901)}' 840_3.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/25819566)}' 841_2.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/40028813)}' 844_1.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/48457425)}' 845_1.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/46857300)}' C57_14.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/37871675)}' C57_15.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/29033776)}' C57_27.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/42387931)}' C57_29.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/47392942)}' C57_31.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/44104825)}' C57_32.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/33604354)}' C57_35.chr.srtd.5Kb.bed) <(awk '{print 1000000*($4/39415987)}' C57_36.chr.srtd.5Kb.bed) | grep -v '^chrM' | column >> CTR_edw23_0003_MeDIPS_Region_Coverage_Combo_5Kb.table
echo "Region 158_2 158_4 158_5 159_6 159_7 161_4 161_5 163_1 163_2 164_8 164_9 165_7 167_3 167_5 172_2 173_1 838_5 838_6 839_1 840_1 840_3 841_2 844_1 845_1 C57_14 C57_15 C57_27 C57_29 C57_31 C57_32 C57_35 C57_36" > CTR_edw23_0003_MeDIPS_Region_Coverage_Combo_25Kb.table
paste <(awk '{print $1"_"$2"_"$3" "1000000*($4/50272322) }' 158_2.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/28810372)}' 158_4.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/46790548)}' 158_5.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/43140479)}' 159_6.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/54859314)}' 159_7.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/46015733)}' 161_4.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/45711768)}' 161_5.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/44752121)}' 163_1.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/48818758)}' 163_2.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/51364916)}' 164_8.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/46571493)}' 164_9.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/53139675)}' 165_7.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/38319249)}' 167_3.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/36954384)}' 167_5.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/38843136)}' 172_2.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/42669423)}' 173_1.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/45680116)}' 838_5.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/43350784)}' 838_6.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/38286077)}' 839_1.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/41052419)}' 840_1.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/62768077)}' 840_3.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/25373431)}' 841_2.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/39337404)}' 844_1.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/47577726)}' 845_1.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/46008168)}' C57_14.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/37209452)}' C57_15.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/28507167)}' C57_27.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/41633077)}' C57_29.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/46528685)}' C57_31.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/43268331)}' C57_32.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/33000160)}' C57_35.chr.srtd.25Kb.bed) <(awk '{print 1000000*($4/38699703)}' C57_36.chr.srtd.25Kb.bed) | grep -v '^chrM' | column >> CTR_edw23_0003_MeDIPS_Region_Coverage_Combo_25Kb.table
echo "Region 158_2 158_4 158_5 159_6 159_7 161_4 161_5 163_1 163_2 164_8 164_9 165_7 167_3 167_5 172_2 173_1 838_5 838_6 839_1 840_1 840_3 841_2 844_1 845_1 C57_14 C57_15 C57_27 C57_29 C57_31 C57_32 C57_35 C57_36" > CTR_edw23_0003_MeDIPS_Region_Coverage_Combo_100Kb.table
paste <(awk '{print $1"_"$2"_"$3" "1000000*($4/51207410) }' 158_2.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/29389131)}' 158_4.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/47816337)}' 158_5.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/44134124)}' 159_6.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/55944880)}' 159_7.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/46912038)}' 161_4.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/46683474)}' 161_5.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/45651502)}' 163_1.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/49786697)}' 163_2.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/52477146)}' 164_8.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/47491500)}' 164_9.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/54282088)}' 165_7.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/39135968)}' 167_3.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/37708993)}' 167_5.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/39668356)}' 172_2.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/43839995)}' 173_1.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/46592414)}' 838_5.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/44236151)}' 838_6.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/39061941)}' 839_1.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/41835829)}' 840_1.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/63976375)}' 840_3.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/25919915)}' 841_2.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/40164082)}' 844_1.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/48498639)}' 845_1.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/46993085)}' C57_14.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/38054347)}' C57_15.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/29156445)}' C57_27.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/42555602)}' C57_29.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/47932679)}' C57_31.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/44263303)}' C57_32.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/33690579)}' C57_35.chr.srtd.100Kb.bed) <(awk '{print 1000000*($4/39649118)}' C57_36.chr.srtd.100Kb.bed) | grep -v '^chrM' | column >> CTR_edw23_0003_MeDIPS_Region_Coverage_Combo_100Kb.table