Benjamin Buchfink
Benjamin Buchfink
`-b 0.001` is not a good setting for this parameter, why are you using it? How much memory does your machine have?
Maybe the target has no taxonomic assignment or something, but I don't know how Megan works. You can ask your question here: https://megan.cs.uni-tuebingen.de/
Sorry that this completely got lost! Since v2.1.0 you can now use `--header simple` to just print the field keywords as headers.
> Is there any update on this, or suggestions for a workaround? You can use the diamond database format which supports taxonomic information.
> Hi! Also I find that `samtools fastq` also doesn't work, with the error > > ``` > [E::sam_hdr_create] Invalid header line: must start with @HD/@SQ/@RG/@PG/@CO > Failed to read...
That should not happen, please let me know your diamond version and command line to reproduce the problem.
I have tried but can't reproduce the problem, would it be possible to send me your fasta input files?
Yes that would probably make sense changing that, I'll look into it.
There's the `--taxonlist` and `--taxon-exclude` options (see wiki).
It will not search the respective database sequences, but it will not remove anything from your query file, if that's what you mean.