Saskia-Oosterbroek / decona

fastq to polished sequenses: pipeline suitable for mixed samples and long (Nanopore) reads
MIT License
43 stars 12 forks source link

Issue with Medaka #44

Open lewhiteside opened 2 years ago

lewhiteside commented 2 years ago

Hi Saskia, I have an issue that Im not sure how to solve.

After installing decona I ran the following code:

decona -d -l 300 -m 2000 -q 10 -c 0.80 -n 100 -M

The output file was as follows (error can be found at the bottom):

Filtering data...
Data filtered with NanoFilt
total raw sequences = 303896
total filtered sequences = 224966
Demultiplexing...
Multi threading is not yet supported in epi2mme mode. Falling back to using a single thread.
Adapters detected in 3288 of 224966 reads
  PBK004/LWB001   2552: |                      |   1.13 %
  RAB204/RAB214    736: |                      |   0.33 %
           none 191852: |    ################# |  85.28 %
Barcodes detected in 3288 of 224966 adapters
           BC01    124: |                      |   0.06 %
           BC02    163: |                      |   0.07 %
           BC03    103: |                      |   0.05 %
           BC04    158: |                      |   0.07 %
           BC05    186: |                      |   0.08 %
           BC11      2: |                      |   0.00 %
      barcode01     35: |                      |   0.02 %
      barcode02    133: |                      |   0.06 %
      barcode03   2251: |                      |   1.00 %
      barcode04     70: |                      |   0.03 %
      barcode05     63: |                      |   0.03 %
           none 191852: |    ################# |  85.28 %
29826 reads were skipped due to the min. length filter.
Demultiplexing finished in 260.19s
Data demultiplexed, working directory changed to: 
/mnt/scratch2/users/40266190/Decona/decona/demultiplexed_data
Fastq reads are being transformed to fasta
cat: barcode01: Is a directory
cat: barcode02: Is a directory
cat: barcode03: Is a directory
cat: barcode04: Is a directory
cat: barcode05: Is a directory
Transforming fastq to fasta Complete
Clustering reads...
Clustering barcode01/...
Clustering barcode02/...
Clustering barcode03/...
Clustering barcode04/...
Clustering barcode05/...
Clustering result/...
Clustering complete.
Done
Done
cat: consensus_medaka_159-60.fa: Is a directory
cat: consensus_medaka_191-37.fa: Is a directory
cat: consensus_medaka_239-64.fa: Is a directory
Aligning and making draft assembly of 0-consensus_medaka_159-60.fa.fa...
tail: error reading ‘0-consensus_medaka_159-60.fa.fa’: Is a directory

The version of decona I am running in 0.1.3,

Any help with this would be really appreciated,

Thank you.