我是Perl编程的初学者。我现在正在研究的问题是如何从文本文件中获取基因长度。文本文件包含基因名称(第10列),起始位点(第6列),结束位点(第7列)。长度可以从第6列和第7列的差异中得出。但是我的问题是如何将基因名称(来自第10列)与从第6列和第7列的差异中得出的相应差异进行匹配。非常感谢! / p>
open (IN, "Alu.txt");
open (OUT, ">Alu_lengthsubfam.csv");
while ($a = <IN>){
@data = split (/\t/, $a);
$genelength = $data[7] - $data[6]++;
$subfam = $data[10]
}
foreach $subfam (0...$#subfam){
print OUT "$subfam, $genelength{$subfam}\n";
}
close (IN);
close (OUT)
答案 0 :(得分:1)
perl -F'\t' -lane'print $F[10], "\t", $F[7]-$F[6]' \
< Alu.txt > Alu_lengthsubfam.csv