SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
coveragebed output result problem edge Bioinformatics 2 07-20-2014 11:46 PM
Some problem in the output of CNVnator yuhao Bioinformatics 1 07-08-2012 11:31 PM
Problem with Cufflinks output admiral RNA Sequencing 0 06-01-2011 02:19 PM
SHRIMP 2 SAM output ilyagoldin Bioinformatics 1 10-20-2010 06:07 AM

Reply
 
Thread Tools
Old 10-09-2012, 11:08 AM   #1
rahilsethi
Member
 
Location: Pittsburgh, PA

Join Date: May 2010
Posts: 22
Unhappy SHRiMP 2 output problem

Hi,

I used SHRiMP 2.2.2 for mapping (gapped mapping) SOLiD Paired End data to hg19 genome reference. The SAM output file does not have any value at QV column, instead has "*" as a placeholder. For example:

271_822_1591_F 163 chrM 1 255 35M = 146 194 GATCACAGGTCTATCACCCTATTAACCACTCACGG * AS:i:310 NM:i:0 CS:Z:T12321112012233211002330301311222130 CM:i:2 XX:Z:GATCACAGGTCTATCACCCTATTAACcACTCaCGG
86_1077_1408_F 163 chrM 1 255 3H32M = 126 175 GATCACAGGTCTATCACCCTATTAACCACTCA * AS:i:320 NM:i:0 CS:Z:T33102321112012233211002330301011221 CM:i:0 XX:Z:GATCACAGGTCTATCACCCTATTAACCACTCA

I understand this because SHRiMP 2 never asks for .qual filename/s when put on run. Now almost all the variant calling algorithms such as SNVer, VarScan, GATK etc. require base/color quality values for evaluation of their variant call and end up giving error. Even the mpileup/bcftools seem to be giving abnormal results with almost all of the variant calls as In-Dels. Is there anyway I can make SHRiMP2 give out value to QV column in output so that I can run the variant calling tools on it?

My SHRiMP 2.2.2 run was:

gmapper-cs \
--max-alignments 800 \
-N 12 \
--pair-mode opp-in \
--isize 0,250 \
--all-contigs \
--single-best-mapping \
--local \
--sam \
-1 trimmed_read_F3_50bp.csfasta -2 read_F5.csfasta \
hg19.ucsc.fa > output_trimmed_50bp.sam 2> output_trimmed_50bp.log

Thanks in advance
rahilsethi is offline   Reply With Quote
Old 10-09-2012, 04:06 PM   #2
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

The *csfasta files do not have any quality information, you probably need to input the *qual files as well.
nilshomer is offline   Reply With Quote
Old 10-10-2012, 06:05 AM   #3
rahilsethi
Member
 
Location: Pittsburgh, PA

Join Date: May 2010
Posts: 22
Default

Hi nilshomer,

Thanks for the suggestion. Although SHRiMP, in general, does not give any separate option for .qual files in the run, I did try to incorporate qual files along with csfasta; but it prompts and error, indicating that it does not recognize the .qual format.
One suggestion has been recommended to me by SNVer authors is to include a dummy value in QUAL field of SHRiMP SAM output such as "IIIII" (size of SEQ length) for all mapped reads so as to generate variant calls without runtime error. But then, that would affect the filter process of the variants.
I have also come across the variant calling tool VARiD that shows an example of working on SHRiMP mapping output. I might give a try to VARiD. Do you know how good is VARiD tool for variant calling?
rahilsethi 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:27 AM.


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