Merging paired fastq read files with small overlap region.
I have some data from a lab member who sequenced a crispr library plasmid pool that they made. The forward and reverse reads overlap by 6 bp's and I was wondering how I could merge these fastq files together based upon knowing this overlap size. When I use bbmerge or other merging tools, few reads (less than 5%) are merged presumably because of the very short overlap region.
• 1,763 views
•
link
updated
by
Biostar
•
written
by
dbready2 •
•
written
by
dbready2 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
merging .fastq.gz files results in corrupted files
written by rDNA •
-
Novaseq read coverage needed for demultiplexing
written by cassie.bishop •
-
BBMerge / Tadpole error correction
written by bio23 •
-
How Should I Deal with Paired-End Shotgun Metagenomic Reads for DIAMOND Analysis?
written by ian.petersen •
-
Merge paired end reads without overlap
written by ahramkim1128 •
-
bbmerge merged read smaller that original reads
written by snishtala03 •
-
Insert size of ~20bp with ATAC-seq?
written by a.rex
-
Deletions in a consensus sequence.
written by kspata
-
Metabarcoding amplicon size too long for paired-end sequencing
written by lvogel •
-
Should I merge forward and reverse reads before mapping
written by qiyunzhu
Yes, that is in fact a short overlap, and due to the limited size of overlapping bases the confidence to decide whether the overlap is true and reliable is actually limited. Wouldn't it be simpler to trim one of the reads back a few bases?
Giving this a try a now.
let me try to understand, it is always 6bp?
Yes, the amplicon library they prepared is from a CRISPR library pool in which the only thing that varies is what is contained in the 20 bp gRNA. I inspected a couple read pairs in the fastq to confirm.
could you post 5 sequences or so?
If these are amplicons and you are sure they should overlap then try the following option with
bbmerge.sh. You will need enough sequence data for this to work. Set this to 10 and see if that works.