Extension of reads
1
0
Entering edit mode
4.4 years ago
madhu.9124 ▴ 50

Hi, I just want to extend the read by 20bp to the 3' or 5' end of the read, with high quality score, based on the reference genome sequence. Kindly suggest command or tool for doing it.

sequence samtool BBMap • 1.1k views
ADD COMMENT
1
Entering edit mode
4.4 years ago

BBMap does not do reference-based extension, but you can do coverage-based extension like this:

tadpole.sh in=reads.fq out=extended.fq el=20 er=20 mode=extend k=62

ADD COMMENT
0
Entering edit mode

thanks, is there any tool which can do reference based extension. How do i also extend with high quality score?

ADD REPLY
0
Entering edit mode

I don't know of a tool like that.

The quality score Tadpole assigns for extension is arbitrary; I think it's Q30.

ADD REPLY
0
Entering edit mode
How do i also extend with high quality score?

Why are you so interested in getting high quality scores on the alignment (considering questions you have posted so far)? We generally want high Q scores for raw data.

ADD REPLY
0
Entering edit mode

I am testing BBmap application and the requirement says high quality score.

ADD REPLY
0
Entering edit mode

What requirement are you talking about?

ADD REPLY
0
Entering edit mode

the functional requirement to test the application.

ADD REPLY

Login before adding your answer.

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