Centrifuge Classifying Reads, but the Report is Empty
0
0
Entering edit mode
15 months ago
echolley ▴ 20

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 • 448 views
ADD COMMENT

Login before adding your answer.

Traffic: 3280 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6