Aligning Reads with Overhangs on the Ends (Edges) of Reference Contigs or Scaffolds
0
1
Entering edit mode
4.9 years ago
FatihSarigol ▴ 240

Hi all!

Does anyone know any option in any software that could say:

"In the case a read aligns to the beginning or end of a contig/scaffold of the reference but has an overhang, as long as the matching part is longer than some threshold, do not think of the overhang as mismatch, and calculate mapping just based on the part that actually aligns"

A short example:

REFERENCE CONTIG=                  AGAGAGAGTTTCCCGAGTATATAGAGACACAGAGTGGTGAGGAT
READ WITH OVERHANG=    TTATAGGTTCCGAGAGAGAGTTTCGC
another read=                                   CGAGTATATAGAGACACA
READ WITH OVERHANG=                                            ACAGAGTGGTGAGGGTGGAGAT
ALIGNED READ WITH 2 SNPS=          AGAGAGAGTTTCGCGAGTATATAGAGACACAGAGTGGTGAGGGT

Many Thanks!

alignment SNP non-model mismatch Assembly • 1.7k views
ADD COMMENT
0
Entering edit mode

Hi Durham,

Such type of reads must be filtered out in the unaligned reads. You can start fishing from there and that is the best starting point.

ADD REPLY
0
Entering edit mode

Hi, thanks, I dont want to filter them out, I want to keep them. Otherwise the mismatch caused by the overhang already excludes these reads from the alignment. I have short contigs in the reference and short reads, and this gives me very bad alignment results at the beginnings and ends of the reference contigs. Thanks

ADD REPLY

Login before adding your answer.

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