Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Bisulfite sequencing - filtering by min. conversion rate mixter Epigenetics 3 06-09-2012 02:33 AM
bowtie2 read length in score-min okyoon Bioinformatics 0 04-09-2012 09:53 AM
min-frags-per-transfrags option in cufflinks Ramprasad Bioinformatics 0 11-10-2011 04:36 AM
Min FPKM Value for Diff Exp pcg Bioinformatics 0 03-01-2011 12:53 PM
tophat -F/--min-isoform-fraction option wenhuang Bioinformatics 0 03-31-2010 01:02 PM

Thread Tools
Old 06-06-2012, 05:03 PM   #1
Junior Member
Location: Vancouver BC. Canada

Join Date: Jun 2012
Posts: 2
Question --min-avg-qual of VarScan

I'm using VarScan.v2.2.11.

The problem I'm reporting here is something about VarScan's parameter --min-avg-qual when using command mpileup2indel. I'm setting it to be default, which is 15.

In the input pileup file:
CHROMOSOME_I 221154 g 42 .,,,,,,,+1c,+1c,,+1c,+1c,+1c,+1c,+1c,+1c,+1c,+1c,+1c,+1c,+1c.+1C,+1c.+1C,+1c,+1c,+1c,+1c.+1C.+1C,+1c,+1c,+1c.+1C,+1c,+1c,+1c.+1C.+1C.+1C.+1C.+1C 1C?C<FFFFDHJEAJDJIJG<JJG'GIDIJBD9JDDDG?FDF

The corresponding output:
CHROMOSOME_I 221154 G +C */+C:41:8:34:45.33%:4.4044E-13 Pass:11:30:10:24:1E0 0 1 0 0 */+C:41:8:34:45.33%:4.4044E-13

The red colored part of the input won't get passed the --min-avg-qual filteration, because it's base quality is less than 15. The prblem is that in the output, we see the Cov=41 (which is correct), Reads1=8 (which is also correct), but Reads2=34 (which should be 33, because one of the insertion has mapping quality less than 15).

Anyone has a similar problem like me? Could anyone help? I think this might be a bug of VarScan.
AlbertZANG is offline   Reply With Quote

base quality, mpileup, varscan

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 05:56 PM.

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