User:Lindenb/Notebook/UMR915/2012/02/22: Difference between revisions

From OpenWetWare
< User:Lindenb‎ | Notebook‎ | UMR915‎ | 2012‎ | 02
Jump to navigationJump to search
(Autocreate 2012/02/22 Entry for User:Lindenb/Notebook/UMR915)
 
Line 8: Line 8:
==Daily==
==Daily==
* [[User:lindenb/Notebook/UMR915/20100607|20100607]]
* [[User:lindenb/Notebook/UMR915/20100607|20100607]]
<pre>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</pre>





Revision as of 08:40, 22 February 2012

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>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</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