User Tools

Site Tools


deconseq_for_decontaminating_reads

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
deconseq_for_decontaminating_reads [2017/08/24 12:57] 129.173.90.165deconseq_for_decontaminating_reads [2018/06/28 15:58] (current) 134.190.133.209
Line 1: Line 1:
-**DECONSEQ**+====== DECONSEQ ====== 
  
 Documentation by Sarah Shah, Shelby Williams Documentation by Sarah Shah, Shelby Williams
Line 13: Line 14:
 bwa64 index -p prefix_of_your_choice_for_bacteria_index -a bwtsw fasta_file_of_your_suspected_bacteria >bwa.log 2>&1 & bwa64 index -p prefix_of_your_choice_for_bacteria_index -a bwtsw fasta_file_of_your_suspected_bacteria >bwa.log 2>&1 &
 </code> </code>
 +
 +You can check the status of your bwa index by using "tail -f bwa.log".
  
 2. Move the 8 index files to the folder "db". 2. Move the 8 index files to the folder "db".
  
-3. Edit the config file (DeconSeqConfig.pm) as such:+3a. Edit the config file (DeconSeqConfig.pm) as such:
 <code> <code>
 'prefix_of_bacteria_index' => {name => 'Nameofyourchoice',  'prefix_of_bacteria_index' => {name => 'Nameofyourchoice', 
                         db => 'prefix_of_bacteria_index'},                         db => 'prefix_of_bacteria_index'},
 </code> </code>
-Do not edit anything else.+ 
 +3b. Also edit the following lines in the config file: 
 +<code> 
 +use constant DB_DIR => '/path/to/deconseq/db/'; 
 +use constant PROG_DIR => '/path/to/deconseq/'; 
 +</code> 
 + 
 +This is because of a bug (see https://www.biostars.org/p/182337/). Do not edit anything else.
  
 4. Write a shell script (example below) and qsub it. 4. Write a shell script (example below) and qsub it.
deconseq_for_decontaminating_reads.1503590251.txt.gz · Last modified: by 129.173.90.165