Comments (1)
Hello,
I would like to use -run_mode 3, but i'm not clear about which files to use. You provide this as your example:
python3 MitoZ.py all2 --genetic_code 5 --clade Arthropoda --outprefix test
--thread_number 12 --fastq1 clean.1.fq.gz --fastq2 clean.2.fq.gz
--fastq_read_length 150 --insert_size 250
--run_mode 3
--filter_taxa_method 1
--requiring_taxa 'Arthropoda'
--quick_mode_seq_file quickMode.fa
--quick_mode_fa_genes_file quick_mode_fa_genes.txt
--missing_PCGs ND4L ND6 ND2
--quick_mode_score_file work71.hmmtblout.besthit.sim.filtered.high_abundance_10.0X.reformat.sorted
--quick_mode_prior_seq_file work71.hmmtblout.besthit.sim.filtered.faBut I don't know how to write the command with these output files:
IC11_AAH1795_Hymenoptera_Ichneumonidae_Campodorus_ultimus.cds
IC11_AAH1795_Hymenoptera_Ichneumonidae_Campodorus_ultimus.fasta
IC11_AAH1795_Hymenoptera_Ichneumonidae_Campodorus_ultimus.misc_feature
IC11_AAH1795_Hymenoptera_Ichneumonidae_Campodorus_ultimus.rrna
IC11_AAH1795_Hymenoptera_Ichneumonidae_Campodorus_ultimus.trna
IC11_AAH1795_Hymenoptera_Ichneumonidae_Campodorus_ultimus_mitoscaf.fa.gbf
IC11_AAH1795_Hymenoptera_Ichneumonidae_Campodorus_ultimus_mitoscaf.fa.sqn
IC11_AAH1795_Hymenoptera_Ichneumonidae_Campodorus_ultimus_mitoscaf.fa.tbl
IC11_AAH1795_Hymenoptera_Ichneumonidae_Campodorus_ultimus_mitoscaf.fa.val
README.txt
circos.png
circos.svg
errorsummary.val
summary.txt
work71.hmmtblout.besthit.sim.filtered.high_abundance_10.0X.reformat.sorted.Not-picked
work71.hmmtblout.besthit.sim.filtered.high_abundance_10.0X.reformat.sorted.Not-picked.fa
work71.hmmtblout.besthit.sim.filtered.low_abundance
work71.hmmtblout.besthit.sim.filtered.low_abundance.fasta
work71.mitogenome.fa
work71.most_related_species.txt
- which mitogenome fasta file to use?
- which work71.hmmtblout.besthit.sim.filtered files to use?
- there are no clean.{1,2}.fq.gz files. Should I use
python3 MitoZ.py all
instead?This is my attempt:
python3 MitoZ.py all2 --genetic_code 5 --clade Arthropoda --outprefix test
--thread_number 12 --fastq1 clean.1.fq.gz --fastq2 clean.2.fq.gz
--fastq_read_length 150 --insert_size 250
--run_mode 3
--filter_taxa_method 1
--requiring_taxa 'Arthropoda'
--quick_mode_seq_file work71.mitogenome.fa
--quick_mode_fa_genes_file work71.mitogenome_genes.txt
--missing_PCGs ATP6 ATP8 COX3 ND2 ND3 ND4 ND4L ND5
--quick_mode_score_file work71.hmmtblout.besthit.sim.filtered.high_abundance_10.0X.reformat.sorted.Not-picked
--quick_mode_prior_seq_file work71.hmmtblout.besthit.sim.filtered.high_abundance_10.0X.reformat.sorted.Not-picked.faWhere
work71.mitogenome_genes.txt
contains the following text:scaffold41 ND1 CYTB ND6
C5251 COX2 COX1This is from the summary.txt
#Seq_id Length(bp) Circularity Closely_related_species
scaffold41 2499 no Diadegma semiclausum
C5251 2608 no Diadegma semiclausum#Seq_id Start End Length(bp) Direction Type Gene_name Gene_prodcut Total_freq_occurred
scaffold41 21 930 910 + CDS ND1 NADH dehydrogenase subunit 1 1
scaffold41 933 1000 68 - tRNA trnS(uga) tRNA-Ser 1
scaffold41 1032 2191 1160 - CDS CYTB cytochrome b 1
scaffold41 2183 2486 304 - CDS ND6 NADH dehydrogenase subunit 6 1
C5251 92 156 65 - tRNA trnD(guc) tRNA-Asp 1
C5251 155 227 73 - tRNA trnK(cuu) tRNA-Lys 1
C5251 227 903 677 - CDS COX2 cytochrome c oxidase subunit II 1
C5251 962 2501 1540 - CDS COX1 cytochrome c oxidase subunit I 1Protein coding genes totally found: 5
tRNA genes totally found: 3
rRNA genes totally found: 0
Genes totally found: 8Potential missing genes:
#Gene total_missing_number
ATP6 1
ATP8 1
COX3 1
ND2 1
ND3 1
ND4 1
ND4L 1
ND5 1
l-rRNA 1
s-rRNA 1
tRNA-Ala 1
tRNA-Arg 1
tRNA-Asn 1
tRNA-Cys 1
tRNA-Gln 1
tRNA-Glu 1
tRNA-Gly 1
tRNA-His 1
tRNA-Ile 1
tRNA-Leu 2
tRNA-Met 1
tRNA-Phe 1
tRNA-Pro 1
tRNA-Ser 1
tRNA-Thr 1
tRNA-Trp 1
tRNA-Tyr 1
tRNA-Val 1
The missing genes might be foud from the
'.high_abundance' and '.low_abundance' files!
Hi Doug,
Sorry for the confusion.
In short, you can look for the files required by -run_mode 3
in your quick-mode-run tmp
directory.
which mitogenome fasta file to use?
you can use the the file work71.mitogenome.fa
. this is a file that should only contain correct mitochondrial sequences.
This file is also in tmp/ZZZ.assembly
which work71.hmmtblout.besthit.sim.filtered files to use?
In the tmp/ZZZ.assembly
directory.
there are no clean.{1,2}.fq.gz files. Should I use python3 MitoZ.py all instead?
Generally, before NGS assembly, we have to filter our raw reads (provided by sequencing companies) and then use the clean reads to perform assembly.
-
Suppose now you have the raw reads, then you can run MitoZ with the
all
subcommand and with quick mode. In this case, MitoZ has output the clean reads and assembly under thetmp/ZZZ.cleandata
directory, -
If you already have your own 'clean' fastq files, you can just provide them to MitoZ for any run mode.
When you finish, you may also found some good results in work31.mitogenome.fa
and work91.mitogenome.fa
(in the result directory and/or tmp/ZZZ.assembly2
) if you are lucky.
Please just ask if you still get into troubles.
guanliang
from mitoz.
Related Issues (20)
- MitoZ annotate error? HOT 1
- Error with filter_by_abundance() HOT 3
- Arguments in findmitoscaf HOT 3
- Use for plant mitochondria annotation HOT 4
- Nanopore sequence for Mitochondrial genome assembly? HOT 4
- Is longer Illumina reads or shorter reads preferred? HOT 3
- mitoz visualize -h HOT 3
- 2 Mitogenomes in One Sample HOT 2
- Other animal taxa? HOT 3
- Error encountered in mt_assembly using mitoz all command HOT 2
- How to introduce breakpoints? HOT 1
- overlap_information file empty and mitogenomes not circularised HOT 1
- How to make cirular map to linear? HOT 1
- I got an error when I install the taxdump.tar.gz into sqlite3, why ? HOT 2
- Open Reading Frame HOT 1
- Missing genes with Udocker HOT 4
- mitoz3.4 annotations draw circle diagrams all have a gap, but in other annotation platforms are circular HOT 13
- Annotation of Gene on Both Sides of Start and End Position Got Cut Off HOT 1
- Can not find assembly_file: HOT 13
- Incorrect stop codon HOT 8
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from mitoz.