Thread: bwa alignement
View Single Post
Old 08-24-2020, 11:32 AM   #1
Junior Member
Location: Europe

Join Date: Jul 2018
Posts: 6
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