I have sequencing data from NextSeq and MiSeq, and the GC content increasing rapidly at 3' end base:
Is there something wrong with library? How did this arise?
1 answer
You didn't mention which one is NextSeq and which one is MiSeq, but I can guess that above is NextSeq and below is MiSeq.
Why NextSeq produces more G?
It's caused by the two-colour chemistry system. Four bases with two different colors:
- Green + Red = A
- Green = T
- Red = C
- None = G
As sequencing goes to last cycles, the signal strength decreases. So more bases will be detected as G incorrectly. That's why NextSeq produces more G.
MiSeq or NextSeq is almost the same, and the figure is coming from fastp exactly, it's good to use and I'll use -t 1. Do you have any idea what does last cycle do?
Hi chen, is this a related issue: Can phasing or pre-phase during basecall cause indel? ?
Thanks for taking a look at least. I figured that the N+1 cycle was related to the calibration of errors in the context of phasing and pre-phasing (I still believe that this is related to this N+1 cycle).
Log in to answer this question.
-
FastQC report
written by Bushra Bibi
-
Increasing per base G content in QC sequencing files
written by robertsr •
-
High A in "Per base sequence content" of fastQC report
written by DVA
-
FastQC report explanation on example
written by dominik.saburo •
-
Content of T was larger than others in NSR-RNAseq data, why?
written by wm
-
SNPedia Database question
written by ahm3dhany •
-
Strange Pattern in Bam File
written by Can Holyavkin
-
Are There Some Nextera Barcodes / Adaptors Left?
written by mazzottidr •
-
What Is The Major Difference Between The Domainogram Images
written by cheater149 •
-
How To Put Bootstrap And Posterior Probability Values Together In The Same Phylogenetic Tree?
written by cacaucenturion •
Do you see this also on miseq? On NextSeq this could mean that you sequenced through the adapter and nothing is left to sequence. It's two colour chemistry, so nothing means G. We see polyG stretches at the end. Can't explain it for miseq though, since AFAIK that still uses 4 colours.
Yes, miseq is the same. Insert sise is about 250, and PE150 is used, as you can see above (maybe not clearly), it is caused by the missing of A rather than ployG, G/C/T all goes up except A.
If the base quality at the 3' end is poor, then do not even consider this finding as valid. Have you produced a FastQC report of the reads? The 3' end of Illumina reads is notorious for being poor quality.
Yes, the quality of last cycle is quite poor, I have posted the quality distribution, do you have any idea what does last cycle do?