humanlongevity / HLA

xHLA: Fast and accurate HLA typing from short read sequence data
Other
101 stars 52 forks source link

problem when run get-reads-alt-unmap.sh #47

Open Huangpq2019 opened 5 years ago

Huangpq2019 commented 5 years ago

When I run get-reads-alt-unmap.sh in.bam out.bam to process my data, there is an error message:

WARNING: Query ERR091574.99978736 is marked as paired, but its mate does not occur next to it in your BAM file. Skipping. ... WARNING: Query ERR091574.9999745 is marked as paired, but its mate does not occur next to it in your BAM file. Skipping.

sambamba 0.7.0 by Artem Tarasov and Pjotr Prins (C) 2012-2019 LDC 1.13.0 / DMD v2.083.1 / LLVM7.0.1 / bootstrap LDC - the LLVM D compiler (0.17.6)

[bam_header_read] EOF marker is absent. The input is probably truncated. [bwt_restore_sa] SA-BWT inconsistency: primary is not the same. Abort! [bam_header_read] invalid BAM binary header (this is not a BAM file). [main_samview] fail to read the header from "-". sambamba-sort: not enough data in stream rm: cannot remove 'out.bam.full.bam': No such file or directory rm: cannot remove 'out.bam.full.bam.bai': No such file or directory

How can I solve this problem?