How get read group information from fastq files in SRA
1
0
Entering edit mode
3.6 years ago

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?

read group fastq SRA Bwa • 2.2k views
ADD COMMENT
0
Entering edit mode
3.6 years ago
ATpoint 81k

This option is not mandatory, just leave it blank. There is no formal definition of a read group, therefore this information is not in the SRA. Only some tools like GATK actually require it.

ADD COMMENT
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

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.

ADD REPLY

Login before adding your answer.

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