Entering edit mode
3.6 years ago
zamani_2012
•
0
Hi everyone, I want to align some paired fastq files to canine reference genome using BWA-MEM tool. I downloaded fastq files from SRA. I need to specify read group information using -R option in BWA-MEM. I dont have such information. How can I access to read group information for each fastq file? Are these information available in SRA? any help will be greatly appreciated?
Many thanks for your comment. If I do not enter this information, will it not have a negative effect on subsequent analyzes? I work on WGS data and I use GATK for variant calling.
If you are going to use GATK then you will have to use read groups. You can make up basic RG information based on type of sample you are working with.