KChen-lab / Monopogen

SNV calling from single cell sequencing
GNU General Public License v3.0
68 stars 16 forks source link

nothing in the germline.vcf file #21

Open hanliu0119 opened 9 months ago

hanliu0119 commented 9 months ago

Case: using our bam file, the result of char20.germine.vcf contains nothing Description: To test your code, we used the same region.lst under the ./example folder, and used our own bam file in the bam.lst file. but the result is not expected. Please let us know if there is anything that we can improve to work it out.

Here are the additional information you might want to check out.

  1. the input
  2. the output

I appreciate it a lot if you would like to take a look

jinzhuangdou commented 9 months ago

Could you check whether there are reads in your input bam file? You can use the whole sequencing file from chromosome 20 (our example only includes reads from chr20:0-2Mb)

hanliu0119 commented 9 months ago

I checked the bam file, it did have reads and we also used the whole sequencing file. Would you mind helping us check why our bam has no output?