SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
variant calling using SeqIO ksw9 Bioinformatics 2 11-01-2014 03:35 PM
Reliable variant calling Mudal Pacific Biosciences 3 05-12-2014 12:12 AM
Doubt about .abi files and .fasta files fc35802 Bioinformatics 2 02-28-2012 10:05 AM
variant calling kjaja Bioinformatics 1 11-04-2011 07:16 AM

Reply
 
Thread Tools
Old 07-21-2017, 12:03 AM   #1
dimo
Member
 
Location: Australia

Join Date: Mar 2012
Posts: 10
Default Variant calling from fasta files.

Greetings and salutations.

I am trying to analyze some old data which consists of several hundred fasta sequences from individual tumour samples which I believe came from some sort of Sanger sequencing. Is there tool that will align these to a genome (BLAST) and call nucleotide and amino acid variants (similar to a NGS analysis)?

Thanks in advance.
dimo is offline   Reply With Quote
Old 07-21-2017, 12:45 AM   #2
vivek_
PhD Student
 
Location: Denmark

Join Date: Jul 2012
Posts: 164
Default

You could likely use BWA-Mem to do the alignment and call variants using the GATK's UnifiedGenotyper module.
vivek_ is offline   Reply With Quote
Old 07-21-2017, 01:19 AM   #3
dimo
Member
 
Location: Australia

Join Date: Mar 2012
Posts: 10
Default

Quote:
Originally Posted by vivek_ View Post
You could likely use BWA-Mem to do the alignment and call variants using the GATK's UnifiedGenotyper module.
Interesting, any idea how GATK would handle a lack of quality scores? I would also basically have only 1 read per sample as well so no idea if it would even flag.
dimo is offline   Reply With Quote
Old 07-21-2017, 02:50 AM   #4
vivek_
PhD Student
 
Location: Denmark

Join Date: Jul 2012
Posts: 164
Default

Quote:
Originally Posted by dimo View Post
Interesting, any idea how GATK would handle a lack of quality scores? I would also basically have only 1 read per sample as well so no idea if it would even flag.
If it complains about quality scores, you can just generate a file with some dummy scores. I'm not sure how it would handle one read per sample though. You can also try samtools for variant calling and modify the min. reads flags.
vivek_ is offline   Reply With Quote
Old 07-21-2017, 06:55 AM   #5
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

You can use BBMap's reformat.sh to convert fasta to fastq, and assign whatever fake quality score you want (the default is Q30, though you can adjust that with the flag "qfake=X"). I like 30 because it ends up with quality scores that are all "?", which is appropriate.

Also, I feel like 30 is a good insurance against bad trimming - I would never, ever recommend trimming to Q30 or above (meaning bases with a claimed quality below 30 are discarded), since it is virtually always a bad idea if you have accurate quality scores, but some pipelines do that automatically. So, a default of Q30 seems reasonable.

Last edited by Brian Bushnell; 07-24-2017 at 08:36 AM.
Brian Bushnell is offline   Reply With Quote
Old 07-24-2017, 06:57 AM   #6
colindaven
Senior Member
 
Location: Germany

Join Date: Oct 2008
Posts: 415
Default

I've done this once using freebayes. It will happily call variants in VCF format from a BAM file with a coverage of 1. Took me a while to find this solution.
colindaven is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 08:16 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO