SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Srst2 virulence genes detection issue aoi Bioinformatics 0 04-30-2015 01:54 AM
Resequencing & variant calling issue Jluis Bioinformatics 5 02-05-2013 02:12 AM
Variant Detection Software Review twaddlac Bioinformatics 0 03-23-2012 11:58 AM
Intra-host variant detection? kga1978 Bioinformatics 0 11-08-2011 06:39 PM
has anyone tried variant detection with RNA-seq using GATK? mard Bioinformatics 0 07-28-2011 12:13 AM

Reply
 
Thread Tools
Old 08-09-2015, 10:31 PM   #1
anyone1985
Member
 
Location: shanghai, chia

Join Date: Mar 2009
Posts: 67
Default issue about microbial variant detection

I use gatk, freebayes and samtools to calling the variant from microbial sequence. According to the common sense, I think I should use the haploid model like gatk or freebayes. However, several literature say that it's good to use samtools to call variant, which was designed for polyploid. I'd like to know should I use the haploid model of gatk and freebayes or just use the default parameters of the three methods to call variant for microbial. I'm not familiar with the algorithm,
anyone1985 is offline   Reply With Quote
Old 08-11-2015, 12:57 AM   #2
SylvainL
Senior Member
 
Location: Geneva

Join Date: Feb 2012
Posts: 179
Default

Hi, do you know how many copies of chromosome your bacterium has? Usually for microbial SNPs/InDels calling, I use freebayes set to haploid but sometimes, like for a Synechocystis project (around 10 copies of the same chromosome), results were nicer considering 10 copies of the genome...
SylvainL is offline   Reply With Quote
Old 08-11-2015, 05:52 AM   #3
anyone1985
Member
 
Location: shanghai, chia

Join Date: Mar 2009
Posts: 67
Default

I do not know the copies of chromosome. My subject is Mycobacterium tuberculosis. In the following papers, they used samtools. I will compare the results from gatk and freebayes (haploid model) and samtools (default parameters). I'd like to know whether there is anyone has done this before.
http://www.ncbi.nlm.nih.gov/pubmed/26116186
http://www.ncbi.nlm.nih.gov/pubmed/25732036
anyone1985 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 12:11 AM.


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