This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Centrifuge Classifying Reads, but the Report is Empty

Hi all, I've been trying to run Centrifuge for metagenomic analysis. It appears to be classifying reads, but when it comes to the output report, the report is empty. I have my code I'm using and the output files below. Please let me know what I am doing wrong/how to fix this? Thanks!

My input line:

centrifuge --threads 32 --met-stderr -x p+h+v -1 /sciclone/scr10/echolley/igem/myco_meta/fastqs/trimmed/M07CI1_val_1.fq.gz -2 /sciclone/scr10/echolley/igem/myco_meta/fastqs/trimmed/M07CI1_val_2.fq.gz -S test1 --report-file test1_report

Output:

52 [bora] head test1
readID  seqID   taxID   score   2ndBestScore    hitLength   queryLength numMatches
LH00260:52:2255MVLT4:5:1101:29736:1042  unclassified    0   0   0   0   70  1
LH00260:52:2255MVLT4:5:1101:30011:1042  unclassified    0   0   0   0   300 1
LH00260:52:2255MVLT4:5:1101:31127:1042  unclassified    0   0   0   0   246 1
LH00260:52:2255MVLT4:5:1101:32309:1042  unclassified    0   0   0   0   194 1
LH00260:52:2255MVLT4:5:1101:33231:1042  unclassified    0   0   0   0   300 1
LH00260:52:2255MVLT4:5:1101:33053:1042  unclassified    0   0   0   0   212 1
LH00260:52:2255MVLT4:5:1101:34688:1042  unclassified    0   0   0   0   254 1
LH00260:52:2255MVLT4:5:1101:35481:1042  unclassified    0   0   0   0   288 1
LH00260:52:2255MVLT4:5:1101:33345:1042  unclassified    0   0   0   0   218 1

Report:

53 [bora] head test1_report 
name    taxID   taxRank genomeSize  numReads    numUniqueReads  abundance
metagenomics centrifuge

0 answers

No answers yet.

Log in to answer this question.