Thanks this answers my question.
• 0 views
•
link
written
by
Priyanka •
Matching R1 and R2 in paired end reads from Illumina
I have a very basic curiosity.
How does tools match paired end reads?
Suppose I have Demultiplexed paired end reads from R2 taking it as single end read and then want to find the matching pairs for all of the read pairs in R1. So are the header information of pairs matching somehow or do the tools that take in paired end reads take reads one after other in the order they appear?
And if I want to manually find read pair of specific R2 read in R1, is there any way to do so?
• 2,093 views
•
link
written
by
Priyanka •
1 answer
I assume you are asking about Illumina data since that is the tech producing paired-end reads now (other technologies have done so but I am not sure how common they are, BGI may be one). You have likely seen the spec for fastq format spec. Read numbers are found in fastq header.
@EAS139:136:FC706VJ:2:2104:15343:197393 1:Y:18:ATCACG
*1*:Y:18:ATCACG - 1 represents read 1.
Most tools may not check for 1:1 concordance of reads in R1/R2 files, most aligners don't. So don't assume tools are paying attention to headers. Tools will generally expect reads to be in sync in R1/R2 files.
Note: If you suspect that your R1/R2 files are not in sync then a tool like repair.sh from BBMap suite can be used to remove singleton reads to bring them back in sync.
You can find specific reads searching your files with zgrep (gzipped)/grep(plain fastq) -A 3 header_ID yourfile_R1/R2.
Log in to answer this question.
More posts like this
-
Merge paired end reads after mapping in BAM files
written by muhammad_elhossary •
-
How to cat 2 paired end reads files together in one file using loops for multiple libraries ?
written by v.berriosfarias •
-
Paired end data with R1 demultiplexed but not R2
written by singlegenome •
-
How to trim both pair end reads and remove low quality reads
written by Bioinfonext
-
Trim Paired-end Fastq Files
written by yuabrahamliu •
-
How to extract equal number of R1 and R2 read from raw Illumina data
written by Bioinfonext
-
Tool to append Illumina Read 1 to Read 2 for downstream demultiplexing
written by achamess
-
GBS paired reads splits it into several smaller files based on barcode matching
written by nakanomasayuki265 •
-
Calculating number of reads for paired end reads?
written by Genosa
-
samtools + map
written by elfeusilva •
