组装RNA-seq的数据后，不能确定是由于组装产生的还是因为本身发生啦同源重组才导致序列会出现相同。有什么办法可以验证啊？？？

terminate called after throwing an instance of 'mapped_file::ErrorMMap'   what():  [...]

我用以前写的heatmap.2写的热力学图 发现出了一个问题 在基因的个数在100到3000左右的范围内，图像会出现很多空白的横纹   这个是空白横纹如何去掉的 我想要得到一个比较完整的清晰的基因表达谱的热力学聚类图

各位帮帮忙，看看这个调用应该怎么写？这样写第二句，为什么会报错啊？谢谢大家了！！！ my $aln='$query.fa'; my $clustalw ="D:\Program Files\ClustalW2\clustalw2"; my $tempoutput='$query.aln'; my [...]

不知道有没有相关的review?

我想去除一组序列（大约1000条）中指定区段序列，例如 file1中包含原始待处理序列，格式和内容如下： >AT2G01210.1 MLASLIIFVALLCNVTVISGLN >AT2G01820.1 MSNSHLGTLCFIISLLGLANFSLS >AT2G01950.1 MTTSPIRVRIRTRIQISFIFLLTHL file2中包含指定区段信息，格式和内容如下： AT2G01210.1 [...]...

  The 5th Virtual Training Workshop on Bioinformatics Asian Bioinformatics [...]

    测了几个转录组，需要做GO和kegg；但是都不是模式生物，没有本物种的GO数据库，像这种情况怎么去做，希望得到大家的指导。     貌似Interpro可以，那么我需要手动先把我的序列按六种ORF框翻译成蛋白序列么？还是提供个map得到的Gene ID即可？？     例外kegg现在ftp已经收费了，是不是对分析大通量转录组数据很不利？看之前的帖子提到API不受影响，但是我理解API的实现方式也是需要上传到服务器端的数据库进行比对；这种海量数据不现实。不知道我猜测的对不对…… [...]...

各位大侠，小弟打算做一个Moran's I Autocorrelation Index检验，在EXCEl表中输入了成对物种间的遗传距离，但是，不知道该怎么导入R中，我按照普通的read.table()导入，结果，不能显示出数据时成对的矩阵，请各位大侠指点，小弟不胜感激！

RT 我听导师说 我们要比对的DNA序列要比对2 3天的时间 .... 如果有软件能不能给我发一个 我的邮箱是zhaozhiqiangg@126.com 刚涉及到这个 可能问题有点幼稚 请见谅

我通过表达谱芯片拿到一组差异表达的酿酒酵母的gene symbol 想找其在染色体上的start 与 end 例如：   start stop Symbol [...]

举例： a1=t(matrix(letters[1:10],nrow=2)) >a1      [,1] [,2] [1,] "a"  "b" [2,] [...]

利用haploview得到了位点间的r2值，我想做LD 衰减图，是把结果里边所有的r2值都拿来用呢，还是有什么过滤条件之类的？谢谢大家

在本地blast的时候，因为没有输入-m7的选项，结果出来是默认的格式，但是现在想把输出文件改成-m7的格式，不知道怎么实现？因为数据量比较大，不想再跑一次blast，谢谢！