Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

missing avg_distance_to_effective_3p_end #63

Open
stemangiola opened this issue May 9, 2019 · 0 comments
Open

missing avg_distance_to_effective_3p_end #63

stemangiola opened this issue May 9, 2019 · 0 comments

Comments

@stemangiola
Copy link

In analysis some BAM files I get these kind of results (one line shown)

chr17
748892
T
1
=:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00
A:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00
C:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00
G:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00
T:1:0.00:37.00:0.00:1:0:0.00:0.00:0.00:0:0.00:0.00:0.00
N:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00

Apparently none of the key info to use for example sinvict are present.

Is this a pathological result?

I aligned the reads with bwa

M04416:61:000000000-C2G5K:1:1101:20297:9269     77      *       0       0       *       *       0       0       GAGGGGATTTATGTTTGGTGCTATTGCTGGCGGTATTGCTTCTGCTCTTGCTGGTGGCGCCATGTCTAAATTGTTTGGAGGCGGTCAAAAAGCCGCCTCCGGTGGCATTCAAGGTGATGTGCTTGCTACCGATAACAATACTGTAGGCATG >1>A1>AAFFFFFAFFGAGGFFDFAGHGHHGCCEEEGGHFHGHHFA11F1GFFF0F?FEAFFAGHHHFHHGHBGFH0/GGEEGGEGEBHHH/FF@?EACHCC/B?FGHH2?GGB1FFB2FBGHHBGHHFHGGGGGHHGHFHHFGHBGGHHH AS:i:0  XS:i:0
M04416:61:000000000-C2G5K:1:1101:20297:9269     141     *       0       0       *       *       0       0       CAGCGGCAGACTTGCCACCAAGTCCAACCAAATCAAGCAACTTATCAGAAACGGCAGAAGTGCCAGCCTGCAACGTACCTTCAAGAAGTCCTTTACCAGCTTTAGCCATAGCACCAGAAACAAAACTAGGGACGGCCTCATCAGGGTTAGG A1>AA?@DDAFFGGGGGGG00BGFHH1BAFGGFD1AAF0FGFGFFGHHHHCGGECECGBGHDFHFEAEHHHHHG@CGEGDF@D1B1FHHHBGGHGHBGBE>FHBEGHHHH1>FGHHFFHHHGGHHHGHFHHGGGGGGGHHHHHHGHGHHHG AS:i:0  XS:i:0
M04416:61:000000000-C2G5K:1:1101:12893:9276     83      chr17   10123357        60      8S66M11S        =       10123357        -66     CCCCCTCCTCTGTTGCTGCTGGTCCAGGGACCACACTTTAAGAAGCCCTACCCTAGAATCTGAAGGTTGGAAAGCCATAGGATTC   GGGGFHGHHHHHHGHHHHHHHHHHHHHIHHGHHHGHHHHHHHHHHGGGHHGHHHHHHHHHHHHGGGGGGGGGGCFFFFFFCCBBB   NM:i:1  MD:Z:58C7       AS:i:61 XS:i:31
M04416:61:000000000-C2G5K:1:1101:12893:9276     163     chr17   10123357        60      8S66M11S        =       10123357        66      CCCCCTCCTCTGTTGCTGCTGGTCCAGGGACCACACTTTAAGAAGCCCTACCCTAGAATCTGAAGGTTGGAAAGCCATAGGATTC   ABCCCCCCFFFFGGGGGGGGGGGHHHGHGGGGHHHGHHHHHHHHHHHHHHHHHHHGHHHHHHHGHHHHHGHHHHGHHHHHHHHHH   NM:i:1  MD:Z:58C7       AS:i:61 XS:i:31

What is the best command to align the reads before using bam-readcount?

Thanks

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant