Skip to content

Commit

Permalink
update to fit Seqstat.py
Browse files Browse the repository at this point in the history
  • Loading branch information
FlyPythons committed Mar 16, 2018
1 parent 9caf3b4 commit ebd2ddf
Show file tree
Hide file tree
Showing 2 changed files with 13 additions and 16 deletions.
25 changes: 11 additions & 14 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -18,30 +18,27 @@ if the error is raised, install is ok.
## 2. Example
* Get sequences distribution in FASTA/Q format:

for one FASTA/Q file
```commandline
fastaStat.py in.fa > in.fa.stat
fastqStat.py in.fq > in.fq.stat
```
for multi FASTA/Q files
for FASTA/Q files, different formatted files can be stat together

```commandline
fastaStat.py 1.fa 2.fa *.fa > in.fa.stat
fastqStat.py 1.fq 2.fq *.fq > in.fq.stat
Seqstat.py 1.fq *.fa > in.stat
```

use '-c' to speed up the process
```commandline
fastaStat.py -c 10 1.fa 2.fa *.fa > in.fa.stat
fastqStat.py -c 10 1.fq 2.fq *.fq > in.fq.stat
Seqstat.py -c 10 1.fq *.fa > in.stat
```

for one file contains FASTA/Q file paths
```commandline
fastaStat.py -f in.fofn > in.fa.stat
fastqStat.py -f in.fofn > in.fq.stat
Seqstat.py -f -c 10 in.fofn > in.stat
```
for *NGS short reads*, you'd better turn on *'-ngs'* to avoid out of memory

for *NGS short reads*, you'd better turn on *'-ngs'* to avoid meaningless stat
```commandline
fastqStat.py -ngs -c 10 *.R1.fq *.R2.fq
Seqstat.py -ngs -c 10 *.R1.fq *.R2.fq
```

* Split large sequences file in FASTA/Q format into small files:

split by sequences number == {max number} per file
Expand Down
4 changes: 2 additions & 2 deletions test/test.py
Original file line number Diff line number Diff line change
Expand Up @@ -16,7 +16,7 @@ def test_fasta():
with open("test.fasta", "w") as out:
out.write(fasta_record)

os.system("python ../fastaStat.py test.fasta")
os.system("python ../SeqStat.py test.fasta")
os.system("python ../fastaSplit.py -m length -n 1 -o fasta test.fasta")


Expand All @@ -40,7 +40,7 @@ def test_fastq():
with open("test.fastq", "w") as out:
out.write(fastq_record)

os.system("python ../fastqStat.py test.fastq")
os.system("python ../Seqstat.py test.fastq")
os.system("python ../fastqSplit.py -m length -n 1 -o fastq test.fastq")


Expand Down

0 comments on commit ebd2ddf

Please sign in to comment.