Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Downloading Data failure and QV calculation #42

Open
Yoon90 opened this issue Jun 25, 2024 · 1 comment
Open

Downloading Data failure and QV calculation #42

Yoon90 opened this issue Jun 25, 2024 · 1 comment

Comments

@Yoon90
Copy link

Yoon90 commented Jun 25, 2024

I've tried to download your corrected data via:
wget http://complex.zesoi.fer.hr/data/downloads/HG002.experimentalUL.corrected.fasta.gz
But it doesn't seem like work and i could not download the data. Could you check the link again?

And one more thing, How did you calculated the QV of corrected reads? - as far as i know the output format of corrected reads is .fasta which does not have QV. I'm sorry if i'm wrong about this.

@dehui333
Copy link
Collaborator

dehui333 commented Aug 7, 2024

Hi,

Please try again with the following:
wget -O HG002.experimentalUL.corrected.fasta.gz https://zenodo.org/records/13252378/files/HG002.experimentalUL.corrected.fasta.gz?download=1


Regarding QV,

The corrected reads themselves do not carry information about their quality.

We obtained some form of measurement of the quality of reads via two ways:

  1. Aligning the reads to a high quality assembly of the same sample (if one exists) and measure the agreement between these reads and the aligned-to segment (e.g. using bamConcordance (https://github.com/PacificBiosciences/hg002-ccs/blob/master/concordance/bamConcordance) to get some form of Phred-scaled score, or just calculate accuracy from the CIGAR string).
  2. Using yak (https://github.com/lh3/yak) to obtain a QV using illumina short reads from the same sample.

Best,
Dehui

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants