detection of in/del in sanger sequencing
1
0
Entering edit mode
6.1 years ago
Amirosein ▴ 70

Hi I'm working on Sanger sequencing reads to detect the variants and annotate them, i have only one Sanger read from specific regions (exon in example) of the genome of the patient, using sangerseqR package in R to work with it. in some patients sanger data, one of the two copies of the region have some indel so bases shift and signal start to have two peaks. along with some mismatches how to fix the data and find indel variants?

example :

from Mom: A C G T A C G T
from Dad: A C G - A C G T

so the sanger will report:

ACGTACGT
ACGACGTT

and from the 4th base there will be two peaks and basecaller may report this instead:

ACGACGGT
ACGTACTT

any idea to solve this? is there any package or function to handle that? any software who report indels? thanks all

sanger indel variant detection • 2.1k views
ADD COMMENT
0
Entering edit mode

Try polyphred. I used to work on a similar tool ~14 years ago in BGI, but it never got it good enough.

ADD REPLY
0
Entering edit mode

as i know, Polyphred could not handle Indels but only handle Substitutions.

ADD REPLY
0
Entering edit mode

As I remember, it did, but I could be wrong – haven't used it for 10+ years.

ADD REPLY
0
Entering edit mode
6.1 years ago

You don't really fix this in software, you sequence with primers in both directions so you can approach the indel from both sides.

ADD COMMENT
0
Entering edit mode

in some examples, patients have some indel at the beginning and some indel at ending bases so having both directions won't help that much, i designed my own algorithm but i'm searching if there's another idea or tool o compare with.

thanks

ADD REPLY
0
Entering edit mode

and also enabling to fix and handle this issue will let us pay less sequencing both directions,

ADD REPLY
0
Entering edit mode

So use internal primers. Sanger sequencing isn't really designed in a way to meaningfully allow reading over indels.

ADD REPLY
0
Entering edit mode

ok thanks, i'll check again :)

ADD REPLY

Login before adding your answer.

Traffic: 1423 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6