Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,688,990 | 0 | A | C | 28.3% | 100.1 / 10.9 | 46 | T92P (ACC→CCC) | nth | DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III) |
Reads supporting (aligned to +/- strand): ref base A (28/5); new base C (0/13); total (28/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.42e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCG > REL606/1688915‑1689062 | cTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa < 1:1327300/76‑1 (MQ=255) cTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa > 2:447066/1‑76 (MQ=255) ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt < 1:60529/76‑1 (MQ=255) ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt > 2:1871183/1‑76 (MQ=255) gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg < 1:538637/76‑1 (MQ=255) gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg < 2:1360106/76‑1 (MQ=255) gCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg < 1:483719/76‑1 (MQ=255) gTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGcc > 2:733089/1‑76 (MQ=255) ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCg < 2:792062/76‑1 (MQ=255) tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt > 1:1842384/1‑76 (MQ=255) tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt > 2:1958255/1‑76 (MQ=255) gAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTa > 2:207306/1‑76 (MQ=255) gAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTa > 2:662549/1‑76 (MQ=255) ggggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCt > 1:698749/1‑76 (MQ=255) gggTAAAACCCTATATAAAAACGTTTGGGCTTTATAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCtt < 2:713853/76‑1 (MQ=255) gggTAAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCtt < 2:1546248/76‑1 (MQ=255) gTGAAACCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGc < 1:920762/76‑1 (MQ=255) gTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGc > 2:423919/1‑76 (MQ=255) gAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTg < 2:1194961/76‑1 (MQ=255) gAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTg > 2:760836/1‑76 (MQ=255) aaaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTgg > 1:1756158/1‑76 (MQ=255) aaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGa > 2:926542/1‑76 (MQ=255) aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagc > 2:390230/1‑76 (MQ=255) ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagca > 1:1976359/1‑76 (MQ=255) ccTATATCAAAACGATGGGCCTTTATAACAACAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagca < 2:748995/76‑1 (MQ=255) cTATATAAAAACGATTGGGTTTTATAACAGCAAAGAAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGagcag < 1:1871184/76‑1 (MQ=255) tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcagc > 1:1096620/1‑76 (MQ=255) atatCAAAACGATTGGGCTTTATAACAGAAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAgcagca < 2:1974610/76‑1 (MQ=255) atatCAAAACGATTGGGCTTTAAAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAgcagca < 1:1712037/76‑1 (MQ=255) cAAAACGATTGGGCTTTATAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAAt < 2:602834/76‑1 (MQ=255) aTTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg > 2:1643883/1‑76 (MQ=255) ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt > 2:532269/1‑76 (MQ=255) tGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt > 1:537435/1‑76 (MQ=255) gggCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTc < 2:1548298/76‑1 (MQ=255) gCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCg > 1:521317/1‑76 (MQ=255) cTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg > 1:1863729/1‑76 (MQ=255) cTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg > 2:558392/1‑76 (MQ=255) tttATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa < 1:1319116/76‑1 (MQ=255) tttATAACAGAAAAGAAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa < 2:1612224/76‑1 (MQ=255) aCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCg > 2:1484689/1‑76 (MQ=255) aaaGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgctg > 2:1228084/1‑76 (MQ=255) cAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCt > 1:346067/1‑76 (MQ=255) aaaaaaTATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCtt < 1:1858944/74‑1 (MQ=255) aaaTATCAACAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGaa < 1:1295273/76‑1 (MQ=255) aaTATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAg < 1:1722514/76‑1 (MQ=255) atcatcAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGccc > 1:1601830/1‑76 (MQ=255) atcatcAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGccc > 2:892742/1‑76 (MQ=255) tcatcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCt < 2:1228435/76‑1 (MQ=255) catcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg > 2:551784/1‑76 (MQ=255) catcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg > 2:622648/1‑76 (MQ=255) atcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGc < 2:515226/76‑1 (MQ=255) caAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGccc > 1:192706/1‑76 (MQ=255) aaaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg > 2:807930/1‑76 (MQ=255) aaaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGGTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg > 2:1635421/1‑76 (MQ=255) | CTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCG > REL606/1688915‑1689062 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |