Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,492,082 | 0 | T | G | 25.8% | 149.5 / 11.4 | 62 | V490G (GTG→GGG) | glpD | sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base T (8/38); new base G (16/0); total (24/38) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.69e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC > REL606/3492007‑3492157 | acgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggt < 1:1334844/76‑1 (MQ=255) gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtgg < 2:1673127/76‑1 (MQ=255) cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtggg > 1:917902/1‑75 (MQ=255) ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggggg > 1:974888/1‑74 (MQ=255) cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGt < 2:1092766/76‑1 (MQ=255) cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGt < 2:1231062/76‑1 (MQ=255) cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGt < 1:770318/76‑1 (MQ=255) tgtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGta < 1:84683/76‑1 (MQ=255) tgtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGta > 2:91956/1‑76 (MQ=255) gtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtat < 1:579172/76‑1 (MQ=255) tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtata < 1:517543/76‑1 (MQ=255) tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtata < 1:994248/76‑1 (MQ=255) tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGgata > 1:1588769/1‑76 (MQ=255) ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAc > 2:1641172/1‑76 (MQ=255) gCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACg < 1:367624/76‑1 (MQ=255) gCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACg < 1:676482/76‑1 (MQ=255) cGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgc < 1:1101208/76‑1 (MQ=255) gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgca < 2:1646070/76‑1 (MQ=255) ctCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgca < 2:566818/75‑1 (MQ=255) tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcag < 1:711190/76‑1 (MQ=255) cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagc < 2:1272892/76‑1 (MQ=255) cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagc < 1:1725977/76‑1 (MQ=255) cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACgcagc > 2:193643/1‑76 (MQ=255) cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAAgcagcaa > 2:1626938/1‑75 (MQ=255) caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGTATACGCAGCAGaa > 2:238687/1‑76 (MQ=255) aaaaCAAGGCATTTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAg > 2:290394/1‑76 (MQ=255) cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGt > 1:1477310/1‑76 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTc < 2:1868983/76‑1 (MQ=255) aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTc > 1:1913838/1‑76 (MQ=255) aGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCg < 1:808623/76‑1 (MQ=255) ggCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGc < 1:431308/76‑1 (MQ=255) cATGTGGCTAATTGCGGATCAGCAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTg < 2:961936/76‑1 (MQ=255) cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTg > 2:1305033/1‑76 (MQ=255) aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTgg < 1:861004/76‑1 (MQ=255) tgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTGGcgt > 1:1462776/1‑76 (MQ=255) ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGcggc > 1:369358/1‑74 (MQ=255) gCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcg < 2:1763296/76‑1 (MQ=255) cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcgt < 1:1903043/76‑1 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTg > 2:1258058/1‑76 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTg > 2:970438/1‑76 (MQ=255) tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCTGAAGAAGTCGGTGGCGGCGTg > 1:1202333/1‑76 (MQ=255) aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGGATACGCAGCAGAAGCTGTCGCTGGCGGCGTgg > 1:617417/1‑75 (MQ=255) aaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTggt > 1:939087/1‑74 (MQ=255) aTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAtt < 2:1006357/76‑1 (MQ=255) cgCGGATCAACATTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAttt < 1:436447/75‑1 (MQ=255) gATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATg > 2:1002053/1‑76 (MQ=255) gATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTGGCGGCGGGGTTTAATg > 2:547119/1‑76 (MQ=255) aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGt < 1:638308/76‑1 (MQ=255) tCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTa < 1:635318/76‑1 (MQ=255) tCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTa < 1:892178/76‑1 (MQ=255) caacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTaa < 1:192454/76‑1 (MQ=255) caacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTaa < 1:434770/76‑1 (MQ=255) aacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGGGTCGTGATTTAATGTAAg > 2:1684699/1‑76 (MQ=255) aacaaTCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGGGGGGTTTAATGTAAg > 2:53753/1‑76 (MQ=255) tctcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGc < 1:91273/76‑1 (MQ=255) tcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCa < 1:858260/76‑1 (MQ=255) cGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAg > 2:1717900/1‑76 (MQ=255) tgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCa < 1:581798/76‑1 (MQ=255) tgAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGGGATTTAATGTAAGGGGGCCAgggga > 2:1700077/1‑73 (MQ=255) gAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAg < 2:1535423/76‑1 (MQ=255) aGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTAAAGTAAGGGGGGCAAGGCAGa > 2:854790/1‑76 (MQ=255) tCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGAtt < 2:135220/76‑1 (MQ=255) tCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGAtt < 2:871525/76‑1 (MQ=255) cAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGAttt < 2:96645/76‑1 (MQ=255) gTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCa < 1:1471210/76‑1 (MQ=255) tGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGGCAGGATTCaa > 2:1156528/1‑76 (MQ=255) tGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGGCAGATTTCaa > 2:1758570/1‑76 (MQ=255) ggCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAAt < 1:501129/76‑1 (MQ=255) cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt > 2:282001/1‑76 (MQ=255) cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt < 1:857416/76‑1 (MQ=255) tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg > 2:687843/1‑76 (MQ=255) tggAGTATACGCAGCAGAGGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGGCAGATTTCAATCTGAcc > 2:11065/1‑76 (MQ=255) | ACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC > REL606/3492007‑3492157 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |