This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Mira Tcs Output Currently Not Available

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
• 3,360 views • link
updated by Biostar Biostar • written by diltsjeri diltsjeri

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

• 0 views • link
written by Vikas Bansal Vikas Bansal

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

• 0 views • link
written by diltsjeri diltsjeri

3 answers

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.

• 0 views • link
written by Leonor Palmeira Leonor Palmeira

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!

• 0 views • link
written by diltsjeri diltsjeri

yes. He emailed me back and said that the feature wasn't available yet.

• 0 views • link
written by diltsjeri diltsjeri

Log in to answer this question.

More posts like this