This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Align to a single gene -- best procedure

Hello,

I have a set of files derived from a WGS experiment. Using Qualimap, I found a mean coverage of about 30x, which was expected.

I now need to concentrate on a single gene, so I aligned these files against the reference for that gene instead that all the human genome. Now the coverage is about 1000x. Something is not right.

So, my question is: what is the best practice to align for a specific reference? Shall I align the fastq files against the reference of the gene of interest (as I have done now)? Or against the whole human genome (since the main experiment was on human WGS) and extract from there?

Thank you

alignment wgs best practice
• 1,684 views • link
written by marongiu.luigi marongiu.luigi

Something is not right.

Are you just looking for confirmation since you already know what you should do?

Or against the whole human genome (since the main experiment was on human WGS) and extract from there?

• 0 views • link
written by GenoMax GenoMax

Microfuge confirmed my fear: the aligner did not work properly and I should have used the second approach. So essentially, yes I was looking for a confirmation.

• 0 views • link
written by marongiu.luigi marongiu.luigi

1 answer

I would strongly be for aligning to whole genome and not to a specific gene (unless it was targeted sequencing and not WGS). In absence of a complete reference the aligner would force reads from other regions on to partial reference.

• 0 views • link
written by microfuge microfuge

Log in to answer this question.

More posts like this