search

nh13 / DWGSIM

Whole Genome Simulator for Next-Generation Sequencing
GNU General Public License v2.0
92 stars 36 forks source link

contig parsing in mut_txt.c required sorted input #64

Closed fbattke closed 3 years ago

fbattke commented 3 years ago

There seems to be a bug in mutation txt parsing in mut_txt.c, line 53 ("i=0;"). This line should be INSIDE the while loop in line 54.., as otherwise parsing only works if contigs in the txt file are sorted identically to those in the fasta reference

nh13 commented 3 years ago

@fbattke I think folks should have the same contig ordering as the FASTA, and so I have updated the error message in #66. I don't see why the two should not be ordered the same (genomic coordinate order). See other formats like VCF or SAM/BAM