8 个回答
求助mRNA ID 批量+1问题
老师，利用这个脚本修改，awk -F"[\t=]" 'BEGIN{OFS="\t"}{if($3=="mRNA"){print $1,$2,"gene",$4,$5,$6,$7,$8,$9"="$10 ; print $1,$2,$3,$4,$5,$6,$7,$8,$9"="$10"Parent="$10;}else{print $0}}' PO.gene.gff >PO.gene.gff1 可否提示一下这个脚本修改那一部分可以实现mRNAID后面批量加一个.1
回答于 2020-11-13 15:00
MCSCANX共线性图代码报错
我的是family.ctl里面的染色体和基因组染色体的ID号不一样，报错Exception in thread "main" java.lang.NullPointerException at family_circle_plotter.paint(family_circle_plotter.java:197) at family_circle_plotter.main(family_circle_plotter.java:334) ，修改一致后，完美运行，希望能帮到同样问题的同学，也感...
回答于 2020-08-30 00:18
利用gff文件提取基因结构，结果文件是空的
use Getopt::Long; my %opts; use Data::Dumper; GetOptions( \%opts, "in1=s", "in2=s", "out=s", "h" ); if ( !defined( $opts{in1} ) || !defined( $opts{in2} ) || !defined( $opts{out} ) || defined( $opts{h} ) ) { &USAGE; } open( IN1, "$opts{in1}" ) || die "open $opts{in1} failed\n"; ope...
回答于 2020-08-20 02:09