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