-
Notifications
You must be signed in to change notification settings - Fork 18
Expected output log for the test example
Sina Majidian edited this page Sep 14, 2022
·
2 revisions
This is the expected stdout log for the test example.
read2tree --tree --standalone_path marker_genes/ --reads sample_1.fastq sample_2.fastq --output_path output/
--- Load OGs with min 0 species from oma marker gene export! ---
Loading files for pre-filter: 100%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 20/20 [00:00<00:00, 241.13 OGs/s]
2022-09-14 12:11:37,788 - read2tree.OGSet - INFO - --- Load ogs and find their corresponding DNA seq using the REST api ---
Loading OGs: 100%|█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 20/20 [00:04<00:00, 4.49 OGs/s]
2022-09-14 12:11:42,250 - read2tree.OGSet - INFO - sample_1: Gathering of DNA seq for 20 OGs took 4.461216449737549.
--- Generating reference for mapping ---
Loading records: 100%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 20/20 [00:00<00:00, 56833.39 record/s]
2022-09-14 12:11:42,251 - read2tree.ReferenceSet - INFO - sample_1: Extracted 5 reference species form 20 ogs took 0.0009484291076660156
--- Alignment of 20 OGs ---
2022-09-14 12:11:45,821 - read2tree.Aligner - INFO - sample_1: Alignment of 20 OGs took 3.566136598587036.
--- Mapping of reads to reference sequences ---
Mapping reads to species: 0%| | 0/5 [00:00<?, ? species/s]2022-09-14 12:11:45,830 - read2tree.Mapper - INFO - sample_1: --- Mapping of reads to HUMAN reference species ---
[E::idx_find_and_load] Could not retrieve index file for '/tmp/ngm_j0p1u43c/HUMAN_OGs.fa.bam'
2022-09-14 12:11:49,884 - read2tree.Mapper - INFO - sample_1: Mapped 29 / 50000 reads to HUMAN_OGs.fa
2022-09-14 12:11:49,913 - read2tree.Mapper - INFO - sample_1: Mapping to HUMAN_OGs.fa references took 4.08157205581665.
Mapping reads to species: 20%|███████████████████████▌ | 1/5 [00:04<00:17, 4.33s/ species]2022-09-14 12:11:50,159 - read2tree.Mapper - INFO - sample_1: --- Mapping of reads to GORGO reference species ---
[E::idx_find_and_load] Could not retrieve index file for '/tmp/ngm_j0p1u43c/GORGO_OGs.fa.bam'
2022-09-14 12:11:52,815 - read2tree.Mapper - INFO - sample_1: Mapped 27 / 50000 reads to GORGO_OGs.fa
2022-09-14 12:11:52,831 - read2tree.Mapper - INFO - sample_1: Mapping to GORGO_OGs.fa references took 2.6480491161346436.
Mapping reads to species: 40%|███████████████████████████████████████████████▏ | 2/5 [00:07<00:10, 3.38s/ species]2022-09-14 12:11:52,880 - read2tree.Mapper - INFO - sample_1: --- Mapping of reads to RATNO reference species ---
[E::idx_find_and_load] Could not retrieve index file for '/tmp/ngm_j0p1u43c/RATNO_OGs.fa.bam'
2022-09-14 12:11:56,418 - read2tree.Mapper - INFO - sample_1: Mapped 44 / 50000 reads to RATNO_OGs.fa
2022-09-14 12:11:56,447 - read2tree.Mapper - INFO - sample_1: Mapping to RATNO_OGs.fa references took 3.5645296573638916.
Mapping reads to species: 60%|██████████████████████████████████████████████████████████████████████▊ | 3/5 [00:10<00:06, 3.49s/ species]2022-09-14 12:11:56,496 - read2tree.Mapper - INFO - sample_1: --- Mapping of reads to XENLA reference species ---
[E::idx_find_and_load] Could not retrieve index file for '/tmp/ngm_j0p1u43c/XENLA_OGs.fa.bam'
2022-09-14 12:11:59,318 - read2tree.Mapper - INFO - sample_1: Mapped 19 / 50000 reads to XENLA_OGs.fa
2022-09-14 12:11:59,334 - read2tree.Mapper - INFO - sample_1: Mapping to XENLA_OGs.fa references took 2.8225371837615967.
Mapping reads to species: 80%|██████████████████████████████████████████████████████████████████████████████████████████████▍ | 4/5 [00:13<00:03, 3.25s/ species]2022-09-14 12:11:59,371 - read2tree.Mapper - INFO - sample_1: --- Mapping of reads to MNELE reference species ---
[E::idx_find_and_load] Could not retrieve index file for '/tmp/ngm_j0p1u43c/MNELE_OGs.fa.bam'
2022-09-14 12:12:02,773 - read2tree.Mapper - INFO - sample_1: Mapped 3 / 50000 reads to MNELE_OGs.fa
2022-09-14 12:12:02,801 - read2tree.Mapper - INFO - sample_1: Mapping to MNELE_OGs.fa references took 3.4108190536499023.
Mapping reads to species: 100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 5/5 [00:17<00:00, 3.40s/ species]
2022-09-14 12:12:02,834 - read2tree.Mapper - INFO - sample_1: Mapping to all references took 17.006091117858887.
Adding mapped seq to alignments: 100%|██████████████████████████████████████████████████████████████████████████████████████████████████████| 20/20 [00:00<00:00, 570653.61 alignments/s]
Adding mapped seq to OG: 100%|█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 20/20 [00:00<00:00, 590747.04 OGs/s]
--- Add inferred mapped sequence back to OGs ---
Adding mapped seq to OG: 100%|█████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 20/20 [00:00<00:00, 122640.47 OGs/s]
2022-09-14 12:12:02,845 - read2tree.OGSet - INFO - sample_1: Appending 7 reconstructed sequences to present OG took 0.0011239051818847656.
--- Add inferred mapped sequence back to alignment ---
Adding mapped seq to alignments: 100%|███████████████████████████████████████████████████████████████████████████████████████████████████████| 20/20 [00:00<00:00, 19840.61 alignments/s]
2022-09-14 12:12:02,860 - read2tree.Aligner - INFO - sample_1: Appending 1 reconstructed sequences to present Alignments took 0.0012431144714355469.
--- Tree inference ---
2022-09-14 12:12:10,413 - read2tree.TreeInference - INFO - sample_1: Tree inference took 7.521272659301758.
(sample_1:0.0057385847,((HUMAN:0.0041201107,GORGO:0.0272788512):0.0432944811,(XENLA:0.1714987394,MNELE:0.9177713687):0.1141522197):0.0628937266,RATNO:0.0107653458);