Question

Mira Tcs Output Currently Not Available

1

Entering edit mode

12.5 years ago

diltsjeri ▴ 470

Hi,

I'm trying to find coverage and %calls at every bp in my resulting sequencing data. I chose MIRA because of this wonderful TCS file which is suppose to provide all this information, but after running the assemblies my file states "TCS output currently not available, please contact author". I did a mapping assembly. I'm not sure if TCS can be done with mapping assemblies. Has anyone else had this problem? Am I missing an option?

Has anyone tried using caf2tcs to force convert the caf file to tcs? or maf2sam?

Thanks!

mira • 2.5k views

ADD COMMENT • link updated 6.5 years ago by Biostar 20 • written 12.5 years ago by diltsjeri ▴ 470

0

Entering edit mode

I am sure someone can help if you can give details about the version of MIRA and the command you used.

ADD REPLY • link 12.5 years ago by Vikas Bansal ★ 2.4k

0

Entering edit mode

No problem. MIRA 3.4.x . Command-->

mira --project=test --job=mapping,genome,accurate,iontor -MI:sonfs=no -SB:lb=1:sbuip=1:bsn=reference IONTOR_SETTINGS -AS:mrpc=100 >test_log_assembly.txt

ADD REPLY • link 12.5 years ago by diltsjeri ▴ 470

score 2 · Answer 1 · 2012-07-01

Your error message says:

TCS output currently not available, please contact author

Have you tried this? The mira_talk mailing list is very active, and Bastien is always available to provide help on MIRA. Also, are you using the latest MIRA version? This might be already fixed and available in the new version.

score 1 · Answer 2 · 2013-01-24

1

Entering edit mode

11.9 years ago

diltsjeri ▴ 470

Not to bring up an old topic but TCS is still not available. However, after doing some communication with Francisco Martins on the mira_talk mailing list (as suggested above), he suggested I use his script "CAF_to_TCS.py". Worked like a charm!