This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Tool to calculate per-nucleotide confidence score in genome assembly using short-read data

What would be a good tool or collection of tools to calculate a confidence score for each nucleotide in a genome assembly using short-read data, ideally taking into account both the reads pile-up and each read's sequencing quality score, and able to handle SNPs and INDELs?

quality score genome assembly
• 1,015 views • link
updated by Corentin Corentin • written by Sebastien_Vigneau • Sebastien_Vigneau

Hi, I do not really know what you mean by "confidence" score. Would it be a measure of the probability of misassembly at a specific position? Assemblies are often assessed as a whole and not in a per-base manner.

It may not be exactly what you want, but Pilon has a --vcf parameter to produce a .vcf file listing detailed information about base and indel evidence at every base position in the genome. See: https://github.com/broadinstitute/pilon/wiki/Output-File-Descriptions#vcf

• 0 views • link
written by Corentin Corentin

0 answers

No answers yet.

Log in to answer this question.

More posts like this