Hello,
When I am trying to do samtools
on each of sequence, I met the problem of seperating them.
I have file gene.info
which contain contents:
ENA|MBYW01000151|MBYW01000151.1 34978 35427
ENA|MBYW01000151|MBYW01000151.1 36152 36476
what I want to have is them separate,
file1
with content ENA|MBYW01000151|MBYW01000151.1 34978 35427
file2
with contnet ENA|MBYW01000151|MBYW01000151.1 36152 36476
I thought it will be easy to achieve with cut
command, but
cat geneinfo.txt | cut -f 1
ENA|MBYW01000151|MBYW01000151.1
ENA|MBYW01000151|MBYW01000151.1
cat geneinfo.txt | cut -f 2
34978
36152
cat geneinfo.txt | cut -f 3
35427
36476
which is really wield, I don't know why?
Thanks for your helping.
for more information, geneinfo.txt
is obtained by run
blastn -query S288C_YLL052C_AQY2_genomic.fsa -db beer001.fa -outfmt "6 sseqid send sstart" >> geneinfo.txt
head -n 1 gene.info > file1
head -n 2 gene.info > file2
What does that mean?
samtool is a tool to clip certain regions of genes.