使用FGAP进行补洞

1. FGAP简介

FGAP利用BLAST将contigs序列比对到基因组草图序列上,寻找重叠到gap区间的最优序列,从而进行补洞。其参考文献:Piro, Vitor C., et al. “FGAP: an automated gap closing tool.” BMC research notes 7.1 (2014): 371.

2. FGAP下载和安装

FGAP官网:http://www.bioinfo.ufpr.br/fgap/

  1. $ wget http://sourceforge.net/projects/fgap/files/MCR_LINUX64b.tar.gz/download
  2. $ tar zxf MCR_LINUX64b.tar.gz
  3. $ cd MCR_LINUX64b
  4. $ ./installMCR.sh /opt/biosoft/MCR
  5.  
  6. $ wget http://sourceforge.net/projects/fgap/files/FGAP_1_7_LINUX64b.tar.gz/download
  7. $ tar zxf FGAP_1_7_LINUX64b.tar.gz -C /opt/biosoft/

2. FGAP的使用

FGAP的简单使用示例:

  1.  
  1. $ ln -s /opt/biosoft/FGAP_1_7_LINUX64b/sample_data/* .
  2. $ /opt/biosoft/FGAP_1_7_LINUX64b/run_fgap.sh /opt/biosoft/MCR/v717/ \
  3. -d DRAFT_ecoli_hiseq454.fasta \
  4. -a "DATASET_ecoli_454.fasta,DATASET_ecoli_hiseq.fasta"
  5. FGAP的使用参数:
  6.  
  7. -d /--draft-file Draft genome file [fasta format - Ex: "draft.fasta"]
  8. -a /--datasets-files List of datasets files to close gaps [fasta format - Ex: "dataset1.fasta,dataset2.fasta"]
  9.  
  10. -s /--min-score Min Score (raw) to return results from BLAST (integer) - Default: 25
  11. -e /--max-evalue Max E-Value to return results from BLAST (float) - Default: 1e-7
  12. -i /--min-identity Min identity (%) to return results from BLAST (integer [0-100]) - Default: 70
  13.  
  14. -C /--contig-end-length Length (bp) of contig ends to perform BLAST alignment (integer) - Default: 300
  15. -T /--edge-trim-length Length of ignored bases (bp) upstream and downstrem of the gap (integer) - Default: 0
  16. -R /--max-remove-length Max number of bases (bp) that can be removed (integer) - Default: 500
  17. -I /--max-insert-length Max number of bases (bp) that can be inserted (integer) - Default: 500
  18.  
  19. -p /--positive-gap Enable closing of positive gaps (with insertion) (integer [0-1]) - Default: 1
  20. -z /--zero-gap Enable closing of zero gaps (without insert any base) (integer [0-1]) - Default: 0
  21. -g /--negative-gap Enable closing of negative gaps (overlapping contig ends) (integer [0-1]) - Default: 0
  22.  
  23. -c /--gap-char Base that represents the gap (char) - Default: "N"
  24. -b /--blast-path Blast+ package path (only makeblastdb and blastn are needed, version 2.2.28+ or higher) - Default: ""
  25. -l /--blast-alignment-parameters BLAST alignment parameters (opengap,extendgap,match,mismatch,wordsize) - Default: "1,1,1,-3,15"
  26. -r /--blast-max-results Max results from BLAST for each query (integer) - Default: 200
  27. -t /--threads Number of threads (integer) - Default: 1
  28.  
  29. -m /--more-output More output files with gap regions after and before gap closing (integer [0-1]) - Default: 0
  30. -o /--output-prefix Output prefix [File or folder - Ex: "out" or "out/" ] - Default: "output_fgap"
  31. -h /--help This help message

原文来自:http://www.chenlianfu.com/?p=2333

 

发表评论

匿名网友

拖动滑块以完成验证
加载失败