dwgsim_eval: read not generated by dwgsim?

Question

[fiona.newberry]

I am trying to evaluate mapping accuracy of BWA (by using false positive and true positive).

I want to use dwgsim_eval to get this information from my files. Here is the link to the manual: https://github.com/nh13/DWGSIM/wiki/Evaluating-Mappings-from-DWGSIM

It states if the alignments were made from BWA to use "-b", which I am using but I am getting an error message.

Here is my code:

dwgsim_eval -b sorted_sub_human_BWA.bam

And here is the error message:

0************************************************************
In function "process_bam": Fatal Error[OutOfRange]. Variable/Value: simulated.7796575.
Message: [dwgsim_eval] read was not generated by dwgsim?.
 ***** Exiting due to errors *****
************************************************************

Not sure if there is maybe something wrong with my files? Any help would be greatly appreciated.

[Pierre Lindenbaum]

how did you generate sorted_sub_human_BWA.bam ?

[fiona.newberry]

Using this code:

for filename in ./seqtk_2/subsample_1/*_1.fq.gz;
do file=`echo $filename|sed 's/_1.fq.gz//'`;
filenopath=`basename $file`;
outputpath=./BWA/seqtk_2/subsample_1;
bwa mem -v 3 ./combine_reference.fa.gz ${file}_1.fq.gz ${file}_2.fq.gz | sambamba view -S -f bam /dev/stdin | sambamba sort /dev/stdin -o ${outputpath}/sorted_${filenopath}_BWA.bam;
done

Answer 1

[Pierre Lindenbaum]

how did you generate /seqtk_2/subsample_1/*_1.fq.gz ?

(reads must be generated by dwgsim.... )

This utility evaluates the mappings from reads produced by dwgsim

[fiona.newberry]

I used a simulator called Mason. I see the problem now. I was assuming it meant reads could be simulated by any program but have to specify if it generated by BWA. Are there any other programs that I can use to determine false positive and true positive rate? I want to make roc curves for several different aligners I am trying

[Pierre Lindenbaum]

why don't you just use the fastq generated by dwgsim ?

[fiona.newberry]

I selected the simulator based on the features. I am building an artificial gut community model (with 1000's of sequences) and dwgsim did not have some of the features I needed