Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
  • Oct 18 17:59
    muhammed-ali commented #85
  • Oct 09 17:29
    dfermin commented #85
  • Oct 08 13:14
    dfermin commented #85
  • Apr 09 00:16

    BenLangmead on master

    embed parsed_md (compare)

  • Apr 08 21:26

    BenLangmead on master

    make more portable (compare)

  • Apr 08 21:18

    BenLangmead on master

    some attempts to make this scri… (compare)

  • Dec 03 2018 15:25
    gianmaz edited #88
  • Dec 03 2018 14:25
    gianmaz opened #88
  • Mar 15 2018 20:05
    ChristopherWilks opened #87
  • Mar 15 2018 20:01

    ChristopherWilks on master

    switched to use sratoolkit 2.8.… (compare)

  • Mar 04 2018 22:50

    nellore on master

    patches bowtie2-build in travis… (compare)

  • Mar 04 2018 22:11

    nellore on master

    uses bowtie2 2.3.4.1 (compare)

  • Mar 04 2018 21:57

    nellore on master

    specifies samtools version to i… (compare)

  • Mar 04 2018 21:47

    nellore on master

    updates dependencies Merge branch 'master' of https:… (compare)

  • Mar 04 2018 21:39

    nellore on master

    quote rules Merge pull request #86 from Ben… (compare)

  • Mar 04 2018 21:39
    nellore closed #86
  • Feb 04 2018 01:11

    nellore on master

    enables --keep-alive for covera… (compare)

  • Jan 10 2018 13:09
    BenLangmead opened #86
  • Dec 29 2017 05:47

    nellore on master

    fixes unit test in bed_pre (compare)

  • Dec 29 2017 03:11

    nellore on master

    writes coverage tsv, optionally… (compare)

abhinav
@nellore
can you view /disc2/rail-rna_logs/precoverage/dp.reduce.log/6.0.log?
Vy T. Nguyen
@vnguyen0009
That would be great! What do you need from me?
abhinav
@nellore
that log should have more information about the error
Vy T. Nguyen
@vnguyen0009
Traceback (most recent call last):
File "app_main.py", line 75, in run_toplevel
File "/home/ubuntu/raildotbio/rail-rna/rna/steps/coveragepre.py", line 180, in <module>
for
, _, uniqueness, diff in diffs:
ValueError: expected length 4, got 8
abhinav
@nellore
that's no good and somewhat opaque
do you have /disc2/91_190506_TruSeq/scratch -S 307200 -k1,1 -k2,3 -t$'\t' -m /disc2/rail-rna_logs/precoverage/dp.tasks/6.*?
err
that is
/disc2/rail-rna_logs/precoverage/dp.tasks/6.*?
Vy T. Nguyen
@vnguyen0009
yes I have that file.
abhinav
@nellore
how many such files are there?
Vy T. Nguyen
@vnguyen0009
There are seven 6.* files
abhinav
@nellore
are they large?
Vy T. Nguyen
@vnguyen0009
no, between 40-42M
abhinav
@nellore
so those files contain no identifiable info EXCEPT for possible deletion signal -- they simply contain just enough information to construct coverage bigwigs specifying how many primary alignments map across each base of the genome
if you're comfortable sending them to me
i can take a look
i suspect it'll be hard to diagnose this issue without studying those files
Vy T. Nguyen
@vnguyen0009
I think it should be fine. Let me double check with my supervisor really quick.
Yes I should be able too. Sorry, just needed to double check.
Vy T. Nguyen
@vnguyen0009
Trying to send them through gitter. The files are too big for email. ;(
Okay, sent you a link to share the files through onedrive. Let me know if you get it. I sent it to your email.
abhinav
@nellore
i see it!
thanks
taking a look
abhinav
@nellore
@vnguyen0009 i know what's happening, and i can tell you a quick fix. i do not know why it's happening
in a previous step, a newline is not being printed
for whatever reason
if you edit the file 6.1
the line that reads chr2;23995 000119979136 5 1 1chr17;8056 000040283230 7 1 -1
insert a newline between 1 and chr17
then resume the pipeline
it will work
this is an I/O bug that is likely platform-specific
Vy T. Nguyen
@vnguyen0009
Okay thank you! I will do that and let you know what happens!
abhinav
@nellore
great!
Vy T. Nguyen
@vnguyen0009
which platform are you talking about, if you don't mind me asking. My Ubuntu instance?
And if I run into this issue again, would it be the same fix?
Vy T. Nguyen
@vnguyen0009
Whenever I resume flow from the run command give to me in the flow.2019-10-30T13\:56\:17.966831.log, the line I edited goes back to the way it started.
abhinav
@nellore
wow
i'll have to take another look in the next couple days
sorry you're experiencing this issue
Vy T. Nguyen
@vnguyen0009
It's all good! Thanks for taking the time. We'll be trying to do our own troubleshoot so hopefully we can figure it out together!
abhinav
@nellore
thank you!!!
Vy T. Nguyen
@vnguyen0009
Quick question, which script makes the precoverage/dp.task files?
abhinav
@nellore
coverage_pre.py
Vy T. Nguyen
@vnguyen0009
Thank you!
Vy T. Nguyen
@vnguyen0009
Hello! Just running into some issues. We were able to run rail through our data sets but now we need to compile them into gene count tables. Is there an easy way of doing this? I took a look at the tutorial you put a link to on the rail homepage (http://jtleek.com/protocols/bigwig_gene_counts/) but I am having some trouble figuring out how to do it for all chromosomes.
Vy T. Nguyen
@vnguyen0009
Hello, I am still trying to obtain gene counts from my data obtained from Rail-RNA. I managed to convert the BigWig files to Wig files and was wondering if you could clarify what the values in the last column. It would really help me be able to figure out how to get to gene counts.

bedGraph section chr1:0-631636

chr1 0 14678 0
chr1 14678 14689 1
chr1 14689 14690 6
chr1 14690 14692 7
chr1 14692 14697 10
chr1 14697 14754 22
chr1 14754 14765 21
chr1 14765 14766 16
chr1 14766 14768 14
chr1 14768 14771 12
chr1 14771 14772 5