เรียกใช้คำสั่งบนหลาย ๆ ไฟล์และบันทึกผลลัพธ์สำหรับแต่ละไฟล์ในไฟล์ที่แตกต่างกันอย่างไร [ปิด]


-1

ฉันมี vcf-file (file.vcf) และไฟล์ประชากรหลายไฟล์(HA1.txt, HA1.txt, HA2.txt,...,HA28.txt)ที่มี ID แต่ละรายการ

ไฟล์ประชากรแต่ละไฟล์มีลักษณะดังนี้ (แต่ละบรรทัดเป็นหนึ่งไฟล์):

$ cat HA1.txt
    QQ48_SD1A-37    
    ED19_SD1A40-3_357
    TT335_SD1A-20
    HH356_SD1A-7
    Q029B_SD1A38    
    HT73_SD1A-28    
    HT288_SD1A-24
    Q004B_SD1A-1    
    Q027_SD1A-4
    Q096_SD1A-40

ฉันต้องการแยกความลึกของการครอบคลุมสำหรับแต่ละบุคคลในแต่ละประชากร ฉันจะรันคำสั่งนี้ได้อย่างไร: grep -v "^#" file.vcf | cut -f 10 | cut -d ':' -f2สำหรับแต่ละคนในไฟล์ประชากรและบันทึกผลลัพธ์ของประชากรแต่ละคนในไฟล์แยกกัน

ผลลัพธ์ที่ฉันต้องการสำหรับประชากรแต่ละคนมีลักษณะดังนี้:

.   6   4   6   .   5   .   10  .   7   .
.   9   16  8   3   8   .   16  9   22  .
.   8   11  8   .   8   .   16  8   18  .
6   20  12  20  12  28  3   24  4   26  14
6   25  15  24  13  32  3   25  3   25  15

แต่ละคอลัมน์ในหนึ่งบุคคล!


ใช่ฉันเห็นด้วยและพยายามถามคำถามของฉันในทางที่ดีขึ้น!
Anna1364

โปรดแสดงเวอร์ชันที่เป็นตัวอย่างของคุณfile.vcfโปรดที่นี่ฉันไม่สามารถจำแนกได้ว่าผลลัพธ์ที่คุณคาดหวังจะสร้างด้วยอินพุตตัวอย่าง HA1.txtซึ่งไม่เกี่ยวข้องกับเอาต์พุตและคำสั่งที่คุณใช้อยู่
αғsнιη

Chrom1 1308 A 6622.86 T 0448 GT: DP: DPR: RO: QR: AO: QA: GL 0/0: 8: 8,0: 8: 292: 0: 0: 0, -2.40824, -26.6303 0/0: 6: 6,0 : 6: 247: 0: 0: 0, -1.80618, -21.0485 0/1: 18: 18,6: 12: 438: 6: 211: -13.9004,0, -32.5893
Anna1364

คำตอบ:


โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.