Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Cleaning multiple sequences alignement Flox1991 Bioinformatics 0 04-20-2020 02:27 PM
Mini cluster specifications for sequence alignement Won-Sul Bioinformatics 2 05-02-2016 06:43 AM
calculate alignement differences atma_weapon Bioinformatics 0 04-11-2013 03:22 PM
local mode alignement of reads StephaniePi83 Bioinformatics 0 06-07-2012 07:18 AM
bwa alignement outside of the reference sequence cklopp Bioinformatics 1 09-14-2011 08:25 AM

Thread Tools
Old 08-24-2020, 10:32 AM   #1
Junior Member
Location: Europe

Join Date: Jul 2018
Posts: 7
Default bwa alignement

I am trying to align viral RNA sequences against reference fasta file.
It looks like it is not working properly

my index code looks like that

./bwa index reference.fasta

but I keep getting output like

[M::bwa_idx_load_from_disk] read 0 ALT contigs
[M:rocess] read 136386 sequences (10000131 bp)...
[M:rocess] read 136392 sequences (10000087 bp)...
[M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (0, 24370, 7, 0)
[M::mem_pestat] skip orientation FF as there are not enough pairs
[M::mem_pestat] analyzing insert size distribution for orientation FR...
[M::mem_pestat] (25, 50, 75) percentile: (138, 175, 224)
[M::mem_pestat] low and high boundaries for computing mean and (1, 396)
[M::mem_pestat] mean and (186.45, 62.40)
[M::mem_pestat] low and high boundaries for proper pairs: (1, 482)
[M::mem_pestat] skip orientation RF as there are not enough pairs
[M::mem_pestat] skip orientation RR as there are not enough pairs
[M::mem_process_seqs] Processed 136386 reads in 3.278 CPU sec, 3.070 real sec

Should I include options like -p and -a? I have no idea what exactly should be written in "Prefix of the output database" after index command
Kaskere is offline   Reply With Quote
Old 08-24-2020, 10:54 AM   #2
Senior Member
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,101

What is not working properly? Indexing genome and alignment are two independent steps. You completed the indexing of the genome. Now you need to align using `bwa mem`+ this index + your reads => Aligned SAM file => Aligned BAM file.

You could give your database a name that is what the "prefix" refers to. If you don't provide a name then your index is going to use the same "basename" as your reference file. In this case, "reference.fasta".
GenoMax is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 02:02 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO