Trimmomatic trimming principle
0
0
Entering edit mode
7.6 years ago
shuchen ▴ 10

Why the whole read has to be clipped when the technical sequence overlap or partial overlap with the 5' end of the read?

Figure 1 illustrates the alignments tested for each technical sequence. The process begins with a partial overlap of the 3′ end of the technical sequence with the 5′ end of the read, as shown in (A). Testing proceeds by moving the putative contaminant toward the 3′ end of the read. In both the partial overlap (A) and complete overlap at the 5′ end (B) scenarios, the entire read will be clipped. If the contaminant is found within the read (C), the bases from the 5′ end of the read to the beginning of the alignment are retained. The testing process continues until only a partial alignment on the 3′ end of the read remains (D).

Figure 1 illustrates the alignments tested for each technical sequence. The process begins with a partial overlap of the 3′ end of the technical sequence with the 5′ end of the read, as shown in (A). Testing proceeds by moving the putative contaminant toward the 3′ end of the read. In both the partial overlap (A) and complete overlap at the 5′ end (B) scenarios, the entire read will be clipped. If the contaminant is found within the read (C), the bases from the 5′ end of the read to the beginning of the alignment are retained. The testing process continues until only a partial alignment on the 3′ end of the read remains (D).

trimmomatic read trimming • 2.6k views
ADD COMMENT
0
Entering edit mode

I am unsure what you mean but if a part of the read is recognized as adapter, any sequence that follows must be part of the sequencing adapter and therefore be removed.

ADD REPLY
0
Entering edit mode

I meant why not only clip off the 5' end that was recognized as adapter sequence and retain the rest of read?

ADD REPLY
0
Entering edit mode

You keep the sequence "left" of it (left=before 5' end of the adapter sequence) but the following should just be the remainder adapter sequence.

ADD REPLY
0
Entering edit mode

Thank you very much for your reply. This is exactly where I got confused. Why the following will be the remainder adapter sequence? Is it possible that a read is consist of the 3' end of the adapter sequence and valid sequence as below?

enter image description here

ADD REPLY
0
Entering edit mode

If the read starts with the adapter sequence, then your insert size is 0 and you had 2 adapters joined together.

ADD REPLY
0
Entering edit mode

So you meant that if after the adapter sequence, it is detected again a match with the technical sequence, the read is supposed to be a adapter dimer and that is why it is removed. Right?

ADD REPLY
0
Entering edit mode

That is correct.

ADD REPLY

Login before adding your answer.

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