Advertisement
ProzacR

Vardict_paleisti

Jun 29th, 2018
559
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.21 KB | None | 0 0
  1. http://genome.ucsc.edu/cgi-bin/hgTables?hgsid=678698205_C0Y8esueodF1NMKYgNEW3NyIUbUo&clade=mammal&org=Human&db=hg19&hgta_group=allTracks&hgta_track=knownGene&hgta_table=0&hgta_regionType=genome&position=chr21%3A33%2C031%2C597-33%2C041%2C570&hgta_outputType=primaryTable&hgta_outFileName=
  2.  
  3. #!/bin/bash
  4. #
  5. #for i in *.bam
  6. #
  7. #do
  8. #
  9. #echo "Indexing: "$i
  10. #
  11. #samtools index $i $i".bai"
  12. #
  13. #done
  14.  
  15. cd /home/farmacija/bin/vardict/VarDict-master
  16. AF_THR="0.01" # minimum allele frequency
  17. #./vardict -G ~/Downloads/hg19/hg19.fa -f $AF_THR -N test -b my.bam -c 1 -S 2 -E 3 -g 4 my.bed | ./teststrandbias.R | ./var2vcf_valid.pl -N test -E -f $AF_THR
  18.  
  19.  
  20. #nohup sh -c "./vardict -G ~/Downloads/hg19/hg19.fa -f $AF_THR -N test -b /home/farmacija/HKMP-dysl/IonCode_0101_R_2018_04_19_05_33_53_user_S5-00727-7-cardio-dysl_Auto_user_S5-00727-7-cardio-dysl_140.bam -c 1 -S 2 -E 3 -g 4 /home/farmacija/Downloads/knownGene.bed | ./teststrandbias.R | ./var2vcf_valid.pl -N test -E -f $AF_THR" &
  21.  
  22. nohup sh -c "./vardict -G ~/Downloads/hg19/hg19.fa -f $AF_THR -N test -b /home/farmacija/HKMP-dysl/IonCode_0102_*.bam -c 1 -S 2 -E 3 -g 4 /home/farmacija/HKMP-dysl/kmp-dysl.bed | ./teststrandbias.R | ./var2vcf_valid.pl -N test -E -f $AF_THR" &
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement