Open AndAvia opened 2 months ago
Hi Chuang Ma
Your result sounds quite different. You could run a tool like Vibrant to get closer to the set of 4422 vOTUs used in the analysis (see article "Furthermore, we also parsed vOTUs using VIBRANT (v.1.2.1, default settings) and found a 4240/4422 (96%) viral prediction agreement.").
hi, I download the vOTU.fna (There are 5251 vOTUs in there.)data in the text and used the MGV pipeline(https://github.com/snayfach/MGV/tree/master/viral_detection_pipeline), like this
During the run no error was reported, all steps were 5251, after running the last step python viral_classify.py, finally only got 2897 viruses, in the literature the result was 4422, where is the problem, I see viral_classify.py used to a rules file classification_rules.tsv, is that file different from what you guys were using. Thanks for answering the questions!