VariantCalling with single sanger read
5
0
Entering edit mode
5.1 years ago
Amirosein ▴ 70

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

VariantCalling snps indel blast • 1.4k views
ADD COMMENT
1
Entering edit mode
5.1 years ago
DG 7.2k

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

ADD COMMENT
1
Entering edit mode

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 REPLY
1
Entering edit mode

Absolutely second this.

ADD REPLY
1
Entering edit mode
5.1 years ago

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 COMMENT
1
Entering edit mode
5.1 years ago

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 COMMENT
0
Entering edit mode
2.7 years ago
jawhar.saks ▴ 10

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 COMMENT
0
Entering edit mode

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 REPLY
0
Entering edit mode
2.7 years ago
trausch ★ 1.6k

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

ADD COMMENT

Login before adding your answer.

Traffic: 2022 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