Mapping amplicon to a reference sequence with iteration / elongation
1
0
Entering edit mode
22 days ago

Hi, I try to reconstruct a sequence of 700 bp (barcode) with 250x pair-end Illumina reads.

I have several barcodes which are multiplexed in an Illumina Index file. Each final barcode have a length of 700bp. For the sequencing needs, each barcode was aleatory cut. After sequencing, I juste know the tag+primer nucleotide sequences, which permit to find the beginning of my barcode (250bp). Then I would like map all the reads of my Index file to this beginning read.

I currently work with Geneious Prime, which permit to do : "Map to reference with iteration" It's work great, but take a lot of time to do it manually.

Have you an idea of how i can do that ? The only solution, I have found, is to use minimap 2 and export the consensus alignment from sam files to fasta, and do it again until my consensus > 700bp Thanks

Barcoding NGS Mapping • 117 views
ADD COMMENT
0
Entering edit mode
2 days ago

Hi, Finally, I solve my problem with several customs scripts.

ADD COMMENT
0
Entering edit mode

Your problem description is not clearly stated. I have tried to read it a few times but can't figure out exactly what you want/need to do. Since your custom scripts are not going to be available to people who find this thread in future if you are able to provide a step-wise logical description of what your scripts are doing perhaps an alternate solution can be found.

ADD REPLY

Login before adding your answer.

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