找到约 15 条结果

问题 关于基因受选择计算

...体中是否基因受到选择,我用到的命令是 vcftools --vcf sample.vcf --weir-fst-pop population_1.txt --weir-fst-pop population_2.txt --out P_1_2 --fst-window-size 10000 --fst-window-step 1000 然后选取结果中top1%的区段和我基因区段做overlap,然后去看是否受到选...

问题 请问差异代谢物数量过多能否调高筛选条件?最大能提高多少?

利用VIP&gt;1;p<0.05筛选出的差异代谢物数量太多了,能不能吧VIP的阈值提高至1.3或者1.5呢?

问题 想请教一下,如果在stringtie定量后,得到的某基因有多转录本,那么它的定量结果,即FPKM值应该如何取舍?

...exon_number "2"; cov "14.629032"; 比如上面FAN_iscf00000040.1基因有3转录本,最后的FPKM值也有三,像这种情况该如何取舍?谢谢!

文章 必须收藏!比较基因组学分析必备神器-OrthoVenn2

...通用性更强,单次在线分析的上限物种数量提高了1倍达12。OrthoVenn2还提供了本地版本下载,可以处理更大数据量的分析,无物种数限制,分析环境封装于docker镜像中,人安装docker后直接下载OrthoVenn2镜像即可使用,还未安装d...

文章 第11期《基因家族分析实操培训班》邀您参加!

基因家族(Gene family),是来源于同一祖先,由一基因通过基因重复而产生两或更多的拷贝而构成的一组基因,它们在结构和功能上具有明显的相似性,编码相似的蛋白质产物。 基因家族分析,一老生常谈的分析内容...

问题 老师,我想问一下,我在预测顺式作用元件时,在提取出来的基因上游1500bp的文件中,有的基因是这样的

&gt;PGSC0003DMG400010704 nnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnntttttttttttaaaaaaaTGTTTAAG 前面出现的这些n正常吗?然后后面这种有些大写有些小写的也正常吗?

文章 MUMmer + SVMU : CNVs

...ta .  echo "lastz ref.fa[multiple]  $query.fa[multiple] --chain --format=general:name1,strand1,start1,end1,name2,strand2,start2,end2 --ambiguous=iupac  &gt; $query.lastz.txt 2&gt; $query.lastz.log" &gt;&gt;lastz.cmds.sh  echo "svmu $query.delta ref.fa $query.fa h $query.lastz.txt $query &amp;&gt...

问题 脚本如果直接从ensamble下载cds和pep还需要去提取吗,是不是可以直接转换格式

for species in Anser_brachyrhynchus Equus_caballus Homo_sapiens Mus_musculus Oryctolagus_cuniculus Ovis_aries_rambouillet Sus_scrofa.Sscrofa;do      genome=$datadir/$species.genome.fa.gz   gff=$datadir/$species.gff.gz   #保留蛋白编码基因   agat_sp_filter_feature_by_attribute_value...

文章 MySQL报错:ERROR 2002 (HY000): Can't connect to local MySQL server through socket '/tmp/mysql.sock'

...ERROR 2002 (HY000): Can't connect to local MySQL server through socket '/tmp/mysql.sock' (2) 出现问题原因: 有可能是 my.cnf 配置文件中设置了 [mysqld] 的参数 socket ,而没有设置[client]的参数socket mysql.sock 文件有什么用: mysql 支持 socket 和 TCP/IP 连...

文章 深入理解Linux shell中2>&1的含义 标准错误输出与输入

1 和 2 在 Linux 中代表什么 在Linux系统中0 1 2是一文件描述符 从上表看的出来,我们平时使用的 echo "hello" &gt; t.log 其实也可以写成 echo "hello" 1&gt; t.log 关于2&gt;&amp;1的含义 关于输入/输出重定向本文就不细说了,不懂的...

问题 fst_pi_select_sweep.r软件不可用

[root@eff2c0ebb6b8  08:44:02 /work/3.map/vcf]# fst_pi_select_sweep.r --fst F-G.fst  --pi1 FS_F.pi --pi2 FS_G.pi --zscore --log2 -A FS_F -B FS_G -c 0.05 -n fst-pi.F-vs-G  -f pdf bash: /work/3.map/scripts/fst_pi_select_sweep.r: /biosoft/miniconda/bin/Rscript: bad interpreter: No such file or direc...

文章 AI助手帮我解决如何修改kvm虚拟机root密码

...问题为例,看下AI的解决效果。 对话如下 :https://www.doubao.com/thread/w986bf96b90261a31 这里给出了3种方法,我们选择使用第一种方法: 其中用到了virt-edit命令,我们顺便问下他的下载和使用: 我试了一下发现密码没有修改成...

文章 linux系统中ls命令的常用用

...   每行只列出一文件 -h              用"K","M","G"来显示文件和目录的大小。 这就是比较常用的ls命令。

问题 微生物16S分析时候,两端测序数据拼接

...数据双端合并的时候,执行这项代码:  qiime vsearch merge-pairs \     --i-demultiplexed-seqs primer-trimmed-demux.qza \     --p-threads  7 \     --o-merged-sequences demux-joined.qza后,查看合并结果,发现合并前后数据少了特别多。例如:合并前正...

问题 KEGGGO富集分析报错问题

P和Q值设置为1,依旧提示报错:Failed to create bus connection: No such file or directory --&gt; No gene can be mapped.... --&gt; Expected input gene ID: LOC105382119,LOC105392385,LOC105383189,LOC105384039,LOC105384627,LOC105384151 --&gt; return NULL... Error in barplot.default(kegg, ...