Hello all,
I have been doing some sequencing with Nanopore and during the library prep for nanopore, we add in a DNA Control strand (DCS) to each of our samples as a QC to ensure that our library prep is done right. Now that I have my fastq.gz files, I assume the DCS is still in my sequences. What tools can I use to remove the DCS from my fastq.gz files? (the DCS is given as fasta sequences in the nanopore website).
Thank you so much for the help.
•
written
by
janny.lau •
2 answers
Two ways:
a) Run alignment as usual, and add the sequence as an extra chromosome to the reference genome fasta. Then you can later easily remove the reads aligning to that extra chromosome from the bam file. That's probably the (imo) best way as it does not really requires much custom code other than adding the sequence to the genome and index that.
b) Use something like https://bioinf.shenwei.me/seqkit/usage/#grep to find the reads with a match for this sequence. Extract the read names for the matches, and then do some Unix-fu to remove them. Probably this can also be done with some seqkit magic, it's a quite powerful toolkit.
Log in to answer this question.
-
How to normalise Nanopore mRNA sequencing data between two cell lines
written by Mo •
-
Out of Sequencer Nanopore Raw Data
written by SomeOne •
-
For people doing high molecular weight DNA sequencing library prep, do you see a steep drop off in …
written by Mark •
-
Trimming adaptors and primers for RNAseq reads
written by Jennifer •
-
Clipping primer after alignment
written by finswimmer
-
How to remove contamination reads from fastq flle
written by jeccy.J
-
Can we get absolute abundance of TCR using TCR-Seq?
written by CY
-
how to add inline barcode to raw DNA fragment without PCR
written by liaoyunshi •
-
Consequences of degraded DNA for Sequencing
written by firestar
-
Forum: PacBio Library Prep workshop and PacBio SMART-Portal bootcamp - UC Davis - April 2015
written by DNATECH •
