Hi all,
As you may see from the picture I have this QC from all R2 reads of my Paired End sequenced samples. Due to a technical error during the sequencing I am ending up with 30+ R2 reads with serious errors in the middle of the sequence.
Do you know any way to mask (or to allow mismatch at) a specific number of bases (2-3) at a specific position WITHIN the fragment length prior to alignment? Biostrings is an option that I would prefer not to use for reasons of speed.
Can what you propose applied selectively to only one of the two reads in the paired end samples?
It would be ideal if this could be directly applied directly with Bowtie like the trimming left/right that already exists as an inherent option.
As you may see from the picture I have this QC from all R2 reads of my Paired End sequenced samples. Due to a technical error during the sequencing I am ending up with 30+ R2 reads with serious errors in the middle of the sequence.
Do you know any way to mask (or to allow mismatch at) a specific number of bases (2-3) at a specific position WITHIN the fragment length prior to alignment? Biostrings is an option that I would prefer not to use for reasons of speed.
Can what you propose applied selectively to only one of the two reads in the paired end samples?
It would be ideal if this could be directly applied directly with Bowtie like the trimming left/right that already exists as an inherent option.