Inserting N bases at 5' end of clipped reads
0
0
Entering edit mode
3.0 years ago

I have been dealing with some clipping of sequencing reads such that the reads align nearly identically, but some reads are trimmed at the 5' end. I want to be able to insert Ns at bases that are missing alignment at the start, what is the best way to do this?

    Seq 1: ATGCTAG
    Seq 2: TGCTAG
    Seq 3: GCTAG

Becomes:

    Seq 1: ATGCTAG
    Seq 2: NTGCTAG
    Seq 3: NNGCTAG

During processing, I have been using Samtools / Bowtie2 / cutadapt and several R libraries, so a solution using one of them would be ideal but open to other suggestions.

GenomicAlignments Samtools cutadapt Bowtie2 • 429 views
ADD COMMENT
0
Entering edit mode

out of curiosity, why would you want to do that ?

ADD REPLY

Login before adding your answer.

Traffic: 2378 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