Answers.Body: <blockquote>
  <p>maybe you can open the demovir.R file with Rstudio and execute line by line to check where the failure is. - @zorbax</p>
  </blockquote>
  <p>Thanks for your suggestion. I eventually got a correct final output already.</p>
  <p>The key issue was Prodigal would add extra strings after my original contig IDs, so it would cause R not recongnising my viral contigs. Therefore, I added a flag <code>-trunclabels</code> in Usearch (after e-value flag) to truncate my contig names. Finally, I re-ran the demovir.sh in ubuntu VM and it has worked well.</p>
  <p>Ernie</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2021-09-25T11:35:02.923Z
Answers.LastActivityDate: 2021-09-29T04:27:48.183Z
Answers.LastEditDate: 2021-09-29T04:27:48.183Z
Answers.UserDisplayName: Ernie Hsieh
Body: <p>I am new to bioinformatics and am working on viral metagenomic (virome) datasets and so would like to use <a href="https://github.com/feargalr/Demovir" rel="nofollow noreferrer">Demovir*</a> for doing taxonomic annotations on my viral contigs.</p>
  <p>First, I installed all prerequisites in my macOS but I found that Usearch 32-bit version could not run on macOS Big Sur, so I installed all and re-ran <code>demovir.sh my.viral.contig.fasta threads</code> in ubuntu VM. I could therefore successfully obtained trembl_ublast.viral.u.contigID.txt as an input for R, but I would always get an empty output after performing the demovir.R command of demovir.sh in ubuntu VM.</p>
  <p>Thus, I took these back to my macOS, used conda to create a virtual environment for R-3.5.1 and installed R using <code>conda install -c r r=3.5.1</code>, then ran <code>Rscript demovir.R trembl_ublast.viral.u.contigID.txt</code>. Again, I got an empty output with only showing the headers of five columns.</p>
  <p>I also tried to modify the demovir.sh script by removing all previous parts I've completed and retaining the last command <code>Rscript $DIR/demovir.R</code>then ran <code>modified.demvir.sh my.viral.contigs.fasta threads</code>again, but still got the same result.</p>
  <p>Hopefully anyone can do me a favour to check and resolve this issue.</p>
  <p>Many thanks.</p>
  <p>Ernie</p>
  <hr />
  <p>Quoting the GitHub repository:</p>
  <blockquote>
  <p>&quot;Democratic taxonomic classification of viral contigs to Order and
  Family level</p>
  <p>When performing metagenomic sequencing of Viral-Like Particle (VLP),
  the majority of returned sequences often bare little to no homology to
  reference sequences - Viral Dark Matter. Frequently it may be useful
  to know which viral taxonomic group these novel viruses are likely to
  belong to as this will give information about nucleic acid type, size
  and behaviour.&quot;</p>
  </blockquote>
  
Tags: <sequence-annotation><covid-19><metagenome><taxonomy><virus>
Title: Demovir produces an empty output
ThreadText: # Demovir produces an empty output
  
  (Asked by: username_0 on 2021-09-22)
  
  I am new to bioinformatics and am working on viral metagenomic (virome) datasets and so would like to use [Demovir\*](https://github.com/feargalr/Demovir) for doing taxonomic annotations on my viral contigs.
  
  First, I installed all prerequisites in my macOS but I found that Usearch 32-bit version could not run on macOS Big Sur, so I installed all and re-ran `demovir.sh my.viral.contig.fasta threads` in ubuntu VM. I could therefore successfully obtained trembl\_ublast.viral.u.contigID.txt as an input for R, but I would always get an empty output after performing the demovir.R command of demovir.sh in ubuntu VM.
  
  Thus, I took these back to my macOS, used conda to create a virtual environment for R-3.5.1 and installed R using `conda install -c r r=3.5.1`, then ran `Rscript demovir.R trembl_ublast.viral.u.contigID.txt`. Again, I got an empty output with only showing the headers of five columns.
  
  I also tried to modify the demovir.sh script by removing all previous parts I've completed and retaining the last command `Rscript $DIR/demovir.R`then ran `modified.demvir.sh my.viral.contigs.fasta threads`again, but still got the same result.
  
  Hopefully anyone can do me a favour to check and resolve this issue.
  
  Many thanks.
  
  username_0
  
  ---
  
  Quoting the GitHub repository:
  
  > "Democratic taxonomic classification of viral contigs to Order and
  > Family level
  >
  > When performing metagenomic sequencing of Viral-Like Particle (VLP),
  > the majority of returned sequences often bare little to no homology to
  > reference sequences - Viral Dark Matter. Frequently it may be useful
  > to know which viral taxonomic group these novel viruses are likely to
  > belong to as this will give information about nucleic acid type, size
  > and behaviour."
  
  **Tags:** sequence-annotation, covid-19, metagenome, taxonomy, virus
  
  
  ---
  
  ## 1 answers
  
  ### Answer by username_0 on 2021-09-25. Score: 2 | Accepted answer
  
  > maybe you can open the demovir.R file with Rstudio and execute line by line to check where the failure is. - @zorbax
  
  Thanks for your suggestion. I eventually got a correct final output already.
  
  The key issue was Prodigal would add extra strings after my original contig IDs, so it would cause R not recongnising my viral contigs. Therefore, I added a flag `-trunclabels` in Usearch (after e-value flag) to truncate my contig names. Finally, I re-ran the demovir.sh in ubuntu VM and it has worked well.
  
  username_0
  
  
  ---
  
  