Question: Align Reads from PacBio
0
gravatar for pthom010
2.7 years ago by
pthom0100
pthom0100 wrote:

I have a gene of interest (about 3000 bps) in several genotypes of a plant I did PCR on. The PCR products were then sequenced using PacBio sequencing and I received my sequencing reads in a file containing multiple fasta files for each individual genotype. I want to align these reads to get a consensus sequence for the genotype and then compare the consensus sequences with each other. What would be the best program to use and would it be possible to use something with a GUI?

sequencing alignment • 1.5k views
ADD COMMENTlink modified 2.7 years ago by WouterDeCoster45k • written 2.7 years ago by pthom0100
2

Hi, recent issue of nature has very nice article talking about 3rd generation sequencing (e.g PacBio, Oxford nanopore etc. ) and available tools for various bioinformatics analysis. Look in to the Table2 of the article. Probably, you will find the solution.

ADD REPLYlink modified 2.7 years ago • written 2.7 years ago by Chirag Parsania1.9k

+1 for sharing the article Chirag. Thanks!

ADD REPLYlink written 2.7 years ago by lakhujanivijay5.3k
0
gravatar for toralmanvar
2.7 years ago by
toralmanvar900
toralmanvar900 wrote:

You can align your pacbio data with 'bbmap mapPacBio.sh' or 'bwa mem' to reference gene sequence. However, mapPacBio can very well handle raw PacBio reads with high error rates. Have a look at this previous post.

ADD COMMENTlink modified 2.7 years ago • written 2.7 years ago by toralmanvar900
0
gravatar for WouterDeCoster
2.7 years ago by
Belgium
WouterDeCoster45k wrote:

I suggest minimap2: https://github.com/lh3/minimap2

and would it be possible to use something with a GUI?

Nope, but a command line is not that bad and it's better to bite the bullet now and spend some time with it to get comfortable. You'll thank yourself in the future.

ADD COMMENTlink written 2.7 years ago by WouterDeCoster45k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 985 users visited in the last hour