I would like to make a fake quality file for a Velvet de novo assembly. What is the average quality of a velvet assembly?
I've been told that Velvet assemblies are basically infallible. If they are, then for an assembly I might say that there is less than one error in the entire genome,
Where G is the size of the genome in nucleotides. If the size of my bacterial genome is 2 MB, then the quality of the assembly would be greater than -10log(1/(210^6)), or 63.
So, does anyone have a good idea what a Velvet quality file would look like? Do you agree with my "63" assessment?