I am currently dealing with the problem of reading in a Fastq-File with "readDNAStringset", trimming the Sequences and then writing them in to a new fastq-file.
The reading of the fastq-file with "readDNAStringSet" is working just fine.
I am then trying to trim a fixed length from the left side of the Sequence (ex.10 Bases). Right now I use "subseq(my_Stringset, start =10)".
But when I try to write the fastq-file it seems like I have not trimmed the QualityScores and get the Error message "'x' and 'quality' must have the same width".
I have been looking for a while now, how to trim the sequence aswell as the Q-Score, but I just can't find a Solution (Which should be also fast, since I am working with very big Files).
Is there an correct way to do, what I am trying to do?
Thank you in advance.
I am using R 4.1.0