This is a beta test.
Question: best practices to use GATK with PacBio
1
Entering edit mode

My question Is some how related to this previous post with some differences

Order Of Gatk Commands

I have 2 samples each contains 4 libraries and each library formed of 8 illumina lanes,

add to this that also each sample have PacBio library.

my question is,

How can I make use of PacBio reads:

  • shall I treat the pacbio library like any other library except for the step of Realign Indels I shall skip it in the pacbio reads, but I shall do it for the rest of libraries normally!
  • Or correct the pacbio reads with illumnia for example, then use only the corrected pacbio reads with GATK?

regarding if I will treat PacBio like just another library shal I use the corrected PacBio reads Or the subreads

any previous experience Or Idea?

ADD COMMENTlink 3.6 years ago Medhat 8.3k • updated 14 months ago prasundutta87 • 330
1
Entering edit mode

GATK has not been developed keeping long reads (such as PacBio) in mind. It should only be used for short reads (such as illumina). For variant calling from PacBio reads, its better to use their own variant caller- GenomicConsensus (https://github.com/PacificBiosciences/GenomicConsensus)

ADD COMMENTlink 14 months ago prasundutta87 • 330
Entering edit mode
1

The question was around 3 years ago, there are tools now available for LR. This is an example: https://github.com/aquaskyline/Clairvoyante .

Thanks for suggestion

ADD REPLYlink 14 months ago
Medhat
8.3k

Login before adding your answer.

Powered by the version 1.6