Quantcast
Channel: Post Feed
Viewing all articles
Browse latest Browse all 3231

join non-overlapping paired-end reads

$
0
0

I have a batch of illumina paired end reads with a shorter than anticipated reverse sequence for a ~550bp amplicon.

The joining tools which I have used before (i.e. USEARCH, ea-utils) talk about merging the overlaps and discarding other sequences but is there a way to join forwards and reverses leaving a gapped interval?

Expected output:

>seq1

ATGCAGTCGATCGATCGACTGCATGCATCGATCGAATCGATCG--------------TGCAGTCGATCGTACG

>seq2

CAGTCTAGCTACGATCGATCGACTGCATACGTACGTACG-------------------CTGCATGCATCGTAG

>seq3

CTAGACTGATCGCTAGCCTAGCTACGATCGATCGACACTGCTGTCGACTGACTGATCGATCGATCGACTGCAT

 

Thanks


Viewing all articles
Browse latest Browse all 3231

Trending Articles