|
awk -F "\t" –v S=”Sample” -v OFS="\t" '{ |
|
if($12=="-" && $25=="-" && $36=="-" && $13==0 && $26==0 && $37==0) print $4,$5,$6,"-",S,$10,$12,$13,$23,$25,$26,$34,$36,$37,$NF; |
|
else if($12=="-" && $25=="-" && $13==0 && $26==0 && $37!=0) print $4,$5,$6,"-",S,$10,$12,$13,$23,$25,$26,$34,$36,$37,$NF; |
|
else if($12=="-" && $36=="-" && $13==0 && $26!=0 && $37==0) print $4,$5,$6,"-",S,$10,$12,$13,$23,$25,$26,$34,$36,$37,$NF; |
|
else if($12=="-" && $13==0 && $26!=0 && $37!=0) print $4,$5,$6,"-",S,$10,$12,$13,$23,$25,$26,$34,$36,$37,$NF; |
|
else if( $25=="-" && $36=="-" && $13!=0 && $26==0 && $37==0) print $4,$5,$6,"-",S,$10,$12,$13,$23,$25,$26,$34,$36,$37,$NF; |
|
else if( $25=="-" && $13!=0 && $26==0 && $37!=0) print $4,$5,$6,"-",S,$10,$12,$13,$23,$25,$26,$34,$36,$37,$NF; |
|
else if( $36=="-" && $13!=0 && $26!=0 && $37==0) print $4,$5,$6,"-",S,$10,$12,$13,$23,$25,$26,$34,$36,$37,$NF; |
|
else print $4,$5,$6,"+",S,$10,$12,$13,$23,$25,$26,$34,$36,$37,$NF; |
|
}' snippet.hmr.bed.5prime.3prime.cloestKnownGFeneTSS1000 > snippet.hmr.overlap.tss.cp |