one question in GATK

[evonne16]

hi，
There ‘s error message ,when I use GATK. Please give me some advice,thanks

ERROR MESSAGE: SAM/BAM file SAMFileReader{/home/evonne/work/exome_data/result/GATA1_3_1.sorted.bam} is malformed: Read HWI-ST1106:141:C09YVACXX:3:1308:12485:83335 is either missing the read group or its read group is not defined in the BAM header

[dpryan]

You can use picard tools to add the read group information to reads lacking a read group.

[evonne16]

thanks,but i am a freshman,so how will i know what information of read group is i need?

[dpryan]

Well, the information in the read group describes the sample, the platform used to sequence it, as well as a library identifier. You should be able to ask your grad student/post-doc/PI (or whomever you're directly working under) what he or she would like listed there.

The utility of read groups becomes more apparent with some downstream analyses, particularly if multiple files are merged together. For example, if you need to remove reads which are PCR duplicates, it's nice to have a designation of which library individual reads belong to. GATK can make use of some of the read group information to structure it's analyses, which makes life easier for us users!

[evonne16]

thanks so much for your help！