Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,318,527 | +AC | coding (1393/1446 nt) | tldD ← | predicted peptidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,318,524 | 1 | . | C | 85.7% | 167.8 / 23.4 | 49 | coding (1396/1446 nt) | tldD | predicted peptidase |
Reads supporting (aligned to +/- strand): ref base . (5/2); new base C (26/16); total (31/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01 | |||||||||||
* | REL606 | 3,318,524 | 2 | . | A | 85.7% | 157.8 / 21.7 | 49 | coding (1396/1446 nt) | tldD | predicted peptidase |
Reads supporting (aligned to +/- strand): ref base . (5/2); new base A (26/16); total (31/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01 |
TATCCTTCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGG > REL606/3318449‑3318588 || taTCCTTCTGAGAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGc‑‑c > 1:53552/1‑76 (MQ=255) taTCCTTCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGc‑‑c > 1:284207/1‑76 (MQ=255) ttCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑Acac < 1:210273/76‑3 (MQ=255) tCTGACAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑Acacg > 2:1351015/1‑73 (MQ=255) tCTGAAAATATTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑Acacg > 1:1927666/1‑73 (MQ=255) tCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑Acacg > 2:403725/1‑73 (MQ=255) tCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑Acacg > 2:1351011/1‑73 (MQ=255) tCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑Acacg > 1:1310224/1‑73 (MQ=255) aCTNAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGc > 1:1838405/1‑76 (MQ=255) gAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCaaa > 2:714562/1‑76 (MQ=255) aaaTTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAAc < 1:704175/76‑1 (MQ=255) attattACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACtt > 1:1401622/1‑76 (MQ=255) tgttaAGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAAACGGCAAACTTTg < 2:85731/74‑1 (MQ=255) tattaCGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTg < 2:1866457/76‑1 (MQ=255) tattaCGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTg > 1:216467/1‑76 (MQ=255) tagtaCGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTg < 2:1345736/73‑1 (MQ=255) taCGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGccc > 2:1841521/1‑76 (MQ=255) aCGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCct > 2:701884/1‑76 (MQ=255) cGCAGTACCGCCAACAGTCAGGTGATGNACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCctt > 2:1778/1‑76 (MQ=255) cAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCcttct > 1:411025/1‑76 (MQ=255) cAGNACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCcttct > 1:1834714/1‑76 (MQ=255) aGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCcttctt > 2:804133/1‑76 (MQ=255) aGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCcttctt > 2:1311553/1‑76 (MQ=255) gTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCttt < 1:1069199/76‑1 (MQ=255) gTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCttt < 2:1491489/76‑1 (MQ=255) gTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑ACANGCCGACCGGCAAACTTTGCCCTTCttt < 2:2619/76‑1 (MQ=255) tACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTa > 2:673420/1‑76 (MQ=255) ccGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCGTTAcc > 2:78680/1‑76 (MQ=255) cGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCg > 2:1684363/1‑76 (MQ=255) gCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGc < 2:1087375/76‑1 (MQ=255) aCAGTCAGGTTATCAACGGTCAGCATAGGCTGGCTAA‑‑G‑‑AC‑‑GCCGACCGCCAAACTTTGCCCTACTTTAGCGTAGac > 1:2041942/1‑76 (MQ=255) aGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGacac > 2:1849255/1‑76 (MQ=255) cAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCa > 1:1213665/1‑76 (MQ=255) cAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCa > 1:1477615/1‑76 (MQ=255) ggTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACg < 1:460068/76‑1 (MQ=255) ggTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACg < 2:530792/76‑1 (MQ=255) gTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGc < 1:1311554/76‑1 (MQ=255) gTTATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGc < 1:804134/76‑1 (MQ=255) ttATCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGcc > 2:1504264/1‑76 (MQ=255) aTCAACTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGt > 1:659542/1‑76 (MQ=255) aaCTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTAt < 2:776341/76‑1 (MQ=255) aaCTTTCAGCGTTGGCTCGC‑‑C‑‑CAGAG‑‑GCCAACCGGCAAANTTTGCACTTCTTTACNGCAGACCCCCACGCCGTTAt > 2:2019605/1‑76 (MQ=255) aCTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATc > 2:801577/1‑76 (MQ=255) cTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATcc > 1:1054326/1‑76 (MQ=255) cTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATcc < 1:209357/76‑1 (MQ=255) cTTTCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATca < 2:848952/76‑2 (MQ=255) tttCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCa < 2:876869/76‑1 (MQ=255) ttCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAg > 1:691041/1‑76 (MQ=255) tCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGt > 2:1172161/1‑76 (MQ=255) tCAGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGt > 1:580537/1‑76 (MQ=255) aGCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGttt > 2:1834165/1‑76 (MQ=255) gCGTTGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTc < 1:1103189/76‑1 (MQ=255) ttGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACGCCCACGCCGNTATCCAGTTTCAgg < 1:9528/76‑1 (MQ=255) ttGGCTGGC‑‑CCAC‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAgg > 1:1226596/1‑76 (MQ=255) || TATCCTTCTGAAAATACTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCAACTTTCAGCGTTGGCTGGC‑‑C‑‑C‑‑AC‑‑GCCAACCGGCAAACTTTGCCCTTCTTTACCGCAGACACCCACGCCGTTATCCAGTTTCAGG > REL606/3318449‑3318588 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |