Hey guys,
I am quite new to this field and recently I have been asked for an exercise in which I get two different files from an NGS machine corresponding to the forward and the reverse sequence into a fastq format. The exercise consists of finding the mate reads in both files.
Should I just match the whole name on the different files using the sequence identifier, i.e, @SEQ_ID? Or is it a faster way to do it in which I do not need to match the complete name and only the run ID possibly (even if I do not know how the sequence ID is going to be given)?
•
written
by
artist01 •
1 answer
The mate/paired reads can simply be identified by order. E.g., the 3rd read in the R1 file and 3rd read in the R2 file are mates, same for the 1000th, 15347th read etc. No need to match the read IDs. There is some convention that mates are named readid123/1 and readid123/2, but it is not always the case.
Log in to answer this question.
-
Alignment algorithm
written by artist01 •
-
Splitting multifasta by string match
written by Winston •
-
grep command or to select specific protein ID with sequences from a big file...
written by sherafzalk769 •
-
Extract fasta sequnces not matching a list
written by genomes_and_MGEs •
-
Read a fasta file | Java
written by ozdavidd •
-
filtering some sequences from txt file in python
written by ashkan
-
Paired End FASTQ data; separate files for reverse forward strand reads?
written by thenameisizaak •
-
Fastq File Split Using Corresponding Id Identifier In Tag File, Not Barcode Sequence
written by Tonyzeng •
-
Library, Sample Id Tag In Sam Format
written by Ashutosh Pandey
-
Bwa Alignment Of The Same Bam File However Got Two Different Flagstat Results
written by Houkto •
