如果有基因对信息的话,最好还是绘制出来。 关于scaffold序列长度信息,通过代码“ ## 获得染色体长度 seqkit fx2tab --length --name ../01.data_prepare/Arabidopsis_thaliana.TAIR10.dna.toplevel.fa.gz|awk '{print $1,$NF}' >genome.len ” 是可以得到的,如果scaffold序列是在Arabidopsis_thaliana.TAIR10.dna...
回答于 2025-06-05 16:53
这个基因没有在染色体上的位置信息,只有在scaffold上有位置信息。如果想展示这个基因的话,只能将该scaffold也跟染色体一样去准备长度信息。在结果图中,将scaffold及在scaffold上的基因信息绘制出来。
回答于 2025-06-03 16:10
正确的测序应该是有比对结果的。先根据流程看上一步生成的结果有没有问题,运行脚本时看有没有报错日志。看看map结果中每个样本的比对率(*summary文件)是否正常,正常结果中有90%以上的reads会比对到基因上
回答于 2025-06-03 11:43
具体绘图可以参考操作笔记:https://www.omicsclass.com/article/1269 如果报错的话可以把报错发出来看一下
回答于 2025-06-03 09:15