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
•
written
by
Amirosein
1 answer
You don't really fix this in software, you sequence with primers in both directions so you can approach the indel from both sides.
So use internal primers. Sanger sequencing isn't really designed in a way to meaningfully allow reading over indels.
Log in to answer this question.
-
Can sangeranalyseR be changed to read SCF files?
written by Pedro A D P •
-
SNP filtering
written by rheab1230 •
-
Figure out genotype from VCF line
written by imbiling •
-
Documentation on VCI files? Or convert VCI to VCF?
written by QVINTVS_FABIVS_MAXIMVS
-
Creating A New File by Combining Two Separate Files
written by Tastulek •
-
What does ALT 0 mean in a VCF?
written by emyli •
-
VariantCalling with single sanger read
written by Amirosein
-
How to get minor allele frequency of a list of variants; for example NM_004401 (c.T122>G) from 1000…
written by Sudheer Gupta •
-
Python - finding INDELS and SNPS
written by will.duncan26 •
-
Is This An Error Of Ref Dbsnp?
written by Tonyzeng •
Try polyphred. I used to work on a similar tool ~14 years ago in BGI, but it never got it good enough.
as i know, Polyphred could not handle Indels but only handle Substitutions.
As I remember, it did, but I could be wrong – haven't used it for 10+ years.