10 SeparateChromosomes2划分连锁群不合理怎么解决

我分析的群体是果树的F1,具有25条染色体,因此我想要分成25个连锁群。利用zcat p.call.gz | java -cp /mnt/e/QTL/20240813HLVCF/my_genetic_map/01.PrepareData/Lep-MAP/bin SeparateChromosomes2 sizeLimit=30 numThreads=10 data=-  lodLimit=$lod usePhysical=1  0.000000000001 > LOD.$lod.map 2> LOD.$lod.log   

done将2000多个(总计4000多个bin)bin能划分为25个连锁群,且进一步查看不同染色体之间划分是否合理发现,同一条染色体上的bin划分在各种不同的连锁群。这个问题该怎么解决呢。attachments-2024-08-dW7Ph3py66c1bea686040.png

请先 登录 后评论

2 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

你这个不是染色体,是contig,这些短的contig组装不可靠可以删除一下;

请先 登录 后评论
做个梦给你

我是在按照参考基因组的命名方式来命名为contig的,所以这应该就是我的染色体attachments-2024-08-8NSnIaBb66c2e7c94eda7.png

请先 登录 后评论