Variant calling with duplicates
0
0
Entering edit mode
3.3 years ago
storm1907 ▴ 30

Hello, I would like to check how viral RNA Seq called variants look before read duplicate removal. However, when I try to filter seq files with duplicates, I only have empty output

Are variant callers supposed to work with duplicates?

Thank you!

next-gen • 858 views
ADD COMMENT
0
Entering edit mode

What are you filtering? Not sure if the question is clear for every one.

But the answer is yes, variant callers can work with duplicate reads. Some times it is the default that they remove the duplicates but the variant callers do that based on a tag in the sam/bam file. If you don't add that tag it will use all reads for sure. Maybe others can give a better answer if you tell which variant caller you are using.

ADD REPLY

Login before adding your answer.

Traffic: 1772 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6