Cannot open configure file in msisensor-pro 1.2.0
Dichromatism opened this issue · 3 comments
Hi,
I have tried msisensor-pro (both installed by conda and pulled from docker). When I was using hg19 as a reference and tried to make a baseline using my .bam files,
docker run pengjia1110/msisensor-pro msisensor-pro baseline -d microsatellites_new.list -i normal_2.tsv -o baseline
the software returned following errors:
`Check for the environment ...
Current work path: /
Microsatellites file: /microsatellites_new.list
Configure file path: /normal_2.tsv
Output path:/baseline
/baseline/ is not existing! now make it! OK!
/baseline/detail is not existing! now make it ! OK!
Load files ...
Open configure file failed, please provide valid configure file ! `
I am sure that my .bam files were indexed correctly and the names were *.sorted.dedup.bam
and *.sorted.dedup.bam.bai
.
When I using version 1.0.2 it returned as cannot find the index files.
So what should I do to figure this out? Thanks!
when you use docker version, you must use -v option in docker run to convert you local path to docker environment path, I will give you a small expamle.
In you computer, suppose you put your bam file at /home/user1/data/baselinebams/,
suppose you have two bam file:
/home/user1/data/baselinebams/sample1.bam and /home/user1/data/baselinebams/sample2.bam
the microsatellie list: /home/user1/data/baselinebams/microsatellites_new.list
write you configurefile (/home/user1/data/baselinebams/normal_2.tsv) as follow:
sample1 /output/sample1.bam
sample2 /output/sample1.bam
PS: the prefix /output is at the docker environment
now you can run msisensor-pro with docker as follow:
docker run -v /home/user1/data/baselinebams:/output pengjia1110/msisensor-pro msisensor-pro baseline -d /output/microsatellites_new.list -i /output/normal_2.tsv -o /output/baseline
Thanks to your reply!
I ran msisensor-pro successfully using pro
method. The baseline were constructed using 33 normal bams. I tested it using one of my paired data. When I using msisensor msi
(is this one equal to msisensor-pro msi
?) to check msi sites with my normal.bam and tumor.bam, the results returned as 76 out of 196 sites were MSI. However when I using msisensor-pro pro
with my tumor.bam and the baseline which I just constructed, the results returned as 5 out of 180 sites were MSI. Which of those results should I believe?
I noticed the msisensor-pro pro
had different defaults in -m
and -w
to msisensor msi
, so I changed those values but the results were same as before.
Hope to your further reply!
I thinks it may be caused by the threshold, you can try to change the -f in msi command. You can also plot a scatter diagram to see if these two results are linear correlation?