I have separate read 1 and read 2 fastq whole genome data and want go for their alignment. Just want to clear, shall I align them separately or first I will concatenate them and then make the alignment together?
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Hi dear,
The data you have is the paired end data, with read 1 and read 2 fastq files are the short read data files with some insert size, most of the times it is 200-400 bps. Ask for the insert size of your data. You can use bowtie2, MAQ or BWA softwares for mapping the data onto the reference genome, choice of mapping software will depend on the read length/insert size. If you want to run bowtie, you may consider the following commands:
You will first need to create a bowtie index of your "reference genome". Run bowtie-build on the genome fasta file:
bowtie2-build Genome.fa Genome.build
This will create six new files which constitute the bowtie index necessary for bowtie:
Genome.build.1.ebwt Genome.build.4.ebwt
Genome.build.2.ebwt Genome.build.rev.1.ebwt
Genome.build.3.ebwt Genome.build.rev.2.ebwt
We can now map paired reads in fastq format to the Genome.build reference sequence:
bowtie2 -S -q --solexa1.3-quals -p 1 -I 100 -X 600 --fr Genome.build -1 read1.fastq -2 read2.fastq OUT.sam
-S: output will be in SAM format.
-q: Quality scores of your data.
-p: Number of processors you want to use.
-l: Minimum insert size.
-X: maximum possible insert size.
-fr: read files are in forward-reverse order.
-1: read 1 fastq file
-2: read 2 fastq file
OUT.sam is the output file
Hope this will be helpful.
Best wishes,
RahulRahul Sharma,
Ph.D
Frankfurt am Main, Germany
-
Yes, there is no need to concatenate the read files, you will lose the pair information. And many of the reads will not map, as they will have insert in-between. Look at the following example:
|=================================| (Reference genome)
(1)-------> <-------- (2)
(1)--------> <--------(2)
(1)--------> <--------(2)
=== is the reference genome.
(1)-----> read 1, your read 1.fastq file will contain all the (1)-----> reads.
<-------(2) read 2, your read 2.fastq file will contain all the <------(2) reads.
Please refer these terms: Paired end data, mate pairs, insert size, read quality scores, read coverage. It would help you in your analysis.
Best wishes,
RahulRahul Sharma,
Ph.D
Frankfurt am Main, Germany
Comment
-
Oh sorry my mistake, from your last post I thought you wanted to actually physically join each read.
If you have many 30 fastq files then yes you can just concatenate them together to create one fastq read library for each read pair.
ie.
Code:cat *R1.fastq > all.reads.R1.fastq cat *R2.fastq > all.reads.R2.fastq
Comment
-
Hi,
Ya you can concatenate the files, but not reads. Please check the Id's in all of your files,
they should be unique before "cat" command.
You may also consider some data preprocessing methods: To trim the adapters and primers, discard low quality reads(with its pair), discard reads with more than 5%-10% N's in it. This can improve your analysis. Please check the quality of your reads with FASTQC or FASTX tools.
Best wishes,
RahulRahul Sharma,
Ph.D
Frankfurt am Main, Germany
Comment
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...-
Channel: Articles
Yesterday, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
55 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
51 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
45 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
55 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Comment