Question: Align to a single gene -- best procedure
2
gravatar for marongiu.luigi
11 days ago by
Germany, Mannheim, UMM
marongiu.luigi420 wrote:

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

best practice alignment wgs • 58 views
ADD COMMENTlink written 11 days ago by marongiu.luigi420

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?

ADD REPLYlink modified 11 days ago • written 11 days ago by genomax74k

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.

ADD REPLYlink written 11 days ago by marongiu.luigi420
2
gravatar for microfuge
11 days ago by
microfuge1.4k
microfuge1.4k wrote:

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.

ADD COMMENTlink written 11 days ago by microfuge1.4k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1498 users visited in the last hour