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.
-
Developments in sequencing technologies and methodologies have transformed the field of epigenetics, giving researchers a better way to understand the complex world of gene regulation and heritable modifications. This article explores some of the diverse sequencing methods employed in the study of epigenetics, ranging from classic techniques to cutting-edge innovations while providing a brief overview of their processes, applications, and advances.
Methylation Detect...05-31-2023, 10:46 AM -
Differential Expression and Data Visualization: Recommended Tools for Next-Level Sequencing Analysis
After covering QC and alignment tools in the first segment and variant analysis and genome assembly in the second segment, we’re wrapping up with a discussion about tools for differential gene expression analysis and data visualization. In this article, we include recommendations from the following experts: Dr. Mark Ziemann, Senior Lecturer in Biotechnology and Bioinformatics, Deakin University; Dr. Medhat Mahmoud Postdoctoral Research Fellow at Baylor College of Medicine;...05-23-2023, 12:26 PM -
Continuing from our previous article, we share variant analysis and genome assembly tools recommended by our experts Dr. Medhat Mahmoud, Postdoctoral Research Fellow at Baylor College of Medicine, and Dr. Ming "Tommy" Tang, Director of Computational Biology at Immunitas and author of From Cell Line to Command Line.
Variant detection and analysis tools
Mahmoud classifies variant detection work into two main groups: short variants (<50...05-19-2023, 10:03 AM