Hi
I wonder whether anyone has a solution for a restriction digestion script that works on fastq files (Illumina reads), splitting any reads that contain the restriction site into two separate reads (and writing it back into the fastq file)?
Actually I am not so much interested in restriction as such - I have reads from co-ligated and nebulized PCR fragments, so I'd like to separate those reads at the ligation junctions, and I thought I could tweak a restriction script (preferentially, Perl) into doing that.
Any hints appreciated.
regards,
Berthold
I wonder whether anyone has a solution for a restriction digestion script that works on fastq files (Illumina reads), splitting any reads that contain the restriction site into two separate reads (and writing it back into the fastq file)?
Actually I am not so much interested in restriction as such - I have reads from co-ligated and nebulized PCR fragments, so I'd like to separate those reads at the ligation junctions, and I thought I could tweak a restriction script (preferentially, Perl) into doing that.
Any hints appreciated.
regards,
Berthold
Comment