SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Best way to detect SNV / InDels against reference genome? camhabib Bioinformatics 2 10-22-2015 11:47 AM
Calling Indels IonTom Ion Torrent 2 05-03-2014 07:11 AM
Calling SNPs and InDels simultaneously or seperately? skblazer Bioinformatics 0 08-16-2012 10:31 AM
SNP calling from a reference sequence blackrabite Genomic Resequencing 2 05-21-2011 08:48 PM
Calling indels from samtools-pileup kakar_nipun Bioinformatics 2 01-21-2011 02:31 AM

Reply
 
Thread Tools
Old 01-26-2016, 07:36 PM   #1
sjuzhet
Junior Member
 
Location: Australia

Join Date: Jul 2014
Posts: 9
Default Calling indels against a very small reference sequence?

Hi y'all,

So I'm trying to detect a CRISPR-induced indel in a specific region of the genome, so I only sequenced that particular region and I was planning on using that specific region as a reference sequence, rather than the whole genome. However, most workflows I've researched involved calling indels against large genomes - which is probably gross overkill for my application.

Does anyone have any thoughts on what would be the best way to detect what percentage of my reads contain an indel within a 500-bp region?
sjuzhet is offline   Reply With Quote
Old 01-27-2016, 04:38 AM   #2
HESmith
Senior Member
 
Location: Bethesda MD

Join Date: Oct 2009
Posts: 506
Default

It would be easier to provide suggestions if you described 1) how you selected the region of interest (amplicon? oligo capture?) and 2) the sequencing platform and type of data you generated (Sanger? Illumina PE-250?).
HESmith is offline   Reply With Quote
Old 01-27-2016, 05:43 AM   #3
sjuzhet
Junior Member
 
Location: Australia

Join Date: Jul 2014
Posts: 9
Default

They're PCR amplicons, yes, and sequenced on a MiSeq. Several different cutting conditions were multiplexed, and I was hoping to compare them all against the reference sequence.
sjuzhet is offline   Reply With Quote
Old 01-27-2016, 05:45 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,059
Default

Try BBMap and align to the region of interest (you will need to make BBMap index). Using appropriate maxindel= (and any additional parameters that make sense).
GenoMax is offline   Reply With Quote
Old 01-27-2016, 06:01 AM   #5
HESmith
Senior Member
 
Location: Bethesda MD

Join Date: Oct 2009
Posts: 506
Default

Or, assuming the paired-end reads overlap, you could use BBMerge (part of BBMap) and filter by sequence length. Any merged reads that are shorter than the predicted amplicon length contain deletions.

But BBMap alignment is fairly quick and easy to use, so that's a valid option.
HESmith 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 04:15 PM.


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