Closed NikoLichi closed 6 months ago
Hi Wouter,
I gave it a try again with the feather files and have exactly the same errors. All the other files are produced but not the final compiled file.
How this can be solved?
This is the code I am using after feather:
NanoComp -t 32 --verbose -f pdf --feather $FILEIN -p NovoVsGeneC -o NovoVsGenC/qualCont/nanopack --names $PREFIX
Thanks and all the best, Niko
Hi Niko,
That is quite remarkable. Could you share the NanoStats.txt file?
Thanks, Wouter
Hi Wouter,
If you mean the .log
file after the Nanocomp run, it is enclosed below. Otherwise, please let me know which file you refer to.
All the best, Niko NovoVsGeneC_5don5TimPYFNanoComp_20240322_1426.log
Was there no NanoStats.txt file? That should also be generated in /NovoVsGenC/qualCont/nanopack/
I completely missed the file across all the other files, sorry. Here it is.
Aha now I see! I see the read identifiers on line 21 are like "141:329|2a491583-c20d-48e9-8ccc-49afe630be59". Is this a duplex run?
Oh... interesting No duplex run. This is cDNAseq from an RNA isolation protocol.
This is the output after trimming and finding directionality using Pychopper. They add those identifiers (e.g.,"141:329|") before the actual read identifier.
But... nanocomp is able to process all the metrics with those headers in separate files (HTML and PDF), but the consensus globall output (HTML) fails.
Yes, most of the time it doesn't care about that :
, but it is used when generating the HTML report of the NanoStats file. Let me see if I can come up with an easy fix.
As a workaround, could you see if running with --tsv_stats
fixes things?
Yes! It worked! Thanks!
I ran the command for two data sets already, and it works fine! I'll keep this command trick in mind when using pychopper sequences.
All the best, Niko
Thanks for the feedback!
Hi Wouter,
I used nanocomp to compare different runs on technical samples. It runs well before printing the whole report. Please see the errors below. This may be related to the issue #41 I will give a try to the create_feather script and see how it goes...
Thanks for any other help, Niko