如何用Python实现去除特定长度核酸接头序列的技巧？

本文共计122个文字，预计阅读时间需要1分钟。

输入：input.txtATTCGATTATAAGCTCGATCGATCGATCGATCGATCGATCGATCGATCGATCATTCGATTATAAGCACTGATCGATCGATCGATCGATGCTATCGTCGTATTCGATTATAAGCATCGATCACGATCTATCGTACGTATGCATATCGATATCGATCGTAGTCATTCGATTATAAGCACTATCGATGATCTAGCTACGATCGTAGCTGTAATT

输入

input.txt

ATTCGATTATAAGCTCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATC
ATTCGATTATAAGCACTGATCGATCGATCGATCGATCGATGCTATCGTCGT
ATTCGATTATAAGCATCGATCACGATCTATCGTACGTATGCATATCGATATCGATCGTAGTC
ATTCGATTATAAGCACTATCGATGATCTAGCTACGATCGTAGCTGTA
ATTCGATTATAAGCACTAGCTAGTCTCGATGCATGATCAGCTTAGCTGATGATGCTATGCA

编码

file=open("input.txt")
output=open("trimmed.txt","w")

fordnainfile:
trimmed_dna=dna[14:]#接头特定长度14
trimmed_length=len(trimmed_dna)-1
output.write(trimmed_dna)

print("processedsequencewithlength"+str(trimmed_length))

输出

trimmed.txt

TCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATC
ACTGATCGATCGATCGATCGATCGATGCTATCGTCGT
ATCGATCACGATCTATCGTACGTATGCATATCGATATCGATCGTAGTC
ACTATCGATGATCTAGCTACGATCGTAGCTGTA
ACTAGCTAGTCTCGATGCATGATCAGCTTAGCTGATGATGCTATGCA

作者：Bioinfarmer

若要及时了解动态信息，请关注同名微信公众号：Bioinfarmer。