Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,492,082 | 0 | T | G | 52.6% | 33.3 / 56.8 | 38 | V490G (GTG→GGG) | glpD | sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base T (3/15); new base G (20/0); total (23/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.27e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC > REL606/3492007‑3492157 | acgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggt < 1:772806/76‑1 (MQ=255) acgCACTGTGTCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggt < 1:674444/76‑1 (MQ=255) gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTTAGTGGCtggtgg < 2:409350/76‑1 (MQ=255) cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGCGTGAGTCAGTGGCtgggggggg > 2:39701/1‑73 (MQ=255) ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAc > 2:237163/1‑76 (MQ=255) tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAGTAATCTCGTGTGAGTCAGTGNCTGGTGGAGTATACgcag < 1:1154838/76‑1 (MQ=255) cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACgcagc > 1:480760/1‑76 (MQ=255) aaaaCAAGGCATGTGGCTNAATGCGGATCAACAATCTCGCGTGAGTCAGTGGCTGGTGGGGGATACGCAGCAGAAg > 1:1154305/1‑76 (MQ=255) aaaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCCGCAGAAg > 2:67847/1‑76 (MQ=255) aaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCt < 1:415063/76‑1 (MQ=255) aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTg > 1:1077793/1‑76 (MQ=255) cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGt < 1:748361/76‑1 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGGTGTc > 1:559796/1‑76 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGAGGCTGGGGGGGGATACGCAGCAGAAGCTGTc > 2:200051/1‑76 (MQ=255) aGGCATGTGGCTAAATGTGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCg > 1:433457/1‑76 (MQ=255) aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTgg > 2:1075414/1‑76 (MQ=255) tgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgt < 1:502795/76‑1 (MQ=255) ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtc < 2:762489/76‑1 (MQ=255) ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGGCGATGGcggg > 2:812158/1‑74 (MQ=255) cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGcggcgt > 1:725719/1‑76 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGGGTCGTg < 2:933889/76‑1 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGgcggcggg > 1:991439/1‑74 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCCGGCGTCGTg > 2:1180315/1‑76 (MQ=255) aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGgcggcgggg > 1:616388/1‑73 (MQ=255) aaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTggt > 2:985883/1‑74 (MQ=255) tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAttt < 2:673830/76‑1 (MQ=255) gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTa < 1:834931/76‑1 (MQ=255) cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCCGGCGGCGTGGtttta > 1:909147/1‑74 (MQ=255) ggATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCATAAGCTGTCGCTGGCGTCGTGATTTAAt < 1:1187821/76‑1 (MQ=255) gATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCCGTCGCTGGCGGCGTGATTTAATg > 1:729352/1‑76 (MQ=255) caacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTaa < 1:638207/76‑1 (MQ=255) acaaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGCGGCTGGGCGCGTGGTTTTATGTAAgg > 1:57376/1‑76 (MQ=255) aaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTg < 1:39701/76‑1 (MQ=255) aTCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGGGGCGGGGTTTTATGTAAGgggg > 1:901097/1‑76 (MQ=255) ctcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGcc > 2:412080/1‑76 (MQ=255) tcGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTAATGTAAGGTGGCCa > 1:946862/1‑76 (MQ=255) gtgAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTAATGGAAGGGGGCCAgggc > 2:716820/1‑74 (MQ=255) tCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGAtt < 1:917535/76‑1 (MQ=255) cAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGGGGCGGGATTTTATGTAAGGGGGCCCGGGCCGAttt > 2:436393/1‑76 (MQ=255) tGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGNGATTTAATGTAAGGTGGCCAGGTCAGATTTCaa < 2:1016300/76‑1 (MQ=255) gCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATc < 1:916106/76‑1 (MQ=255) cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt > 2:11901/1‑76 (MQ=255) cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt < 1:622938/76‑1 (MQ=255) tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg > 1:1001731/1‑76 (MQ=255) tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg < 1:760945/76‑1 (MQ=255) tggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAcc < 1:1029966/76‑1 (MQ=255) | ACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC > REL606/3492007‑3492157 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |