Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 607,884 | G→A | 19.1% | R209* (CGA→TGA) | dck ← | Deoxyadenosine/deoxycytidine kinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 607,884 | 0 | G | A | 19.1% | 113.3 / 19.3 | 57 | R209* (CGA→TGA) | dck | Deoxyadenosine/deoxycytidine kinase |
Reads supporting (aligned to +/- strand): ref base G (22/23); new base A (7/4); total (29/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.06e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01 |
GTATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGTAT > USA300TCH1516_ALE/607746‑608023 | gtATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAg < 1:551110/141‑1 (MQ=255) atatGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCa < 1:119433/141‑1 (MQ=255) tgtgAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAtt < 2:468921/141‑1 (MQ=255) tgtgAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAtt < 2:468912/141‑1 (MQ=255) tAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTa < 1:116407/141‑1 (MQ=255) tGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCa > 2:82507/1‑141 (MQ=255) tGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCa > 2:21887/1‑141 (MQ=255) acaGAACTGATTAAGTAATCCATAATATCTACGCAATATACATTGAATATAAACACATACGCTATGTCTTTAGTCTTTCATTGCTTAGCTACTGTTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGAttt < 2:473248/139‑1 (MQ=255) aGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTaa < 1:12291/141‑1 (MQ=255) aGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTaa > 2:450775/1‑141 (MQ=255) aaTTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTTAATAATTAAAGTAATTGTGTGTATTATAGGATTTAAAg > 2:552603/1‑141 (MQ=255) aTTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGGTTTGAAGa > 1:78597/1‑141 (MQ=255) aaaaTAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATAGGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCggg > 1:584946/1‑141 (MQ=255) aataatCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGt < 1:313457/141‑1 (MQ=255) atCTATGGAATATACATAGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATATGATTTAAAGATTCGGGGGCCTTATGGATATc > 2:428276/1‑141 (MQ=255) cTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCat > 1:418288/1‑141 (MQ=255) aTGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCatat < 1:581626/141‑1 (MQ=255) ggAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTc < 2:239363/141‑1 (MQ=255) taCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATGTCATATTCATTAat > 2:594469/1‑141 (MQ=255) aaTATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATa > 1:311972/1‑141 (MQ=255) tataAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACg < 1:144424/141‑1 (MQ=255) tataAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACg < 2:295613/141‑1 (MQ=255) aaCATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACa < 1:239277/141‑1 (MQ=255) aaCATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACa > 1:349405/1‑141 (MQ=255) cATTTACGATATGTCTTTTGTATTTCTTCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAAc < 2:610816/141‑1 (MQ=255) atatACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACt > 1:68349/1‑141 (MQ=255) tataCGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTg > 1:201971/1‑141 (MQ=255) ataCGCTATGTCTTTAGTCTTTTATCGTGTATCTAGTTGTCGATATGTTTGAATAATTAGAGCGATTTTGTTTATCATAGGTTTTAAAGATTCGGGGTCCTTATGAATATCATATTCATTCATATTTATACGTACAACTgg > 2:245842/1‑141 (MQ=255) cGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGTTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACa < 2:394479/141‑1 (MQ=255) cGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACa < 2:536196/141‑1 (MQ=255) cGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACa < 1:318926/141‑1 (MQ=255) agatatgtatgTAGGCTTTTATTGTGTAACTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACa < 1:440718/138‑1 (MQ=255) cTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATg < 1:450775/141‑1 (MQ=255) cTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAAATGTACAtc > 2:320881/1‑140 (MQ=255) aTGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCa > 1:260433/1‑141 (MQ=255) aTGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCa > 1:568668/1‑141 (MQ=255) tGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCAt < 1:106666/141‑1 (MQ=255) gtctttagtctttaATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTGGTGTATCATAGGATGTAAAGATTCGGGGACATCATGGATATCATATTCATTAATATTGATACGTACATCTGGACGTGCAtt > 2:590083/1‑141 (MQ=255) tttAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTaaa < 1:328935/141‑1 (MQ=255) ttATTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAg < 2:506061/141‑1 (MQ=255) cTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTATAGCTAtt > 1:392220/1‑141 (MQ=255) tgtATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATc > 2:147444/1‑141 (MQ=255) ggtATCTACTTGTGGAGATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATc < 1:186041/140‑1 (MQ=255) gtATCTACTTGTCGATTTGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCg > 1:303484/1‑141 (MQ=255) tATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGt < 2:68349/141‑1 (MQ=255) tCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGCCATAGCGttt > 2:95765/1‑141 (MQ=255) tatGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAAt < 2:418288/141‑1 (MQ=255) ggTTTGAATAACTGCAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAg < 1:86598/140‑1 (MQ=255) gAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTCAATAGCtttt > 1:40675/1‑141 (MQ=255) aataatTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGGTATCATATTCATTAATATTGATACGTACAACTGGAGATGCATTCAAGCTATTAATCCAATCGTCATAGTGTTTAAATAGCttttt > 2:267907/1‑141 (MQ=255) ataatTCGAGCAATTTTGTTTATCACAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCttttac < 2:309448/141‑3 (MQ=255) atTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAGTTGCTTTTTCCAg > 2:20458/1‑141 (MQ=255) ttCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCatt > 1:141842/1‑139 (MQ=255) ttCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGt > 1:524394/1‑141 (MQ=255) ttCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGt > 1:71343/1‑141 (MQ=255) tCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAGTCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGTa > 2:269364/1‑141 (MQ=255) cGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGTAt < 2:307888/141‑1 (MQ=255) | GTATATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAGTAT > USA300TCH1516_ALE/607746‑608023 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |