Question: Bug in combine.py for RNA-Seq analysis
0
gravatar for rvjensen
24 days ago by
rvjensen10
rvjensen10 wrote:

Strange bug in the combine.py utility for combining kallisto output from multiple samples:

Works well for combining the data from multiple samples but the target_id's in the first column of the output count.txt file are not in the same order as the target_id's in the input abundance.tsv file.

Any suggestions would be appreciated.

Example:

count.txt

target_id length eff_length X1 X2 X3 BW12616D

chr12:121593689-121594058 369 170 0 0 0 188

chr1:84971674-84972004 330 131 0 0 0 101

chr1:226074938-226075857 919 720 0 0 0 688.182

chr10:55664705-55665110 405 206 0 0 0 148 ...

abundance.tsv

target_id length eff_length est_counts tpm

chr1:65409-65725 316 117 75.5 4.08319

chr1:65731-66073 342 143 104.764 4.6357

chr1:69381-69700 319 120 322.447 17.0026 ...

software error • 97 views
ADD COMMENTlink modified 23 days ago • written 24 days ago by rvjensen10
1
gravatar for rvjensen
23 days ago by
rvjensen10
rvjensen10 wrote:

No bug! It turns out that dictionaries are correctly ordered in Python 3. I was using the old Mac default Python 2.7 (Sigh).

ADD COMMENTlink written 23 days ago by rvjensen10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1653 users visited in the last hour
_