User:Lindenb/Notebook/UMR915/2012/02/22
From OpenWetWare
 Project name
|
<html><img src="/images/9/94/Report.png" border="0" /></html> Main project page <html><img src="/images/c/c3/Resultset_previous.png" border="0" /></html>Previous entry<html> </html>Next entry<html><img src="/images/5/5c/Resultset_next.png" border="0" /></html> |
Daily
VARKIT=${HOME}/src/variationtoolkit/bin
gunzip -c predictions.txt.gz | grep -v "SnpEff" | awk -F' ' '(!($$16=="INTERGENIC" || $$16=="UPSTREAM" || $$16=="DOWNSTREAM"|| $$16=="INTRON" || $16=="UTR_3_PRIME" || $16=="UTR_5_PRIME"))' |\
awk '{printf("%s_%s_%s_%s\t%s\n",$1,$2,$3,$4,$0);}'| sort -t ' ' -k1,1 > _jeter1.txt
gunzip -c all.vcf.gz |\
${VARKIT}/vcf2snpeff |\
awk '{printf("%s_%s_%s_%s\t%s\n",$1,$2,$3,$4,$0);}'| cut -f 1,6-| sort -t ' ' -k1,1 > _jeter2.txt
| |
