| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,688,990 | 0 | A | C | 31.4% | 57.5 / 13.3 | 35 | T92P (ACC→CCC) | nth | DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III) |
| Reads supporting (aligned to +/- strand): ref base A (17/7); new base C (0/11); total (17/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.06e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGG > REL606/1688917‑1689063 | gggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAcc < 2:1917717/76‑1 (MQ=255) cGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGc > 1:446680/1‑76 (MQ=255) ntGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCg > 1:1962683/2‑76 (MQ=255) tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt > 2:176405/1‑76 (MQ=255) gAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTa > 2:1227134/1‑76 (MQ=255) aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt < 1:1659596/76‑1 (MQ=255) aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt > 2:452797/1‑76 (MQ=255) aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt > 2:1838128/1‑76 (MQ=255) gTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGc > 2:643111/1‑76 (MQ=255) tGAAAACCTATATCAAAACGATTGGGCTTTAAAACAGCAAAGCAAAAAATATCATCAAAACCTGCCGTATCTTGCt < 2:958342/76‑1 (MQ=255) aaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGag < 2:993642/76‑1 (MQ=255) aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagc > 1:389885/1‑76 (MQ=255) atatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAgcagca > 1:132556/1‑76 (MQ=255) taacAAAACGATTGGGCTTTATAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCAt < 1:773963/73‑1 (MQ=255) cAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt > 1:759584/1‑76 (MQ=255) cAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt < 2:1024748/76‑1 (MQ=255) cAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt < 1:1402896/76‑1 (MQ=255) aaaaaGGATTGGGCTTAATAACAGAAAAGCAGAAAAAATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAAt < 2:1212658/75‑1 (MQ=255) aaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATgg < 1:1907296/76‑1 (MQ=255) aaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGc < 1:77230/76‑1 (MQ=255) aaCGATTGGGCTTTAAAACAGAAAAGCAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGc < 1:1613392/76‑1 (MQ=255) aTTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg < 2:592401/76‑1 (MQ=255) tggggCTTTATAACAGCAAAGAAGAAAATATCAAAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt < 2:642842/74‑1 (MQ=255) ggggTTTTATAACAGCAAAGAAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt < 1:963927/75‑1 (MQ=255) gCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCAAAATGGCGGGGTTCCg > 2:732247/1‑76 (MQ=255) cTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg > 1:436444/1‑76 (MQ=255) tttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa > 2:1156947/1‑76 (MQ=255) ttATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa < 2:1520844/76‑1 (MQ=255) ttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa > 1:1997944/1‑76 (MQ=255) taaaaCAGAAAAGCAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg < 2:881523/73‑1 (MQ=255) cataaaaTATCATCAAAACCTTCCGTCTCTTGCTGCAGCAGCATCATGGCGAGCATCCGCTACAACGTTCTGCGCt > 2:1809153/4‑76 (MQ=255) caaaaaaTATCCTAAACCCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCt < 1:1247526/73‑1 (MQ=255) aaaaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGCTCCGGAAGATCGACCTGCGCtcca > 1:744707/1‑73 (MQ=255) aaaaTACCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGa < 2:1558374/76‑1 (MQ=255) catcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg > 1:258362/1‑76 (MQ=255) catcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg > 1:310836/1‑76 (MQ=255) catcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg > 1:1206093/1‑76 (MQ=255) caAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGccc < 1:465686/76‑1 (MQ=255) caAAACCTGCCGTATCTTGCTGGAGCAGCATAATGACGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGccc > 2:1811599/1‑76 (MQ=255) aaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCgg > 2:836003/1‑76 (MQ=255) | GGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGG > REL606/1688917‑1689063 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |