Mira Tcs Output Currently Not Available
3
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
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
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
2
Entering edit mode
12.4 years ago

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.

ADD COMMENT
1
Entering edit mode
11.8 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!

ADD COMMENT
0
Entering edit mode

Thanks for the update!

ADD REPLY

Login before adding your answer.

Traffic: 1968 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6