Hi all,
I have checked my RRBS raw data (mouse DNA digested with MSPI and bisulfite converted) with fastQC and I have observed a decrease of the base sequence quality on position 2, it's occur in all of my 20 samples R1/R2, there is probably a link with the % of N content on position 2 also, but what can be the experimental reason of this decrease ? Has anyone obtain this kind of results with FastQC and have an idea ? Did there is something to do to improve the quality of my future rrbs ?
thanks a lot for your help
the link to the fastqc reports :
https://drive.google.com/file/d/1vIScAT6BW8TOAYLtoRZ5ihNFqMDsxX5v/view?usp=sharing https://drive.google.com/file/d/1NiBsoihfqp99kU3k6ZHtaDWunS88R0Iq/view?usp=sharing
1 answer
Setting aside rrbs part this may simply be a problem with cycle 2 for this particular run e.g. a bubble went through the lane. Assuming the rest of the data looks ok (sorry not going to click on a random google drive link) you may simply want to hardtrim the first two-three bases. Alignment may also take care of soft-clipping these bad bases.
Log in to answer this question.
-
Dramatic decrease of mapping efficiency between Bwameth and Bismark
written by Raphael •
-
About R question: how to analyse the correlation between two datasets (corr.test and corrplot)
written by Angelrica •
-
How to calculate the weighted average for each linkage group?
written by wkmustahs21 •
-
Adapter sequences in RNA-seq
written by artgrigorov •
-
How to remove illumina universal adapter
written by luckysardar171 •
-
Run ASTALAVISTA for large number of files
written by waqaskhokhar999
-
kmer content failling with fastqc
written by slimane.khayi •
-
FASTQ- Reverse reads of paired end sequence worse than forward reads
written by dilawerkh4 •
-
Short homopolymer runs overrepresented at 3' ends of R2 Illumina RNAseq reads, breaking tophat alig…
written by ethan.kaufman
-
Sudden quality drop in the middle of HiSeq R1 reads but not in R2?
written by sentausa
