Question: Why do I have displacement in results?
0
gravatar for agata88
3.1 years ago by
agata88770
Poland
agata88770 wrote:

Hi,

I am analysing yeast strains NGS results containing amplicons for fragments of genes. I've mapped reads to whole genome S288C and EC1118. I compared the results. One amplicon were mapped to fragment of gene A in S288C reference. This amplicon matches in reference coordinates for gene A. The lenght of amplicon is 600bp.

I've done the same thing to EC1118 reference and I've discovered that this amplicon were mapped to about 75% of the same fragment of gene A. The mapping results showed that reads were shifted about 150bp. (and cover intron sequence).

Gene A has different coordinates in references and I took this into analysis.

Why do I have displacement in results?

Could it be caused by some kind of rearrangements?

I would really appreciate for help. Best, Agata

dna mapping ngs • 724 views
ADD COMMENTlink written 3.1 years ago by agata88770

Do you have paired end data? Does one end of the amplicon end in intron sequence?

This is just speculation but it is not improbable to imagine that your lab strain may be subtly different from published references because of its age/handling (especially if it has been some years since it was originally acquired).

ADD REPLYlink written 3.1 years ago by genomax65k

Yes, I have PE reads and one end of amplicon ends in intron sequence and second ends in gene.

In that case I would expect the same results for S288C reference, I think. Do you think that this displacement could be caused by insertion? But the amplicon length is 600 in both cases ...so nope. Any other idea?

ADD REPLYlink written 3.1 years ago by agata88770

Ok, But insertion/deletion can be in nearest area, like in other gene! In that case my whole amplicon will move. Because I see that all reads were moved in left I think this may be caused more possibly by deletion. Correct me if I am wrong, please. Agata

ADD REPLYlink modified 3.1 years ago • written 3.1 years ago by agata88770
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: 1033 users visited in the last hour