This is a test version of Biostars. For the public version, visit https://www.biostars.org.
TP, FP, TN and FN in error correction reads

Hello,

I want to know how to calculate the TP, FP, TN and FN in the case of error correction?

Thanks

error-correction assembly reads

1 answer

The same way you would calculate them in any case, with a dataset that has known correct results that's then run through the procedure.

I don't have these results. I have just the Reference genome. Do I use the reference genome?

thanks

Yes, you would have to use a reference genome to generate fake data with known properties what would then be error corrected. Ideally you just sequence the reads, but that assumes that the reference is correct (maybe it's a good assumption, maybe not).

Or the best metric is the mapping of corrected reads to the reference genome. Am I right?

thanks

No, then you're adding another layer of processing and confounding things by how well the aligner is able to deal with the reads and the genome.

I do not have a choice. Otherwise how I validate my corrected reads??

I have the reference genome. So I can map my corrected reads to the reference genome and I see the results of mapping.

Log in to answer this question.