Magdoll / MagPhase

Phasing for metagenomics using PacBio long reads
BSD 3-Clause Clear License
17 stars 3 forks source link

Reads lost after processing #2

Open manuelaraujo1000 opened 2 years ago

manuelaraujo1000 commented 2 years ago

Hi,

I seek help understanding the final output. I have processed a file with around 2000 sequences, with the objective of locating a SNP in a 10 nucleotide region. The software was perfect to locate this SNP, however it seems some reads which spanned this region have been dropped during the analysis. Attached I send two of the output files the pos.txt has the correct number of reads, however in the hap.txt these are reduced and for the read.txt this is the same. I would like to understand why this is happening.

Thanks in advanced. Capture image