保存一个文本字符串，如果它由另一个特定的文本字符串进行？
回答此问题可获得 20 贡献值，回答如果被采纳可获得 50 分。
<p>很抱歉标题不好，我不知道该怎么回答我的问题</p>
<p>我写了一个脚本，从fastq文件（纯文本基因组读取文件）中提取数据。每第一行是一个标题，第二行是一个基本字符串-第三和第四行是不需要的</p>
<pre><code>filename = 'C0_GGCTAC_R1_no_adapter_trimming.fastq'
new_filename = filename[:-9] + '_new.fastq'
with open(filename) as f_obj:
file_contents = f_obj.readlines()
extracted_lines = ''
line_count = 0
# Pull header and base lines
for line in file_contents:
line_count += 1
# Headers
if line_count == 1:
extracted_lines += line
# Reads ending in A
elif line_count == 2 and line[-2] == 'A':
extracted_lines += line
# Reads ending in G
elif line_count == 2 and line[-2] == 'G':
extracted_lines += line
# Reset counter
elif line_count == 4:
line_count = 0
with open(new_filename, 'w') as f_obj:
f_obj.write(extracted_lines)
print(new_filename + " was created.")
</code></pre>
<p>只要基的读取以A或G结束，脚本就会提取每次读取的头和读取中的基字符串。
输入文件的示例如下：</p>
<pre><code>@HWI-D00461:137:C9H2FACXX:3:1101:1239:1968 1:N:0:GGCTAC
NTGTGTAATAGATTTTACTTTTGCCTTTAAGCCCAAGGTCCTGGACTTGAAACATCCAAGGGATGGAAAATGCCGTATAACAGGGTGGAAGAGAGATTTGA
+
#1=BDDFFHHHFHIJJJJJJJJJJJJJJJJJJJJJIJJIJJJJJHJIIJHGIJJJJJJIHJJBGHJHIIJJJHHHHFFFFEEEDD;?BACDDDA?@CDDDC
@HWI-D00461:137:C9H2FACXX:3:1101:1117:1968 1:N:0:GGCTAC
NAAAGTCTACCAATTATACTTAGTGTGAAGAGGTGGGAGTTAAATATGACTTCCATTAATAGTTTCATTGTTTGGAAAACAGAGGTAATTTTTGATACAGA
+
#1=DDDFDFHHHGHIIGJJJJHIJIHHDIHHIJGGEI@GFGHIHIJHEFHIIIIGIJGHHGECFGIDHGIHIIEGIIJHHEEFFF7?ACEECCBBDEDDDC
</code></pre>
<p>输出文件如下所示</p>
<pre><code>@HWI-D00461:137:C9H2FACXX:3:1101:1117:1968 1:N:0:GGCTAC
NAAAGTCTACCAATTATACTTAGTGTGAAGAGGTGGGAGTTAAATATGACTTCCATTAATAGTTTCATTGTTTGGAAAACAGAGGTAATTTTTGATACAGA
@HWI-D00461:137:C9H2FACXX:3:1101:1200:1972 1:N:0:GGCTAC
@HWI-D00461:137:C9H2FACXX:3:1101:1087:1973 1:N:0:GGCTAC
NTAATCCAACTAACTAAAAATAAAAAGATTCAAATAGGTACAGAAAACAATGAAGGTGTAGAGGTGAGAAATCAACAGGATGTTCAGAAGCCTGTGTATGA
</code></pre>
<p>尽管它包含了所需的所有数据，但它会拉出每个标题行（以“@”开头），这是不必要的</p>
<p>如果代码是由以a或G结尾的一串基进行的，如何修改代码以只拉出标题行</p>
请先 登录 后评论
针对此问题: 更多的回答
- 89 关注
- 1 收藏，216 浏览
- 网友 提问于 2天前
相关Python问题
- 如何添加虚拟方法 5 回答
- 如何添加表示整数的擦边字符串？ 10 回答
- 如何添加要在Bokeh中使用的新font.ttf文件？ 3 回答
- 如何添加要显示的矩阵XY轴编号和XY轴 4 回答
- 如何添加计数？ 10 回答
- 如何添加计数器函数？ 1 回答
- 如何添加计数器列来计算数据帧中另一列中的特定值？ 4 回答
- 如何添加计数器来跟踪while循环中的月份和年份？ 8 回答
- 如何添加计数并删除countplot的顶部和右侧脊椎？ 3 回答
- 如何添加计时器wx.应用程序更新窗口对象的主循环？ 2 回答
- 如何添加评论到帖子？PostDetailVew，Django 2.1.5 8 回答
- 如何添加评论拉梅尔亚姆 4 回答
- 如何添加诸如矩阵Python/Pandas之类的数据帧？ 2 回答
- 如何添加谷歌地点自动完成到Flask？ 3 回答
- 如何添加超时、python discord bot 10 回答
- 如何添加超过1dp的检查 1 回答
- 如何添加距离方法 7 回答
- 如何添加跟随游戏的敌人精灵 2 回答
- 如何添加路径以便python可以找到程序？ 5 回答
- 如何添加身份验证/安全性以使用happybase访问HBase？ 3 回答
