serotyping from assemblies as input

[abrozzi]

Dear Friends,

I don't have reads, but assemblies (one or mere contigs) of S. pneumoniae.
How can I assign the serotype?
seroBa doesn't look like to support fasta input right?

Bests,
Alex

[hujingchu]

Hi Alex,
Did you find any solutions to serotyping directly from assemblies directly?

[abrozzi]

Dear Hu,
thank you for reaching out!
No, I didn't. I found a workaround downloading from NCBI Bioproject serotyping information (when present).
It would be nice to have for Pneumonia a tool like Kaptive https://github.com/katholt/Kaptive/wiki/Interpreting-the-results

I will keep you updated if I am aware of a serotyping directly from assemblies directly.

With Regards

Alessandro Brozzi,
GSK Vaccines (Italy)

[hujingchu]

Dear Alex,
Thanks you for the information. I tried by mapping reference sequence of different serotype directly to the assembly genomic sequence via minimap2, then choosing the the best hit (mapping quality, length and diversity rate), after checking couples of assemblies from refseq, it seems works
Best,
JC