Differences

This shows you the differences between two versions of the page.

--- metamos_test_2017 [2017/03/23 09:59] – [run_pipeline_test.sh의 실패 원인은?] hyjeong
+++ metamos_test_2017 [2021/03/17 13:09] (current) – external edit 127.0.0.1
@@ Line 4: / Line 4: @@
 ===== Candidatus Carsonella ruddii 샘플을 이용한 테스트 =====
-MetAMOS의 샘플 데이터는 160 kb의 작은 유전체를 지닌 bacterial endosymbiont "Candidatus Carsonella ruddii PV"의 50x read를 사용하는 것으로 되어있다. fasta 및 fastq 파일과 reference(complete)가 전부 존재한다.
+MetAMOS의 샘플 데이터는 160 kb의 작은 유전체를 지닌 bacterial endosymbiont "Candidatus Carsonella ruddii PV"의 50x read를 사용하는 것으로 되어있다. fasta 및 fastq 파일과 reference(complete)가 전부 존재한다. _pe_가 암시하듯이 interleaved 형태의 paired end sequence이다.
   8000000 bp / 100000 seqs; 80.0 average length -- carsonella_pe_filt.fna
@@ Line 17: / Line 17: @@
 ==== runPipeline의 이해 ====
 -f STRING(force this step)과 -n STRING(step to skip)을 지정하지 않으면 MultiAlign, FunctionalAnnotation, FindRepeats, FindScaffoldORFS 단계를 생략한다. 실제로 사용된 명령어 라인은 output_directory/Logs/COMMANDS.log를 열어보면 그대로 나온다.
+=== runPipeline -d test1만 수행하면? ===
+default assembler인 SOAPdenovo(k-mer 31)에 의한 조립은 되지만 Abundance(no plot), Annotate(100% no hits) Classify (all unknown: 13 = 17 contigs + 6 scaffolds) 결과는 이상하다. Postprocess에 의하여 proba.ctg.fa와 proba.scf.fa가 나온다. 이제 옵션을 하나씩 추가해 보자.
+=== runPipeline -d test1 -k 55 -f Abundance,Annotate,Classify ===
+Abundance(no plot), Annotate(100% no hits) Classify (all unknown)
+=== runPipeline -d test1 -k 55 -c fcp -f Abundance,Annotate,Classify ===
+default annotate program인 FCP를 명시적으로 설정한다. 드디어 Annotate 결과가 나타났고 Classify도 100% Gammaproteobacteria로 나왔다. 그러나 Abundance는 여전히 플롯을 그리지 못한다.
+=== runPipeline -d test1 -k 55 -c kraken -f Abundance,Annotate,Classify ===
+-c kraken을 설정하면 Annotate(100% No hits), Classify(unknown) 그러나 Abundance는 플롯이 없다.
+=== runPipeline -c kraken -d test1 -k 55 -f Assemble,MapReads,FindORFS,Annotate,Propagate,Classify,Abundance ===
+이렇게 하면 run_pipeline_test.sh와 거의 같아진다. 그러나 결과는 나아지지 않았다. **도저히 이해하기 어렵지만 run_pipeline_test.pl의 마지막 라인과 같은 형식으로 하지 않으면 잘 되지 않는다는 뜻이다.**
 ==== run_pipeline_test.sh의 실패 원인은? ====