SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Illumina insertions deletions andreitudor Bioinformatics 2 02-14-2012 06:20 AM
Novoalign parameters madsaan Bioinformatics 1 10-03-2011 05:34 PM
PubMed: Strategy for Robust Detection of Insertions, Deletions, and Point Mutations i Newsbot! Literature Watch 0 03-01-2011 02:30 AM
Coverage and deletions Hena Bioinformatics 0 01-05-2011 10:05 AM
SOAP2 and deletions/insertions joa_ds Bioinformatics 7 09-21-2010 03:59 AM

Reply
 
Thread Tools
Old 05-27-2011, 12:38 PM   #1
12jrowley2
Member
 
Location: Utah

Join Date: Dec 2010
Posts: 14
Default setting parameters for Novoalign 2-6 bp insertions deletions

I am looking for insertions 2-6 bp in length in the native novoalign output - but see that novoalign did not output any insertions greater than 1 bp. Am I looking for this wrong?

awk '$14 ~ /\+([AGTC]){2,}/ { print $0}'
None found
awk '$14 ~ /\+([AGTC])/ { print $0}'
Lots found

these are the changed parameters used in novoalign with 50 bp single end illumina reads:

-p '-F ILMFQ -t60 -r0.2 -q5' -l -w 24 -k 1 -t -o kaplan-em

am I missing them with these parameters?

thanks for any input!
12jrowley2 is offline   Reply With Quote
Old 10-31-2013, 09:48 PM   #2
sparks
Senior Member
 
Location: Kuala Lumpur, Malaysia

Join Date: Mar 2008
Posts: 126
Default

This works...

awk '$14 ~ /\+([AGTC][ACGT])/ { print $0}'
sparks is offline   Reply With Quote
Old 11-01-2013, 08:48 AM   #3
SNPsaurus
Registered Vendor
 
Location: Eugene, OR

Join Date: May 2013
Posts: 451
Default

I think your -t parameter is too low. According to http://www.novocraft.com/wiki/tiki-i...gnment+Scoring there is a gap opening penalty of 45 and 10 for each base. Given normal quality scores those penalties would pass above -t60. And any additional mismatches would make it worse. Is there a reason you chose a low -t parameter? Just get rid of it to use the default if there is no particular reason. A threshold of 140 would allow a 6 bp indel plus one other mismatch (or seq error).
__________________
Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com
SNPsaurus is offline   Reply With Quote
Reply

Tags
indel, microindel, novoalign

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:02 AM.


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