lindenb / jvarkit

Java utilities for Bioinformatics
https://jvarkit.readthedocs.io/
Other
478 stars 132 forks source link

Are there default filters in sam2tsv? (general question; not an issue) #174

Closed kac2053 closed 7 months ago

kac2053 commented 3 years ago

Verify

java version "1.7.0_221" OpenJDK Runtime Environment (rhel-2.6.18.0.el7_6-x86_64 u221-b02) OpenJDK 64-Bit Server VM (build 24.221-b02, mixed mode)

Subject of the issue

Are there default filters in sam2tsv?

We have previously published data using Rsamtools pileup to generate counts for mutations. We are now using sam2tsv to index which read IDs belong to which counts. However when I count up the number of reads to the mutation in the sam2tsv data, the counts are not matching up with pileup counts, despite filtering the sam2tsv data based on pileup default filters. I don't have a bioinformatics background, so I'm sorry if I answered any tech questions below wrong.

Thank you.

Edit: Is base quality in pileup correspond to the read quality column in sam2tsv?

Your environment

$ java -jar bin/jvarkit/dist/sam2tsv.jar --version 2b40bcc06

java version "1.7.0_221" OpenJDK Runtime Environment (rhel-2.6.18.0.el7_6-x86_64 u221-b02) OpenJDK 64-Bit Server VM (build 24.221-b02, mixed mode)

-bash: /opt/common/CentOS_7/java/jdk1.8.0_131: Is a directory

Linux 3.10.0-957.12.2.el7.x86_64 x86_64

Steps to reproduce

N/A because general question.

Expected behaviour

N/A

Actual behaviour

N/A

lindenb commented 3 years ago

Are there default filters in sam2tsv?

there is none while mpileup will remove duplicates, secondary etc...