Entering edit mode
7.6 years ago
rrsowmya
▴
20
Hi,
I am trying to run this perl script - http://www.pnas.org/content/104/9/3318.abstract?tab=ds with an input file that looks like this ta_iwgsc_7al_v1_4507683 2123 1 GTCCGGAGACGCCGGCG I have tried running this several times but the all_sRNA21.... file has same p value for all reads. This is not possible unless the program doesn't recognize my sRNAfile as a valid one. The first column is very different from the one shown in the script by the inventor which is "Chr". Is that why my output file with p0.001....is empty and I egt same pvalue for all reads? On running the script with a mock test file - like the one the inventor indicated - I am able to get reads that are phased.
Thank you very much for your help.
The link to the perl script is here http://www.pnas.org/content/104/9/3318.abstract?tab=ds
Can you share the input file? Also, if you can explain what are you trying to fetch using the script, may be I can help you.