Header Leaderboard Ad

Collapse

bcftools - retain low quality read position support?

Collapse

Announcement

Collapse

SEQanswers June Challenge Has Begun!

The competition has begun! We're giving away a $50 Amazon gift card to the member who answers the most questions on our site during the month. We want to encourage our community members to share their knowledge and help each other out by answering questions related to sequencing technologies, genomics, and bioinformatics. The competition is open to all members of the site, and the winner will be announced at the beginning of July. Best of luck!

For a list of the official rules, visit (https://www.seqanswers.com/forum/sit...wledge-and-win)
See more
See less
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • bcftools - retain low quality read position support?

    Is it possible to have `bcftools call -c` output a DP4 values that add up to the total DP value. I want to call simple consensus sequences based on a basic vote.

    For example, at one position in my VCF file, I have:
    DP=16092 and DP4=74,48,24,600 (alt allele should win 624 to 122)
    Those dont come close to adding up.
    Checking out the position in IGV, the ref allele wins 11355 to 3716.

    I did run samtools mpileup -B to turn off BAQ scoring.
    What am I missing here?

Latest Articles

Collapse

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Today, 07:14 AM
0 responses
4 views
0 likes
Last Post seqadmin  
Started by seqadmin, Yesterday, 01:08 PM
0 responses
6 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-01-2023, 08:56 PM
0 responses
110 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-01-2023, 07:33 AM
0 responses
247 views
0 likes
Last Post seqadmin  
Working...
X