Question: plasmid spades R1/R2 files?
0
gravatar for sarah.goldstein
6 months ago by
sarah.goldstein0 wrote:

Hello,

I know spades innately saves the corrected R1/R2 fastq files for the assembly. However, the plasmid spades program only uses a subset of those reads, so I was wondering if there is a flag or way to pull out just those reads used to build the plasmid assembly and place them in separate R1/R2 files?

I am also open to using other mapping methods (BWA) but I just don't know how to split fastq files into genomic and plasmid DNA.

Thanks!

ADD COMMENTlink modified 6 months ago • written 6 months ago by sarah.goldstein0

If you have the genomic reference available you can use an aligner (look into bbmap.sh from BBMap suite with outu=file.fq.gz option) to collect reads that do not map. Those should be your plasmid reads (assuming there is no other entity expected to be present).

ADD REPLYlink modified 6 months ago • written 6 months ago by genomax80k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 850 users visited in the last hour