freebayes error variant calling
0
0
Entering edit mode
7 months ago
evelyn ▴ 140

Hello,

I am using freebayes for variant calling using:

list=list.txt
freebayes  -f  transcripts.fasta -L $list 
          --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 > result.vcf

But it gives error:

ERROR(freebayes): multiple samples (SM) map to the same read group (RG)

I have used picard to add readgroups to the bam files. Thank you!

snp • 291 views
ADD COMMENT
0
Entering edit mode

As the error mention, your BAM file is listing samples with the same read group, which is not possible, you need different read group per sample.

ADD REPLY
0
Entering edit mode

I checked the readgroups for the two samples it is showing error for using:

samtools view -H sample.bam | grep "SM:"

It is showing different names for both the files based on their respective sample names.

ADD REPLY
0
Entering edit mode

SM is sample, not read group.

ADD REPLY
0
Entering edit mode

Thank you! How can I add the readgroup information to bam files using picard. Where can I find the readgroup information for the data I have.

ADD REPLY
0
Entering edit mode

OP is going off of the command I mentioned in their last post (the command was for a totally different purpose)

ADD REPLY
0
Entering edit mode

You should have different ID: entries for each SM: entry. Did you check that? The file name doesn't matter.

ADD REPLY
0
Entering edit mode

what is samtools view -H sample.bam | grep "RG:" showing?

ADD REPLY
0
Entering edit mode

It shows nothing for read group.

ADD REPLY
0
Entering edit mode

It would be @RG, not RG:.

ADD REPLY
0
Entering edit mode

I understand I need to add read groups information. But I could not find where to look for this information. How to check what is RGID, RGLB, RGPL, RGPU for my samples in order to add this information using picard.

ADD REPLY
0
Entering edit mode

Here is some information about what those variables mean. You need to come up with those yourself from information you know about your data. It should be straight forward to do. Ask if you have questions.

Here is the picard command line example of how you would actually add them to your files.

ADD REPLY
0
Entering edit mode

Thank you! I got the fastq files from a collaborator. I don't have information about these variables. What can be the best way to find this information?

ADD REPLY
0
Entering edit mode

Check the first link. Most of that information is going to come from fastq headers/file names etc.

ADD REPLY
0
Entering edit mode

I have checked the fastq file header:

@S79.1 1 length=150
CAAAAGATTATTTTAGAAAGCACAATGTGAGGACAGAACTCTACAAAGTGCTGAATGCAAAGAAATACAATTCTGCAGATCCAATTCCATTAAGATTTAAGTGATGCCCCTACTTTTCCAAGGTCCTCAAGTTCCTTTGAGATTTGCTCC
+S79.1 1 length=150
A<,<FFF<A<KFKKKAFAK<,,AFKFF<FKKKFAA,7FKFKAFKKKA7AAFK7AAK7AKFFFKKFKFKKKKKAKKKKAFKKKKKK,F<,A<<<,,F<KKK7,,,K,A,,F7AFFAAF7<7<F,,A7F7A<,,,,,<FF,7,7,FA,A,<A

The fastq file names were based on sample names when I got those.

ADD REPLY
0
Entering edit mode

You are going to have to make up some of this information since the fastq headers appear to have been rewritten. Make sure you do it right and use reasonable values. Some would be obvious. SM=Real_SampleName, PL=Illumina etc.

ADD REPLY
0
Entering edit mode

I have changed the SM but I am not sure how to change ID which freebayes is complaining about.

ADD REPLY
0
Entering edit mode

Give it a random value, or even set it so ID=SM. The requirement is that multiple SMs don't share the same ID.

ADD REPLY
0
Entering edit mode

evelyn, you have deleted multiple posts in the past, after you've received comments on them. That is bad etiquette. If that behavior is seen again, your account will be suspended.

You also have a lot of posts that have answers (some of which I am sure work and solve your problem), but none of them have been accepted. You either abandon your posts or delete them - that is an abuse of the site. Please go back and provide feedback on them. Pierre asked you to do this before, but you ignored that. Please consider this your final warning.

ADD REPLY

Login before adding your answer.

Traffic: 2429 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