Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,492,082 | 0 | T | G | 39.1% | 72.2 / 18.9 | 46 | V490G (GTG→GGG) | glpD | sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base T (4/24); new base G (18/0); total (22/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.59e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC > REL606/3492010‑3492157 | cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtgggggg > 2:1452399/1‑75 (MQ=255) ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg < 1:171599/76‑1 (MQ=255) gtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtat < 1:1050298/76‑1 (MQ=255) tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtata < 1:1382022/76‑1 (MQ=255) ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATAc > 2:1046121/1‑76 (MQ=255) ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAc > 1:1202032/1‑76 (MQ=255) gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAAgca > 1:229932/1‑76 (MQ=255) tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcag < 1:620340/76‑1 (MQ=255) tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcag < 2:153396/76‑1 (MQ=255) tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACgcag > 1:556162/1‑76 (MQ=255) gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagca < 1:1253170/76‑1 (MQ=255) cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGccgcag > 1:820373/1‑76 (MQ=255) acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATAAGCAGCCGa > 2:797901/1‑76 (MQ=255) caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGaa < 1:1456559/76‑1 (MQ=255) aaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCt > 2:536934/1‑76 (MQ=255) aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTg < 2:1424684/76‑1 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTc > 2:388905/1‑76 (MQ=255) gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCt < 1:310292/76‑1 (MQ=255) cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTg < 2:1390703/76‑1 (MQ=255) gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGctggcgctggcg > 1:1116038/1‑76 (MQ=255) tgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGGcgg > 2:112937/1‑75 (MQ=255) ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGcggc > 1:1388069/1‑74 (MQ=255) cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGGATACGCAGCAGAAGCTGTCGCTTGcgtcgt > 1:1359523/1‑76 (MQ=255) cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcgg > 1:220228/1‑75 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGGGGCGTg > 2:1193225/1‑76 (MQ=255) aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGa > 2:105533/1‑76 (MQ=255) aaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAt < 2:159531/76‑1 (MQ=255) tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGGGTCGTGGttt > 1:732804/1‑76 (MQ=255) cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTaa < 2:1048930/76‑1 (MQ=255) cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGGCGGCGGGGtttta > 2:1023956/1‑74 (MQ=255) cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGGGGCGGGATTTaa > 2:236898/1‑76 (MQ=255) ggATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAAt < 2:969912/76‑1 (MQ=255) aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGt < 1:11629/76‑1 (MQ=255) tCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTa < 1:688025/76‑1 (MQ=255) aacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAg < 1:1308699/76‑1 (MQ=255) aacaaTCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGGGGTCGTGGTTTAATGTAAg > 2:1410588/1‑76 (MQ=255) aaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTg > 1:927210/1‑76 (MQ=255) aTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTgg < 2:99703/76‑1 (MQ=255) ctcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGcc < 2:1276202/76‑1 (MQ=255) cGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAg < 1:126318/76‑1 (MQ=255) gtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTc < 1:662578/76‑1 (MQ=255) gtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTc < 2:632780/76‑1 (MQ=255) gtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTc < 2:864994/76‑1 (MQ=255) gTCAGTGGCGGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGAt < 1:908455/76‑1 (MQ=255) cAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGAttt < 2:1125800/76‑1 (MQ=255) cAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGAttt < 1:710556/76‑1 (MQ=255) tGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCaa > 2:648799/1‑76 (MQ=255) ggCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGCCAGGTTTCAAt > 1:942190/1‑76 (MQ=255) cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt > 2:165294/1‑76 (MQ=255) ggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGa < 1:867907/76‑1 (MQ=255) ggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGa < 1:740649/76‑1 (MQ=255) tggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCGGATTTCAATCTGAcc > 1:135289/1‑76 (MQ=255) tggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAcc > 2:1367784/1‑76 (MQ=255) | CCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC > REL606/3492010‑3492157 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |