Hi!
I am trying to automate variant calling on human short exons sequenced with Sanger method. I have the ab1 files, and I want to automate the process for clinical geneticists. The program should be free, I understand that is a huge limitation. Also, it needs to register INDELS in addition to SNPs.
I tried working with tracy, and it works fine but fails silently, and fails to index the reference (I have opened github issues for that). TraceTuner doesn't seem to do anything when I use the -het flag, it just outputs the single consensus fasta (no heterozygous calls). They currently work on Sequencher trial version, but this version doesn't do variant calling. Baseling also doesn't mention automated heterozygous calls. Although there are mentions of writing in ambiguous bases which could be a code for a heterozygous, this becomes a problem for INDELs.
Is there a free resource that somebody uses successfully for human diploid data variant calling? Thanks all!
1 answer
The problem with the indexing in tracy was the leading whitespace, right? If there is any other issue please let us know.
tracy decompose -f align -o out -a homo_sapiens -g Homo_sapiens.GRCh38.dna.primary_assembly.fa.gz trace.ab1
Tracy is still work-in-progress so any feedback is very welcome. The above command generates alignments, decomposes InDels, annotates variants with rs identifiers and generates a BCF file for downstream processing. The reference genome index files can be downloaded here and please let us know if your species of interest is missing.
Log in to answer this question.
-
dbSNP VCF file compatible with GRch38.p14
written by Kelly •
-
What is the set of all possible GT calls for a VCF or gVCF for human data.
written by ariel
-
Samtools mpileup - I'm getting different number of base calls compared to number of quality scores
written by zuhaibzulfiqarahmed •
-
Index genome not working in Tracy
written by Joe •
-
VG Toolkit for human WGS analysis
written by venkat •
-
Extracting FILTER field from gnomAD vcf file
written by Damianos P. Melidis •
-
Indel realignment: GATK vs. BCFtools
written by martytriggers •
-
Tool: Hera: A new tool for RNA-Seq analysis
written by sonpham
-
News: Workshop: GATK Best Practices for Variant Discovery, Edinburgh, 17-19 July 2017
written by Bert Overduin
-
Samtools Mpileup And Bcftools Parameter
written by Ashutosh Pandey
Maybe Phred/Phrap can help here, but I am not sure. See this paper.
When I worked on Sanger data I used Phred/Phrap/Polyphred, and it was probably the best tool around. If you are using the results in clinical setting you might investigate the issues brought forward by genomax with your clinical colleagues.
While I understand a desire to use a free solution have you checked applicable local administrative policy for this. You may need (or be required) to use a validated/approved (and thus perhaps not free) application, if the data is going to be used for patient diagnostics.
This is a very valid concern. If all fails, it will be a nice way to screen or double-check.