Finding non ATGCN sequences using python
I want to search for non sequences in a file, to substitute them.
for file in SeqIO.parse('myfile.fasta', 'r'):
sequence = re.search('^ATGCN', record.seq)
output_file=open('files.fasta, 'w+')
• 917 views
•
link
updated
by
Istvan Albert
•
written
by
AveryB •
•
written
by
AveryB •
1 answer
Go to a site like:
to test your regular expressions (select Python as the language). To match characters except ATGCN the pattern would be:
[^ACGTN]
you get quite helpful tips in the sidebar
• 0 views
•
link
written
by
Istvan Albert
Log in to answer this question.
More posts like this
-
UniprotKB BLAST pipeline troubles
written by Aamir Mehmood •
-
When running this code on a fasta file and a max() arg error occurs
written by rt2421 •
-
A3 count (count of A occurring in 3rd codon position) in multifasta file
written by annushreekurmi •
-
Entrez search to CSV file
written by pmramos95 •
-
remove sequences with non-canonical nucleotides from fasta file
written by grant.hovhannisyan
-
FASTA File Parsing with many species
written by adeel.maliks20 •
-
Fasta Sliding Window
written by jordanmwilliams.94 •
-
How can I handle UnknownSeq?
written by fiascox12 •
-
Filter sequences by length and rename them with biopython
written by linlinli99013 •
-
Extracting Fasta Sequence Using Biopython
written by Ma
