Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Low mapping percentage of reads on assembled contigs morning latte Bioinformatics 20 03-23-2014 10:01 AM
Why is bcftools missing back-to-back/high-divergence-region SNPs (BWA mapping)? Genomics101 Bioinformatics 1 07-23-2012 02:08 PM
Too few reads mapping back to contigs cerebralrust Bioinformatics 8 02-26-2012 11:25 PM
question about mapping reads back to an endosymbiont genome NGS_Newbie RNA Sequencing 0 09-26-2011 07:20 PM
mapping transcriptome reads to assembled ESTs MadraghRua Bioinformatics 0 06-28-2011 09:54 AM

Thread Tools
Old 04-04-2016, 09:35 AM   #1
Senior Member
Location: USA

Join Date: Nov 2013
Posts: 182
Default Mapping assembled genome back to reads: why SNPs?

Hi members,

I've paired end, Illumina, Klebsiella Pneumonia.
Read length 251, MiSeq run
Quality remains good until the read end.

Assembled genome's quality is reasonable checked with Quast.

When I map my reads back to assembled genome, and carry out SNP analysis, I find SNP within it.
I set coverage as 10 to call a base SNP.
I fail to understand why a reference has SNP when mapped against itself.

Is it something to do with SPAdes' correction?

Tool used for mapping, and VCF - bowtie, and GATK.
spades version: 3.5.0

Spades run on auto coverage settings, everything else default. Bowtie, and GATK default.

I found similar post at:

Any help shall be appreciated.
bio_informatics is offline   Reply With Quote

mapping, spades

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 09:47 PM.

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