Question: Align to a single gene -- best procedure
2
gravatar for marongiu.luigi
5 months ago by
Germany, Mannheim, UMM
marongiu.luigi490 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 • 121 views
ADD COMMENTlink written 5 months ago by marongiu.luigi490

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 5 months ago • written 5 months ago by genomax80k

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 5 months ago by marongiu.luigi490
2
gravatar for microfuge
5 months ago by
microfuge1.5k
microfuge1.5k 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 5 months ago by microfuge1.5k
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: 1998 users visited in the last hour