How to extract out gene specific reads from concatenated fastq file
2
0
Entering edit mode
9 months ago
Mo ▴ 40

Hi all,

I have expressed a library of a heterologous gene in the HEK293 cell line and performed direct cDNA sequencing using a nanopore direct cDNA sequencing (DCS109) kit. I wanted to know what people do to extract their gene-specific reads from concatenated fastq files.

I am using the blast method but it's giving me some issues, I wanted to know what are the other options.

Apologies for this basic question, your suggestion will be really helpful.

Thanks,

fastq gene-specific RNA-seq nanopore • 553 views
ADD COMMENT
0
Entering edit mode
9 months ago
GenoMax 139k

With nanopore reads you should use minimap2 (LINK) which is the best aligner for nanopore data. Write out the results in SAM format and then extract reads that align to your gene of interest.

ADD COMMENT
0
Entering edit mode

Hi Max, Thanks for your reply. I also use minimap2 but before that, I blast the common region from my gene of interest against all reads in the fastq file. I do this to determine whether the strand sequenced is a forward or reverse strand. Then I convert all the reverse reads to forward reads and concatenate them with the rest of the forward reads. After that, I do minimap alignments. The strand determination and conversion of all reads to forward reads helps me to align them with the gene of interest with barcodes.

Do you think this is also possible with minimap to determine the strand direction? Thanks a lot.

ADD REPLY
0
Entering edit mode
9 months ago

People generally don't parse the raw fastq directly like that. Use some kind of aligner to get assignments for all of the reads in the fastq, then you can pull out what you want.

ADD COMMENT

Login before adding your answer.

Traffic: 1768 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6