A simple bioinformatics package
Project description
pyoinformatics ๐
Examples
Find the reverse complement of all the sequences in a file:
with open('out.fasta', 'w') as f:
for seq in Bio.read_fasta('in.fasta'):
f.writelines(seq.reverse_complement().to_fasta())
Count the number of occurrences of 'ATG' in seq1
seq1.count('ATG')
Count the number of occurrences of 'ATG' in seq1 that differ by <= 1 base.
seq1.count('ATG', 1)
Find the number of occurrences of 'ATG' or 'AAG' in seq1
len(seq1.find('A[AT]G'))
Find the average position of all occurrences of 'ATG' in a fasta file
from statistics import mean
for seq in Bio.read_fasta('in.fasta'):
print(mean(seq.find('ATG')))
ASCI plot the relative nt counts for all the sequences in a file
for seq in Bio.read_fasta('in.fasta'):
counts = seq.counts
print(f">{seq.id}")
for nt in sorted(counts.keys()):
bar = int((counts[nt]/len(seq))*100)
print(f"{nt}: {'โ' * bar}")
>HSBGPG Human gene for bone gla protein (BGP)
A: โโโโโโโโโโโโโโโโโ
C: โโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโ
G: โโโโโโโโโโโโโโโโโโโโโโโโโโโโโโ
T: โโโโโโโโโโโโโโโโโโ
>HSGLTH1 Human theta 1-globin gene
A: โโโโโโโโโโโโโโ
C: โโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโ
G: โโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโโ
T: โโโโโโโโโโโโโโโ
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
pyoinformatics-0.0.2.tar.gz
(4.9 kB
view hashes)
Built Distribution
Close
Hashes for pyoinformatics-0.0.2-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 6132c9479b1b07e7c9a3423043e00e8cbbf94472d8da6aaff8a4c45ec5cfbd12 |
|
MD5 | cc22f524ac2a0c2d8bce1aeb85794634 |
|
BLAKE2b-256 | 23ede2fc0f880dd9350d6b2d91b6f068b8ead74ac2c8549241b624e7d072cd9e |