Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,492,082 | 0 | T | G | 53.3% | 21.9 / 35.6 | 31 | V490G (GTG→GGG) | glpD | sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base T (3/11); new base G (16/0); total (19/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.66e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTG > REL606/3492009‑3492154 | gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtgg < 2:1081916/76‑1 (MQ=255) cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtgggggg > 1:636425/1‑75 (MQ=255) cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtgggggg > 1:527955/1‑75 (MQ=255) tgtgGAGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGCGGGGta > 1:80209/1‑76 (MQ=255) tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtata < 1:599760/76‑1 (MQ=255) ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAc > 1:1080395/1‑76 (MQ=255) gCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACg < 1:815186/76‑1 (MQ=255) gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAATCAGTGGCTGGGGGGGTATACgcagca > 2:65226/1‑76 (MQ=255) caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGTGTATACGCAGCAGaa < 1:432855/76‑1 (MQ=255) aaaCAAGGCCTGTGGCTCGATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGNTANGCAGCAGAAGc > 2:192762/1‑76 (MQ=255) cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGt < 2:803206/76‑1 (MQ=255) cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAAGCAGCAGAAGCTGt > 1:391155/1‑76 (MQ=255) cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAAGCAGCAGAAGCTGt > 2:565041/1‑76 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTc < 2:1027612/76‑1 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTc > 2:377709/1‑76 (MQ=255) aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCCgg > 2:344229/1‑76 (MQ=255) aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGCATCCGGAGCAGAAGCTCGCGCTgg > 2:344242/1‑76 (MQ=255) tgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGcgt > 1:785187/1‑76 (MQ=255) ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCCGTCGCTGGcggc > 2:652135/1‑74 (MQ=255) ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGcggc > 1:643857/1‑74 (MQ=255) gCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcg < 1:729322/76‑1 (MQ=255) cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcgg > 1:942527/1‑75 (MQ=255) aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGTATACGCAGCAGAAGCTGTCGGTGGCGGCGTgg > 2:83345/1‑75 (MQ=255) aaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTGGCGTCGGGAt > 2:139507/1‑76 (MQ=255) aTGCGGATCAANTATCTCGTGTGAGTCAGTGGCTAGCGGCGTATNCGCAGCAGAAGCTGTCGGTGGTGTCGTGAtt < 1:148181/76‑1 (MQ=255) cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTaa < 2:836103/76‑1 (MQ=255) tgATCAACAATCTCGAGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTGGCGGCGGGGTTTTAt > 2:585624/2‑76 (MQ=255) tCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTa < 1:712098/76‑1 (MQ=255) caacaaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTaa > 2:148756/1‑76 (MQ=255) caaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGGAGCAGAAGCTGTCGCTGGAGGCGTGATTTAATGTAAggg > 2:965802/1‑75 (MQ=255) caaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAggg > 2:488852/1‑75 (MQ=255) aaTCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGGGTCGTGATTTAATGTAAGGTg > 2:454487/1‑76 (MQ=255) aTCTCGTGTGAGTCAGTGGCGGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGCGATTTAATGTAAGGTgg < 1:1199571/76‑1 (MQ=255) ctcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGCGATTTAATGTAAGGTGGcc < 1:1025969/76‑1 (MQ=255) gtgAGTCAGTGGCGTGTGACGTATACGCATCAGACGCTGTTGCTGTCGTCGTGATTTAATGTAAGGTAGCCAGGTc < 2:23254/76‑1 (MQ=255) tGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGGCAGATTTCaa > 2:1032254/1‑76 (MQ=255) cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt < 2:1068574/76‑1 (MQ=255) tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGTCCAGGTCAGATTTCAATCTg > 2:969183/1‑76 (MQ=255) | GCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTG > REL606/3492009‑3492154 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |