How to install DAMIAN?

gene_x 0 like s 519 view s

Tags: software

https://sourceforge.net/projects/damian-pd/

1, install PostgreSQL and gem on 18.04

  1. sudo apt-get update
  2. sudo apt install ruby-dev libffi-dev build-essential
  3. sudo apt-get install postgresql postgresql-contrib
  4. sudo apt-get install libpq-dev
  5. sudo apt install default-jre
  6. sudo apt install hmmer
  7. #sudo apt-get install pgadmin3
  9. sudo gem install pg -v 0.19
  10. sudo gem install axlsx
  11. sudo gem install amatch
  13. #interactive: sudo -u postgres createuser --interactive
  14. #not_interactive: https://medium.com/coding-blocks/creating-user-database-and-adding-access-on-postgresql-8bfcd2f4a91e
  15. #sudo -u postgres psql
  16. #postgres=# create database mydb;
  17. #postgres=# create user myuser with encrypted password 'mypass';
  18. #postgres=# grant all privileges on database mydb to myuser;
  20. sudo -u postgres psql
  21. CREATE USER damian_user WITH PASSWORD 'hamburg_uke';
  22. CREATE DATABASE damian_db WITH OWNER damian_user;
  23. postgre=# \q

2, install blast, tax and pfam

  1. cd databases;
  2. ./get_all.sh;
  3. cd tax; ./get_tax.sh;
  4. cd pfam; ./get_pfam.sh;
  6. #Taxonomy
  7. #The following taxonomy files are required:
  8. #ftp://ftp.ncbi.nlm.nih.gov/pub/taxonomy/taxdump.tar.gz (the downloaded file must unpacked using tar as well as decompressed.)
  9. #http://s3.amazonaws.com/matrixsciencemisc/prot.av2taxid.gz
  10. #
  11. #curl ftp://ftp.ncbi.nih.gov/pub/taxonomy/taxdump.tar.gz -O
  12. #curl "ftp://ftp.ebi.ac.uk/pub/databases/Pfam/releases/Pfam32.0/Pfam-A.hmm.gz" | gunzip > Pfam-A.hmm.txt
  13. #curl "ftp://ftp.ebi.ac.uk/pub/databases/Pfam/releases/Pfam32.0/database_files/pfamA_tax_depth.txt.gz" | gunzip > pfamA_tax_depth.txt
  15. #change settings in config.rb
  16. DB_NAME = 'damian_db'
  17. DB_USER = 'damian_user'
  18. DB_PASS = 'hamburg_uke'
  19. ./damian_database.rb  --erase_and_rebuild --names databases/tax/names.dmp --nodes databases/tax/nodes.dmp --hmm databases/pfam/Pfam-A.hmm.txt --taxdepth databases/pfam/pfamA_tax_depth.txt
  21. #### download and update the blast-database ####
  22. cd /mnt/nvme0n1p1/REFs/blast/
  23. #wget ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nt.gz
  24. #wget ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nr.gz
  25. #perl update_blastdb.pl --decompress nt
  26. #perl update_blastdb.pl --decompress nr
  27. ##https://www.ncbi.nlm.nih.gov/books/NBK569850/
  28. #update_blastdb.pl --decompress nt
  29. #update_blastdb.pl --decompress nr
  30. #curl ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nr.gz -O
  31. #curl ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nt.gz -O
  32. ##makeblastdb -in exons_for_blastall.fasta -input_type fasta -dbtype nucl -title exons_for_blastall -parse_seqids -out exons_for_blastall
  33. #makeblastdb -in nt -out nt -parse_seqids -dbtype nucl
  34. #makeblastdb -in nr -out nr -parse_seqids -dbtype prot
  35. ##curl ftp://ftp.ncbi.nih.gov/pub/taxonomy/taxdump.tar.gz -O
  36. ##curl "ftp://ftp.ebi.ac.uk/pub/databases/Pfam/releases/Pfam32.0/Pfam-A.hmm.gz" | gunzip > Pfam-A.hmm.txt
  37. ##curl "ftp://ftp.ebi.ac.uk/pub/databases/Pfam/releases/Pfam32.0/database_files/pfamA_tax_depth.txt.gz" | gunzip > pfamA_tax_depth.txt
  38. ##or:
  39. ##NO_THIS_SCRIPT: ./get_blast.sh
  41. #Standard databases (nr etc.): rRNA/ITS databases Genomic + transcript databases Betacoronavirus
  42. #curl ftp://ftp.ncbi.nlm.nih.gov/blast/db/Betacoronavirus.tar.gz -O
  43. #The contents are the same between https://ftp.ncbi.nlm.nih.gov/blast/db/ and https://ftp.ncbi.nlm.nih.gov/blast/db/v5/ since v5 is the default!
  44. #https://ftp.ncbi.nlm.nih.gov/blast/db/nt-nucl-metadata.json    #158
  45. for no in 000 001 002 003 004 005 006 007 008 009  010 011 012 013 014 015 016 017 018 019  020 021 022 023 024 025 026 027 028 029  030 031 032 033 034 035 036 037 038 039  040 041 042 043 044 045 046 047 048 049  050 051 052 053 054 055 056 057 058 059  060 061 062 063 064 065 066 067 068 069  070 071 072 073 074 075 076 077 078 079  080 081 082 083 084 085 086 087 088 089  090 091 092 093 094 095 096 097 098 099  100 101 102 103 104 105 106 107 108 109  110 111 112 113 114 115 116 117 118 119  120 121 122 123 124 125 126 127 128 129  130 131 132 133 134 135 136 137 138 139  140 141 142 143 144 145 146 147 148 149  150 151 152 153 154 155 156 157; do
  46.   curl ftp://ftp.ncbi.nlm.nih.gov/blast/db/nt.${no}.tar.gz -O
  47. done
  48. #https://ftp.ncbi.nlm.nih.gov/blast/db/nr-prot-metadata.json    #103
  49. for no in 00 01 02 03 04 05 06 07 08 09  10 11 12 13 14 15 16 17 18 19  20 21 22 23 24 25 26 27 28 29  30 31 32 33 34 35 36 37 38 39  40 41 42 43 44 45 46 47 48 49  50 51 52 53 54 55 56 57 58 59  60 61 62 63 64 65 66 67 68 69  70 71 72 73 74 75 76 77 78 79  80 81 82 83 84 85 86 87 88 89  90 91 92 93 94 95 96 97 98 99  100 101 102; do
  50.   curl ftp://ftp.ncbi.nlm.nih.gov/blast/db/nr.${no}.tar.gz -O
  51. done
  52. tar xzf *tar.gz
  53. ##Sind die RNA-data (Transcriptome)? by default is RNA-data
  54. #damian_database.rb  --erase_and_rebuild --names blast_2020_install/taxdump/names.dmp --nodes blast_2020_install/taxdump/nodes.dmp --hmm pfam/Pfam-A.hmm.txt --taxdepth pfam/#pfamA_tax_depth.txt #pfam annotation cannot be updated!!

3, create .ncbirc and setting

  1. #in the file /home/jhuang/.ncbirc
  2. BLASTDB=/mnt/nvme0n1p1/REFs/blast/
  3. #echo "[BLAST]" > /home/jhuang/.ncbirc
  4. #echo "BLASTDB=/media/jhuang/Elements1/BLAST_db_v5/nt_v5/" >> /home/jhuang/.ncbirc
  6. #mv damian_release damian
  7. # add damian into PATH
  8. DAMIAN_LOCATION='/home/jhuang/Tools/damian'
  9. export PATH=$PATH:$DAMIAN_LOCATION

4, generate bowtie2 index and set damian_reference

  1. ##human
  2. ##Using existing index /ref/Homo_sapiens/UCSC/hg38/Sequence/WholeGenomeFasta/genome.
  3. #
  4. ##Horse (equCab2)
  5. #rsync -a -P rsync://hgdownload.soe.ucsc.edu/goldenPath/equCab2/bigZips/chromFa.tar.gz ./
  6. ##Cattle  NCBI Genome ID: 82 (Bos taurus)
  7. #rsync -a -P rsync://hgdownload.soe.ucsc.edu/goldenPath/bosTau8/bigZips/bosTau8.fa.gz ./
  8. ##ftp://ftp.ensembl.org/pub/release-95/fasta/bos_taurus/dna/
  9. #
  10. ##Sheep NCBI Genome ID: 83 (Ovis aries)
  11. #rsync -a -P rsync://hgdownload.soe.ucsc.edu/goldenPath/oviAri4/bigZips/oviAri4.fa.gz ./
  12. #
  13. ##Wild boar  NCBI Genome ID: 84 (Sus scrofa)
  14. #rsync -a -P rsync://hgdownload.soe.ucsc.edu/goldenPath/susScr11/bigZips/susScr11.fa.gz ./
  15. #
  16. ##salmon salar
  17. #https://www.ncbi.nlm.nih.gov/genome/369?genome_assembly_id=248466
  18. #https://www.ncbi.nlm.nih.gov/genome/?term=salmo%20salar
  19. #https://www.ncbi.nlm.nih.gov/assembly/?term=salmon+salar
  20. #rsync -avz /ref/Homo_sapiens/UCSC/hg38/Sequence/WholeGenomeFasta/genome.fa jhuang@10.162.6.119:/home/jhuang/DATA/
  21. #rsync -a -P salmon_salar_assemblies.tar jhuang@10.162.6.119:/home/jhuang/REFs
  22. #
  23. ##Mosquitoes/culex pipiens
  24. #https://www.ncbi.nlm.nih.gov/Taxonomy/Browser/wwwtax.cgi?lvl=0&id=233155
  25. #https://www.biorxiv.org/content/10.1101/240747v1.full
  26. #aedes mascarensis
  27. #Aedes albopictus
  28. #
  29. #Taxonomy ID: 7176 (Culex quinquefasciatus (southern house mosquito))
  30. #https://www.ncbi.nlm.nih.gov/nuccore/?term=C.+pipiens
  31. #https://www.ncbi.nlm.nih.gov/assembly?LinkName=bioproject_assembly_all&from_uid=18751
  32. #https://www.ncbi.nlm.nih.gov/genome/?term=txid7176[orgn]
  33. #https://www.ncbi.nlm.nih.gov/assembly/GCF_000208785.1/
  34. #https://www.ncbi.nlm.nih.gov/genome/?term=txid263438[orgn]
  35. #
  36. #Taxonomy ID: 7175 (C. pipiens) --> no genome
  37. #https://www.ncbi.nlm.nih.gov/genome/?term=txid7175[orgn]
  38. #https://www.ncbi.nlm.nih.gov/assembly/GCF_000209185.1
  39. #rsync -a -P GCF_000209185_1_CulPip1_0_genomic.fna.gz jhuang@10.162.6.119:/home/jhuang/REFs
  41. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Ovis_aries.Oar_v3.1.dna.toplevel.fa .
  42. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Ovis_aries.Oar_v3.1.cdna.all.fa .
  43. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Ovis_aries.Oar_v3.1.ncrna.fa .
  44. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Sus_scrofa.Sscrofa11.1.dna.toplevel.fa .
  45. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Sus_scrofa.Sscrofa11.1.cdna.all.fa .
  46. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Sus_scrofa.Sscrofa11.1.ncrna.fa .
  47. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Bos_taurus.ARS-UCD1.2.dna.toplevel.fa .
  48. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Bos_taurus.ARS-UCD1.2.cdna.all.fa .
  49. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Bos_taurus.ARS-UCD1.2.ncrna.fa .
  50. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Equus_caballus.EquCab3.0.dna.toplevel.fa .
  51. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Equus_caballus.EquCab3.0.cdna.all.fa .
  52. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Equus_caballus.EquCab3.0.ncrna.fa .
  53. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Salmo_salar.fa .
  54. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/GCF_000209185_1_CulPip1_0_genomic.fa .
  55. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Macaca_mulatta.Mmul_8.0.1.dna.toplevel.fa .
  56. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Macaca_mulatta.Mmul_8.0.1.cdna.all.fa .
  57. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Macaca_mulatta.Mmul_8.0.1.ncrna.fa .
  58. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Ovis_aries_musimon.fa .
  59. rsync -a -P jhuang@10.162.6.119:/home/jhuang/REFs/Cervus_elaphus_hippelaphus.fa  .
  61. #damian_reference.rb --add --host hg38 --type both --fasta /mnt/h/jhuang/ref/Homo_sapiens/UCSC/hg38/Sequence/WholeGenomeFasta/genome.fa --primary --description 'Homo_sapiens_UCSC_hg38 (dna)'
  62. #damian_reference.rb --add --host wildboar --type both --fasta /home/jhuang/REFs/susScr11.fa --primary --description 'Wild boar  NCBI Genome ID: 84 (Sus scrofa) (dna)'
  63. #damian_reference.rb --add --host horse --type both --fasta /home/jhuang/REFs/equCab2.fa --primary --description 'Horse equCab2 (dna)'
  64. #damian_reference.rb --add --host salmon --type both --fasta /home/jhuang/REFs/salmon_salar.fa --primary --description 'Salmon salar RefSeq assembly accession: GCF_000233375.1 (dna)'
  65. ##damian_reference.rb --add --host sheep --type both --fasta /home/jhuang/REFs/oviAri4.fa --primary --description 'Sheep NCBI Genome ID: 83 (Ovis aries) (dna)'
  66. ##damian_reference.rb --add --host cattle --type both --fasta /home/jhuang/REFs/bosTau8.fa --primary --description 'Cattle  NCBI Genome ID: 82 (Bos taurus) (dna)'
  67. ##damian_reference.rb --add --host mosquito --type both --fasta /home/jhuang/REFs/GCF_000209185_1_CulPip1_0_genomic.fa --primary --description 'Culex pipiens quinquefasciatus (dna)'
  69. # -- host index anlegen with ensemble-files --
  70. ftp://ftp.ensembl.org/pub/release-95/fasta/ovis_aries/dna/
  71. #human and human3
  72. damian_reference.rb --add  --host human --type both --fasta ./Homo_sapiens.GRCh38.dna.toplevel.fa --primary --description 'Homo sapiens (dna)'
  73. damian_reference.rb --add  --host human --type rna --fasta ./Homo_sapiens.GRCh38.cdna.all.fa --description 'Homo sapiens (cdna)'
  74. damian_reference.rb --add  --host human --type rna --fasta ./Homo_sapiens.GRCh38.ncrna.fa --description 'Homo sapiens (ncrna)'
  75. #human3 (since for some fastqs, human delete too much and too strictly, therefore we genertate human3 for loose filtering of human reads.
  76. damian_reference.rb --add  --host human3 --type both --fasta ./genome.fa --primary --description 'Homo_sapiens_UCSC_hg38 (dna)'
  77. damian_reference.rb --add  --host human3 --type rna --fasta ./Homo_sapiens.GRCh38.cdna.all.fa --description 'Homo sapiens (cdna)'
  78. damian_reference.rb --add  --host human3 --type rna --fasta ./Homo_sapiens.GRCh38.ncrna.fa --description 'Homo sapiens (ncrna)'
  80. #sheep
  81. damian_reference.rb --add  --host sheep --type both --fasta Ovis_aries.Oar_v3.1.dna.toplevel.fa --primary --description 'Ovis aries (dna)'
  82. damian_reference.rb --add  --host sheep --type rna --fasta Ovis_aries.Oar_v3.1.cdna.all.fa --description 'Ovis aries (cdna)'
  83. damian_reference.rb --add  --host sheep --type rna --fasta Ovis_aries.Oar_v3.1.ncrna.fa --description 'Ovis aries (ncrna)'
  84. #pig
  85. damian_reference.rb --add  --host pig --type both --fasta Sus_scrofa.Sscrofa11.1.dna.toplevel.fa --primary --description 'Sus scrofa (dna)'
  86. damian_reference.rb --add  --host pig --type rna --fasta Sus_scrofa.Sscrofa11.1.cdna.all.fa --description 'Sus scrofa (cdna)'
  87. damian_reference.rb --add  --host pig --type rna --fasta Sus_scrofa.Sscrofa11.1.ncrna.fa --description 'Sus scrofa (ncrna)'
  88. #cow
  89. damian_reference.rb --add  --host cow --type both --fasta Bos_taurus.ARS-UCD1.2.dna.toplevel.fa --primary --description 'Bos taurus (dna)'
  90. damian_reference.rb --add  --host cow --type rna --fasta Bos_taurus.ARS-UCD1.2.cdna.all.fa --description 'Bos taurus (cdna)'
  91. damian_reference.rb --add  --host cow --type rna --fasta Bos_taurus.ARS-UCD1.2.ncrna.fa --description 'Bos taurus (ncrna)'
  93. #horse
  94. damian_reference.rb --add  --host horse --type both --fasta ./Equus_caballus.EquCab3.0.dna.toplevel.fa --primary --description 'Equus caballus (dna)'
  95. damian_reference.rb --add  --host horse --type rna --fasta ./Equus_caballus.EquCab3.0.cdna.all.fa --description 'Equus caballus (cdna)'
  96. damian_reference.rb --add  --host horse --type rna --fasta ./Equus_caballus.EquCab3.0.ncrna.fa --description 'Equus caballus (ncrna)'
  97. #salmo
  98. damian_reference.rb --add  --host Salmo_salar --type both --fasta Salmo_salar.fa --primary --description 'Salmo salar (dna)'
  99. #mosquito
  100. damian_reference.rb --add  --host Culex_pipiens --type both --fasta GCF_000209185_1_CulPip1_0_genomic.fa --primary --description 'Culex pipiens (dna)'
  101. #macaque
  102. damian_reference.rb --add  --host macaque --type both --fasta ./Macaca_mulatta.Mmul_8.0.1.dna.toplevel.fa --primary --description 'Macaca mulatta (dna)'
  103. damian_reference.rb --add  --host macaque --type rna --fasta ./Macaca_mulatta.Mmul_8.0.1.cdna.all.fa --description 'Macaca mulatta (cdna)'
  104. damian_reference.rb --add  --host macaque --type rna --fasta ./Macaca_mulatta.Mmul_8.0.1.ncrna.fa --description 'Macaca mulatta (ncrna)'
  106. #mouflon
  107. damian_reference.rb --add  --host mouflon --type both --fasta ./Ovis_aries_musimon.fa --primary --description 'Ovis aries musimon (dna)'
  109. #reddeer
  110. damian_reference.rb --add  --host reddeer --type both --fasta ./Cervus_elaphus_hippelaphus.fa --primary --description 'Cervus elaphus hippelaphus (dna)'
  112. ##icebear
  113. #damian_reference.rb --add  --host polarbear --type both --fasta ./Ursus_maritimus.UrsMar_1.0.dna.toplevel.fa --primary --description 'Ursus_maritimus (dna)'
  115. ##Der Graue Mausmaki (Microcebus murinus) ist eine Primatenart aus der Gattung der Mausmakis innerhalb der Gruppe der Lemuren.
  116. #damian_reference.rb --add  --host lemur --type both --fasta ./Mmur3.0.fa --primary --description 'Microcebus murinus (dna)'

5, install and configure mutt

  1. sudo apt install mutt
  3. #in ~/.muttrc
  4. set imap_user = 'xxx@yyy.com'
  5. set imap_pass = 'xxxx'
  6. set from= $imap_user
  7. set use_from=yes
  8. set realname='XXX YYY'
  9. set folder = imaps://imap-mail.outlook.com:993
  10. set spoolfile = "+INBOX"
  11. set postponed="+[hotmail]/Drafts"
  12. set mail_check = 100
  13. set header_cache = "~/.mutt/cache/headers"
  14. set message_cachedir = "~/.mutt/cache/bodies"
  15. set certificate_file = "~/.mutt/certificates"
  16. set smtp_url = "smtp://$imap_user@smtp-mail.outlook.com:587"
  17. set smtp_pass = $imap_pass
  18. set move = no
  19. set imap_keepalive = 900
  20. set record="+Sent"
  22. Test: echo -e "Hi XXX,\n\nPlease find attached the latest results from our DAMIAN analysis.\n\nBest,\nYYY" | mutt -s "New results from DAMIAN" -- "xxx@googlemail.com"

6, intermediate commands

  1. --1-- hmmsearch --domE 0.00001 -o /dev/null --domtblout /home/jhuang/rtpd_files/HD04_cons/idba_ud_assembly/domain.table --noali --cpu 10 /home/jhuang/Tools/damian/databases/pfam/Pfam-A.hmm.txt /home/jhuang/rtpd_files/HD04_cons/idba_ud_assembly/orfs.fasta
  2. --2-- megablast
  3. --3-- blastn or blastp
  4. /home/jhuang/Tools/damian/3rd_party/ncbi-blast/bin/blastp -task blastp -evalue 10E-2 -num_threads 26 -query /tmp/rtpd__565_20190514-28525-1cqkejq -db nr -outfmt 6 qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore staxids qcovs qcovhsp
  5. /home/jhuang/Tools/damian/3rd_party/ncbi-blast/bin/blastp -task blastp -evalue 10E-2 -num_threads 10 -query /tmp/rtpd__584_20190515-11072-i8ct4h -db nr -outfmt 6 qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore staxids qcovs qcovhsp
  6. /home/jhuang/Tools/damian/3rd_party/ncbi-blast/bin/blastn -task blastn -evalue 10E-2 -num_threads 10 -query /tmp/rtpd__586_20190515-6605-1wfobqe -db nt -outfmt 6 qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore staxids qcovs qcovhsp

7, Verifying the installation

  1. #damian.rb --left selftest/r1.fastq.gz --right selftest/r2.fastq.gz --sample testrun --threads 12
  3. seqtk sample -s100 ./240621_M03701_0312_000000000-GHL9N/p20534/7448_7501_S0_R1_001.fastq.gz 0.1 > R1_0.1.fastq
  4. seqtk sample -s100 ./240621_M03701_0312_000000000-GHL9N/p20534/7448_7501_S0_R2_001.fastq.gz 0.1 > R2_0.1.fastq
  6. cd /mnt/nvme0n1p1/REFs/blast
  7. damian.rb --host human3 --type rna -1 R1_0.1.fastq -2 R2_0.1.fastq --sample p20534_7448_7501_S0_megablast --blastn never --blastp never --min_contiglength 500 --threads 64 --force
  8. damian_report.rb
  9. zip -r p20534_7448_7501_S0_megablast.zip p20534_7448_7501_S0_megablast/
  10. echo -e "Hi XXX,\n\nPlease find attached the latest results from our DAMIAN analysis.\n\nBest,\nYYY" | mutt -a "./p20534_7448_7501_S0_megablast.zip" -s "New results from DAMIAN" -- "xxx@googlemail.com"
  11. damian.rb --host human3 --type rna -1 R1_0.1.fastq -2 R2_0.1.fastq --sample p20534_7448_7501_S0_blastn --blastn progressive --blastp never --min_contiglength 500 --threads 64 --force
  12. damian_report.rb
  13. zip -r p20534_7448_7501_S0_blastn.zip p20534_7448_7501_S0_blastn/
  14. echo -e "Hi XXX,\n\nPlease find attached the latest results from our DAMIAN analysis.\n\nBest,\nYYY" | mutt -a "./p20534_7448_7501_S0_blastn.zip" -s "New results from DAMIAN" -- "xxx@googlemail.com"
  15. damian.rb --host human3 --type rna -1 R1_0.1.fastq -2 R2_0.1.fastq --sample p20534_7448_7501_S0_blastp --blastn never --blastp progressive --min_contiglength 500 --threads 64 --force
  16. damian_report.rb
  17. zip -r p20534_7448_7501_S0_blastp.zip p20534_7448_7501_S0_blastp/
  18. echo -e "Hi XXX,\n\nPlease find attached the latest results from our DAMIAN analysis.\n\nBest,\nYYY" | mutt -a "./p20534_7448_7501_S0_blastp.zip" -s "New results from DAMIAN" -- "xxx@googlemail.com"

like unlike

点赞本文的读者

还没有人对此文章表态

本文有评论

没有评论

看文章，发评论，不要沉默

最受欢迎文章

  1. Why Do Significant Gene Lists Change After Adding Additional Conditions in Differential Gene Expression Analysis?
  2. Motif Discovery in Biological Sequences: A Comparison of MEME and HOMER
  3. Calling peaks using findPeaks of HOMER
  4. Updating Human Gene Identifiers using Ensembl BioMart: A Step-by-Step Guide
  5. Setup conda environments
  6. pheatmap vs heatmap.2
  7. Should the inputs for GSVA be normalized or raw?
  8. File format for single channel analysis of Agilent microarray data with Limma?
  9. Kraken2 Installation and Usage Guide
  10. 洛那法尼(lonafarnib)：从抗癌症到抗病毒的多功能药物

最新文章

最多评论文章

推荐相似文章

Essential Open Source Software for Science (EOSS)

Easyfig usage

Utilizing Random Forest to Differentiate Bacterial vs. Viral Infections via Host Gene Expression Signatures

Should the inputs for GSVA be normalized or raw?

© 2023 XGenes.com Impressum