Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

sys.stderr.write("Error: " + len(vcf.samples) #3

Open
deb0612 opened this issue Mar 15, 2021 · 3 comments
Open

sys.stderr.write("Error: " + len(vcf.samples) #3

deb0612 opened this issue Mar 15, 2021 · 3 comments

Comments

@deb0612
Copy link

deb0612 commented Mar 15, 2021

Dear sir,
I use the command below:
python bin/pyTMB.py -i WES-N0440vs5.gnomad.g1kpon.vcf --bed "/Everythings/bed/hg38/Roche_KAPA_HyperExome_hg38_capture_targets.bed" --dbConfig config/annovar --varConfig conf/mutect2.yml --polymDb 1k,gnomad > TMB_results.log
and got the error

Traceback (most recent call last):
File "bin/pyTMB.py", line 305, in
sys.stderr.write("Error: " + len(vcf.samples) +
TypeError: must be str, not int

The input vcf was generated from gatk 4.1.9 mutect2
reference genome: Homo_sapiens_assembly38.fasta

nservant added a commit that referenced this issue Mar 15, 2021
@nservant
Copy link
Contributor

Hi,
I made a quick fix.
Please note that we shoul release a new version in the coming days.
Best

@deb0612
Copy link
Author

deb0612 commented Mar 16, 2021

Dear sir,
I had tried the new version. It works well.
But if I use mutect2 tumor-normal output as input of TMB.

$ python bin/pyTMB.py
-i WES-N0440vs5.gnomad.g1kpon.vcf.hg38.filtered.vcf"
--bed "/Everythings/bed/hg38/Roche_KAPA_HyperExome_hg38_capture_targets.bed"
--dbConfig config/annovar_102015.yml
--varConfig config/mutect2.yml
--minDepth 100 --filterLowQual --filterNonCoding --filterSplice --filterSyn --filterPolym --polymDb 1k,gnomad --export >test.log

It showeda message:
Error: 2 sample detected. This version is designed for a single sample !

Here is my input file look like:
image

Is there any suggestion to delete another sample column?

@nservant
Copy link
Contributor

you can try to use the --sample parameter with the name of the sample you would like to focus

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants