WGLab / RepeatHMM

a hidden Markov model to infer simple repeats from genome sequences
Other
34 stars 14 forks source link

Suggested settings for CCS reads #47

Open davidlougheed opened 2 years ago

davidlougheed commented 2 years ago

PacBio CCS reads have a very different error profile to standard PacBio CLRs. Are there suggested HMM settings for genotyping CCS reads, in a context where I do not have access to subreads (even though as suggested in the 2017 paper they'd be better-performing)?

liuqianhn commented 2 years ago

@davidlougheed I do not have specific parameter for PacBio CCS reads. But I assume that CCS reads have similar error profile as illumina reads, and you can try to use --SplitAndReAlign 2 --SeqTech "Illumina" to see whether you have different results.