This is a test version of Biostars. For the public version, visit https://www.biostars.org.
freebayes problem in sample name in vcf

Hello,

I am using following command to create a vcf file with freebayes:

freebayes  -f  transcripts.fasta duplicates.bam --max-complex-gap         1 \
          --theta                   0.01 \
          --ploidy                  2 \
          --min-alternate-fraction  0.2 \
          --min-alternate-count     2 \
          --min-coverage            3 \
          --min-mapping-quality     1 \
          --min-base-quality        20 \
          --report-genotype-likelihood-max \
          --throw-away-complex-obs \
          --throw-away-mnps-obs \
          --throw-away-indel-obs > SNPs.vcf

I created the bam file using:

bwa mem transcripts.fastaR1_paired.fastq.gz R2_paired.fastq.gz > paired_bwa.sam
samtools view -S -b paired_bwa.sam > paired_bwa.bam
samtools sort -o paired_bwa.bam -O bam paired_bwa.sorted.bam
samtools index paired_bwa.sorted.bam

but the vcf file does not show sample name in the header, it shows unknown:

#CHROM  POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  unknown

I want the sample name to show instead of unknown. Thank you!

snp

Does the bam file have sample name in the header?

samtools view -H paired_bwa.sorted.bam | grep "SM:"

It shows, @RG ID:4 LB:lib1 PL:ILLUMINA SM:20 PU:unit1

So the sample name is 20. I wonder why freebayes doesn't pick that up.

By the way, did you delete the post?

Hello evelyn!

We believe that this post does not fit the main topic of this site.

OP opened a new post and did not follow up on this one

For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.

If you disagree please tell us why in a reply below, we'll be happy to talk about it.

Cheers!

0 answers

No answers yet.

Log in to answer this question.