Tips for routine tools

1.hmmer

convalign.py stockholm file.aln(file.fasta)

hmmbuild –amino file.hmm file.sth

hmmsearch -o file.hmm.hit –noali –nobias file.hmm targetPep.fasta

 

 

2.blast

短序列
psiblast -query AT4G13850.1 -db ../../../conservation/subject/AT4G13850.1 -out AT4G13850.1.out -num_iterations 3 -evalue 20000 -matrix PAM30 -comp_based_stats 0

psiblast -query $(basename $@) -db /data7T/mercu-b-backup/pu b-data/NCBI-NR/nr -out $@ -outfmt 0 -num_threads 8 -threshold 0.005 -eval ue 10 -max_target_seqs 1000 -num_iterations 10 -inclusion_ethresh 0.002

#file.fasta由clustralw等多序列比对软件生成

psiblast -in_msa file.fasta ../../../conservation/subject/AT4G13850.1 -out AT4G13850.1.out -num_iterations 3 -evalue 20000 -matrix PAM30 -comp_based_stats 0

 

#返回fasta格式的数据

blastdbcmd -entry NP_172083.1 -db /data7T/mercu-b-backup/pub-data/NCBI-NR/nr

#自定义输出格式

blastp -query 1PP2.L.fasta -db /data7T/mercu-b-backup/pub-data/NCBI-NR/nr -out 1PP2.L.fasta.out.0 -num_threads=8 -outfmt=’6 qseqid sseqid pident length slen evalue’

http://www.ncbi.nlm.nih.gov/blast/html/sub_matrix.html

短序列 psiblast -query AT4G13850.1 -db ../../../conservation/subject/AT4G13850.1 -out AT4G13850.1.out -num_iterations 3 -evalue 20000 -matrix PAM30 -comp_based_stats 0

psiblast -query $(basename $@) -db /data7T/mercu-b-backup/pu b-data/NCBI-NR/nr -out $@ -outfmt 0 -num_threads 8 -threshold 0.005 -eval ue 10 -max_target_seqs 1000 -num_iterations 10 -inclusion_ethresh 0.002

#file.fasta由clustralw等多序列比对软件生成

psiblast -in_msa file.fasta ../../../conservation/subject/AT4G13850.1 -out AT4G13850.1.out -num_iterations 3 -evalue 20000 -matrix PAM30 -comp_based_stats 0

#返回fasta格式的数据

blastdbcmd -entry NP_172083.1 -db /data7T/mercu-b-backup/pub-data/NCBI-NR/nr

#自定义输出格式

blastp -query 1PP2.L.fasta -db /data7T/mercu-b-backup/pub-data/NCBI-NR/nr -out 1PP2.L.fasta.out.0 -num_threads=8 -outfmt=’6 qseqid sseqid pident length slen evalue’

http://www.ncbi.nlm.nih.gov/blast/html/sub_matrix.html

3.meme  [http://meme.sdsc.edu/meme/doc/meme.html]

meme F1_MEF.macs_peaks.pv2000.bed.fa -dna -nmotifs 1000 -evt 0.001 -revcomp \
-mod anr -minsites 50 -maxsites 10000 \
-maxsize `cat F1_MEF.macs_peaks.pv2000.bed.fa | wc -c` \
-oc F1_MEF.macs_peaks.pv2000.meme

 

 

4.meme修改

mcast调用mhhmscan出现错误 sh: line 1: 11764 Aborted

FATAL: mhmmscan failed with exit status 134, command run was:

 

Revise:

run with another parameter –blocksize

 

$(cell).macs_peaks.pv2000.homer: $(cell).macs_peaks.pv2000.bed

findMotifsGenome.pl $< mm9 $@ -size given -p 6
pythonmail2.py $@ $@

$(cell).macs_peaks.pv2000.meme-fimo-uniprobe: $(cell).macs_peaks.pv2000.bed.fa
fimo –output-qthresh 0.1 –max-seq-length 2.5e8 \
–oc $@ –max-stored-scores 1000000000 \
~/soft/meme/db/motif_databases/uniprobe_mouse.meme $<
pythonmail2.py $@ $@
tail -n +2 $@/fimo.txt | cut -f 1 | sort | uniq -c | sed ‘s/^ *//’ \
| awk ‘BEGIN{OFS=”\t”;FS=” “;}{print 1}’ | sort -k1,1 \

$@/fimo.txt.motif.pattern.cnt
join -a1 -t ‘ ‘ -o 2.2,1.2 \
$@/fimo.txt.motif.pattern.cnt $(uniprobe_tf) \
$@/fimo.txt.motif.tf.cnt
$(cell).macs_peaks.pv2000.meme-mcast-uniprobe: $(cell).macs_peaks.pv2000.bed.fa
mkdir -p $@
mhmm –type star –keep-unused –verbosity 2 \
/home/chentong/soft/meme/db/motif_databases/uniprobe_mouse.meme \
$@/uniprobe_mouse.meme.mhmm
mhmmscan –fancy –allow-weak-motifs –p-thresh 0.0005 \
–max-gap 50 –e-thresh 10 –eg-cost 1 \
–blocksize \
`gawk ‘{if(x<length()) x = length()}END{print x}’ $<` \
–verbosity 5 –pseudo-weight 4.0 \
$@/uniprobe_mouse.meme.mhmm $< \
$@/mcast.txt
mhmm2html $@/mcast.txt >$@/mcast.html
pythonmail2.py $@ $@

 

5.

6.

7.

8.

9.

10.

11.

12.

13.

14.

15.

16.

17.

18.

19.

20.

CHENTONG
版权声明:本文为博主原创文章,转载请注明出处。
alipay.png WeChatPay.png

CHENTONG

CHENTONG
积微,月不胜日,时不胜月,岁不胜时。凡人好敖慢小事,大事至,然后兴之务之。如是,则常不胜夫敦比于小事者矣!何也?小事之至也数,其悬日也博,其为积也大。大事之至也希,其悬日也浅,其为积也小。故善日者王,善时者霸,补漏者危,大荒者亡!故,王者敬日,霸者敬时,仅存之国危而后戚之。亡国至亡而后知亡,至死而后知死,亡国之祸败,不可胜悔也。霸者之善著也,可以时托也。王者之功名,不可胜日志也。财物货宝以大为重,政教功名者反是,能积微者速成。诗曰:德如毛,民鲜能克举之。此之谓也。

生信宝典文章集锦

### 程序学习心得* [生物信息之程序学习](http://mp.weixin.qq.com/s?__biz=MzI5MTcwNjA4NQ==&mid=2247483927&idx=1&sn=23adf2b9d13400f2081f790e674e...… Continue reading

R统计绘图 - 柱状图

Published on August 12, 2017

R 学习 - 维恩图

Published on August 01, 2017