Question: VariantCalling with single sanger read
0
Entering edit mode

Hi All

I want to do VariantCalling but with a single sanger read in a single specific region (so i can't use routine VariantCalling pipelines). for example i have a sanger read from a region of a gene and i know that gene, now i want to tell Genomic Poisition of Homozygous and Heterozygous events.

i used sangerseqR package for proccessing and working with sanger read, now i have a single sequence(primary sequence, not sanger anymore) and a reference(taht gene), is there any function or package which can compare them and return genomic positions of Homozygous events in that query sequence ?

thanks a lot

ADD COMMENTlink 3.4 years ago Amirosein • 70 • updated 11 months ago trausch ♦ 1.3k
1
Entering edit mode

You probably want to use the MakeBaseCalls function of sangerseqR. That will give you heterozygous events, all other positions are homozygous.

ADD COMMENTlink 3.4 years ago Dan Gaston 7.1k
Entering edit mode
1

Make sure to look at the actual traces at your variant position(s). When it comes to sanger seq, no algorithm i've come accross beats eyes on the trace for working out if the variant is real or not.

ADD REPLYlink 3.4 years ago
i.sudbery
4.7k
Entering edit mode
1

Absolutely second this.

ADD REPLYlink 3.4 years ago
Dan Gaston
7.1k
1
Entering edit mode

You can use the Staden package, but there might be a bit of a learning curve with getting used to the interface:

http://staden.sourceforge.net/manual/mutations_unix_toc.html#SEC1

http://staden.sourceforge.net/

ADD COMMENTlink 3.4 years ago Charles Warden 6.8k
1
Entering edit mode

Developed in my department and great piece of software: http://www.molgen.ua.ac.be/bioinfo/novosnp/

(I wasn't involved in the development, I'm just a user.)

ADD COMMENTlink 3.4 years ago WouterDeCoster 39k
0
Entering edit mode

Hello,

We developed an automatic software working on batches and allowing for optionnal conclusion files based on theorical mutation position following papers rules.

IAGE Tools Analytics : A solution to detect SNP and INDELS from Sanger in CRISPR analysis workflow.

Video of the software.

Is that what you are looking for ?

ADD COMMENTlink 11 months ago jawhar.saks • 10
Entering edit mode
0

By the way, our software outputs for each sanger file (.ab1) a folder with in it :

  • Plus.fasta
  • Minus.fasta
  • Plus_Wt.aln
  • Minus_WT.aln
  • Chromatogram

And the conclusion files for the batch if a theorical position of mutation is given.

ADD REPLYlink 11 months ago
jawhar.saks
• 10
0
Entering edit mode

Indigo outputs a BCF file with variant calls for Sanger sequencing: https://gear.embl.de/indigo/

ADD COMMENTlink 11 months ago trausch ♦ 1.3k

Login before adding your answer.

Powered by the version 1.8