Assemble mitochondrial genomes from short read data

From HPCwiki
Revision as of 13:06, 26 March 2014 by Megen002 (talk | contribs) (Created page with "-------------------------------------------------------------- <source lang='bash'> #!/bin/bash #SBATCH --time=10000 #SBATCH --mem=16000 #SBATCH --ntasks=8 #SBATCH --nodes=1 ...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

<source lang='bash'>

  1. !/bin/bash
  2. SBATCH --time=10000
  3. SBATCH --mem=16000
  4. SBATCH --ntasks=8
  5. SBATCH --nodes=1
  6. SBATCH --constraint=normalmem
  7. SBATCH --output=output_pep_%j.txt
  8. SBATCH --error=error_output_pep_%j.txt
  9. SBATCH --job-name=blastp
  10. SBATCH --partition=ABGC_Research
  11. SBATCH --mail-type=ALL
  12. SBATCH --mail-user=hendrik-jan.megens@wur.nl

bowtie2 --phred$2 --local -p 8 -x mt_pig.fa -1 $3 -2 $4 | head -2 >$1_mito_align.sam bowtie2 --phred$2 --local -p 8 -x mt_pig.fa -1 $3 -2 $4 | awk '$5>0' | head -10000 >>$1_mito_align.sam

java7 -jar /cm/shared/apps/SHARED/picard-tools/picard-tools-1.109/SamToFastq.jar I=$1_mito_align.sam F=fq1.fq F2=fq2.fq INCLUDE_NON_PF_READS=True SOAPdenovo-63mer all -K 63 -p 4 -s soapdenovo_pig.config -o $1_mito_assembly.fa blastn -query $1_mito_assembly.fa.scafSeq -db mt_pig.fa -outfmt 6 </source

Retrieved from "https://wiki.anunna.wur.nl/index.php?title=Assemble_mitochondrial_genomes_from_short_read_data&oldid=1286"