miwipe / ngsLCA

GNU General Public License v3.0
9 stars 5 forks source link

problems running the test data #30

Closed slennart closed 1 year ago

slennart commented 1 year ago

I've been having problems running ngsLCA on my own data so I tried running the test data first, but I also can't get that to work. I downloaded SPL_015_1444.fq.plastids.bam from ERDA, downloaded the ncbi tax_tax_dmp and then ran the following code:

ngsLCA -editdistmin 0 -editdistmax 0 -names ncbi_tax_dmp/names.dmp -nodes ncbi_tax_dmp/nodes.dmp -acc2tax ncbi_tax_dmp/nucl_gb.accession2taxid.gz -bam SPL_015_1444.fq.plastids.bam -outnames outfile.ed0

slennart commented 1 year ago

ngsLCA_error2

miwipe commented 1 year ago

Hi slennart,

Running the latest version of the ngsLCA with the command you provided and downloading the file again as well as the ncbi taxonomy, gives me no errors. Could you check you have downloaded the file correctly, have the recent version of ngsLCA and ncbi taxonomy. Another thing is that the bam/sam files have to be sorted by coordinate e.g. readID.

Mikkel

Screenshot 2023-01-05 at 11 35 30