SEQanswers (
-   Bioinformatics (
-   -   bwa alignement (

Kaskere 08-24-2020 11:32 AM

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::process] read 136386 sequences (10000131 bp)...
[M::process] 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

GenoMax 08-24-2020 11:54 AM

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".

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

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.