Quake Coverage Cutoff output

Question

[Seigfried]

Hello I am trying to run the Coverage cutoff module of Quake Error Corrector. i get output files called cutoff.txt and kmers.hist

My Cutoff.txt output is as follow:

1 zp.copy: 0.8268095 p.e: -50570.72 shape.e: 253.7109 u.v: 73.32117 var.v: 507.014 -18066.49 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

Can anyone help me understand this file? I could not find anything on the manual and my result seems ambiguous. Is my coverage cutoff value 1?

[st.ph.n]

At the bottom of your cov_cutoff_"kmer_size" file you you should see 'Cutoff: #'. That number should be your cutoff value given the value of the kmer you used.

Answer 1

[Seigfried]

the current Cov_model lists the CUtoff at the end on the command line.

In the file produced the cutoff is the "1" the first number in the file.