=====motif反向定位=====

chip项目中，我们给老师做了motif富集，有老师有需求，要找到motif在染色体中位置。

  1./TJPROJ6/RNA_SH/personal_dir/zhangxin/jiaoben/chip/motif/bin/summits_bed_to_get_motif.sh处理峰的文件，作为bed文件
  2./TJPROJ6/RNA_SH/personal_dir/zhangxin/jiaoben/chip/motif/bin/prepare_motif_file.py生成脚本
  
    usage: prepare_motif_file.py [-h] [-I INMOTIF] [-L LIST] [-O OUTDIR] [-F FA]
                             [-B BED] [-S]
    
    prepare motif file
    
    optional arguments:
      -h, --help            show this help message and exit
      -I INMOTIF, --inmotif INMOTIF
                        the matrix of motif, multi file split by ",", such as
                        /PUBLIC/software/RNA/HOMER/data/knownTFs/known.motifs    #手动流程的总背景文件
                        /TJPROJ7/EPI/WORK/Software/miniconda3/envs/motif/share/homer/data/knownTFs/known.motifs   #自动化流程的总背景文件
       -L LIST, --list LIST  the motif name
       -O OUTDIR, --outdir OUTDIR
                        the out dir
       -F FA, --fa FA        the fa file
       -B BED, --bed BED     the bed file
       -S, --same            use same parameter with pipline  添加该参数可保证找到的目标序列位置与富集中的数量一致