Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,688,990 | 0 | A | C | 21.7% | 165.4 / 10.9 | 60 | T92P (ACC→CCC) | nth | DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III) |
Reads supporting (aligned to +/- strand): ref base A (39/8); new base C (0/13); total (39/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.94e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGG > REL606/1688915‑1689063 | cTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCaaaa < 1:1099370/76‑1 (MQ=255) tGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAAc > 2:1234806/1‑76 (MQ=255) tGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAAc < 1:1809576/76‑1 (MQ=255) ngggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAAc > 1:1555110/2‑76 (MQ=255) gggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAAcc > 2:206892/1‑76 (MQ=255) ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCt < 2:1536252/76‑1 (MQ=255) gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTg < 1:1205590/76‑1 (MQ=255) cGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGc > 1:1038719/1‑76 (MQ=255) gTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGcc < 1:2111305/76‑1 (MQ=255) tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt > 2:172830/1‑76 (MQ=255) tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt > 1:212342/1‑76 (MQ=255) aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt < 1:866820/76‑1 (MQ=255) aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt < 2:1120474/76‑1 (MQ=255) aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt > 2:2111894/1‑76 (MQ=255) aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt > 2:588323/1‑76 (MQ=255) tGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCt > 2:1776253/1‑76 (MQ=255) tGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCt > 2:1526921/1‑76 (MQ=255) gAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTg > 1:225749/1‑76 (MQ=255) aaaaCCTATATCAAAACGATGGGGCTTTAT‑AACAGCAAAGCAGAAAAAATCATCAAACCCTGCCGTATCTTGCTgg < 2:1838790/76‑1 (MQ=255) aaaCCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGa > 2:2088375/1‑76 (MQ=255) aaCCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGag > 2:2021108/1‑76 (MQ=255) ccTATATCAAACCGATTGGGCTTTAT‑AACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGagca < 2:1445041/76‑1 (MQ=255) cTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcag > 1:937675/1‑76 (MQ=255) cTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcag > 2:1210024/1‑76 (MQ=255) cTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcag > 1:1674535/1‑76 (MQ=255) cTATATCAAAACGATGGGGTTTTAT‑AACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagcag < 1:486619/76‑1 (MQ=255) atCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATa < 1:339978/76‑1 (MQ=255) atCAAAACGATGGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcagcnca < 2:1497651/76‑4 (MQ=255) aaaaaaCGATTGGGTTTTAT‑AACAACAAAGCAGAAAATATCAACAAACCCTGCCGTATCTTGCTGGAGCAGCATaa < 1:184831/74‑1 (MQ=255) cAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt > 2:301973/1‑76 (MQ=255) aaaaCGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATg < 2:1637636/76‑1 (MQ=255) aCGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCg < 2:788713/76‑1 (MQ=255) gATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAg > 2:800569/1‑76 (MQ=255) gATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAg > 1:1647200/1‑76 (MQ=255) aTTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg > 1:1549975/1‑76 (MQ=255) aTTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg > 1:1701529/1‑76 (MQ=255) ttGGGGTTTAAAAACAGAAAAGAAAAAAATACCAAAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg < 1:927627/76‑1 (MQ=255) tGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt > 2:1981600/1‑76 (MQ=255) gggCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTc > 1:701596/1‑76 (MQ=255) gCTTTAA‑AACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCg < 1:1485142/76‑1 (MQ=255) cTTTAT‑AACAGAAAAGCAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg < 2:1893010/76‑1 (MQ=255) tttAT‑ACCAAAAAAGCAAAAAATTCCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa < 2:746021/76‑1 (MQ=255) ttAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa > 1:225446/1‑76 (MQ=255) ttAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa > 2:2116738/1‑76 (MQ=255) tat‑aACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg > 1:771262/1‑76 (MQ=255) tat‑aACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg > 2:1948025/1‑76 (MQ=255) tat‑aACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg > 1:1909494/1‑76 (MQ=255) at‑aACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGa < 2:321536/76‑1 (MQ=255) t‑aACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGAt > 1:1060735/1‑76 (MQ=255) t‑aACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGAt > 1:1876021/1‑76 (MQ=255) aaaaCAGAAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGa < 1:11936/74‑1 (MQ=255) aaaCAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGAt < 1:1414336/75‑1 (MQ=255) aaCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATc > 1:444500/1‑76 (MQ=255) aCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCg > 1:270410/1‑76 (MQ=255) gCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgc > 1:504776/1‑76 (MQ=255) gCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgc > 1:357409/1‑76 (MQ=255) aaaaCATCACCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGa < 2:1888430/76‑1 (MQ=255) aaaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGaa > 1:1198559/1‑76 (MQ=255) aaTATCAAAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAg < 1:1092425/76‑1 (MQ=255) atatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGc > 1:76167/1‑76 (MQ=255) tatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGcc > 1:744351/1‑76 (MQ=255) atcatcAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGccc > 1:2039337/1‑76 (MQ=255) catcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg < 2:941155/76‑1 (MQ=255) tcaAAACCTGCCGTATCATGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGcc > 2:1019983/1‑76 (MQ=255) caAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGccc > 1:1074592/1‑76 (MQ=255) aaacCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg < 2:1909494/76‑1 (MQ=255) aaacCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg < 2:771262/76‑1 (MQ=255) aaaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg > 2:1129368/1‑76 (MQ=255) aaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCgg > 1:1855970/1‑76 (MQ=255) aaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCgg > 2:1803421/1‑76 (MQ=255) | CTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTAT‑AACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGG > REL606/1688915‑1689063 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |