This is a test version of Biostars. For the public version, visit https://www.biostars.org.
VariantCalling with single sanger read

Hi All

I want to do Variant Calling 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 Position of Homozygous and Heterozygous events.

I used sangerseqR package for processing and working with sanger read, now I have a single sequence(primary sequence, not sanger anymore) and a reference (that 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

snps blast indel variant-calling
• 3,569 views • link
updated by Ram Ram • written by Amirosein Amirosein

5 answers

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

• 1 views • link
written by DG DG

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.

• 1 views • link
written by i.sudbery i.sudbery

Absolutely second this.

• 1 views • link
written by DG DG

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/

• 1 views • link
written by Charles Warden Charles Warden

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.)

• 1 views • link
written by WouterDeCoster WouterDeCoster

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 ?

• 1 views • link
written by jawhar.saks • jawhar.saks

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.

• 1 views • link
written by jawhar.saks • jawhar.saks

Indigo outputs a BCF file with variant calls for Sanger sequencing: https://www.gear-genomics.com/indigo/

• 1 views • link
written by trausch trausch

Log in to answer this question.

More posts like this