Question: genepredtogtf command fails
0
gravatar for toptant
3.3 years ago by
toptant0
toptant0 wrote:

Hi, I am trying to generate a gtf file using the txt gene annotation file of a virus genome and getting this message. This command line worked perfectly with hg19 or hg38 txt files...So I really don't know what I am doing wrong here.

$ cut -f2-11 XGenome.txt|genePredToGtf file stdin XGenome.gtf

"Expecting at least 10 words line 2 of stdin got 1"

txt to gtf conversion • 1.4k views
ADD COMMENTlink modified 2.7 years ago by Tye Kahn10 • written 3.3 years ago by toptant0

What are the first few lines of XGenome.txt?

ADD REPLYlink written 3.3 years ago by genecats.ucsc570
1
gravatar for Tye Kahn
2.7 years ago by
Tye Kahn10
Tye Kahn10 wrote:

I had the same error, and from what I can see, the error comes from the input file. That's because your initial file does not have at least 12 words/columns in all lines, and therefore, after using cut, the resulting words used as input to genePredToGtf are less than 10 (in your case, that happens in the second line).

Therefore, I would recommend you to check that line in your file.

ADD COMMENTlink written 2.7 years ago by Tye Kahn10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1499 users visited in the last hour
_