| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,819,959 | 0 | A | C | 32.7% | 113.8 / 10.0 | 52 | T121P (ACC→CCC) | gdhA | glutamate dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base A (27/8); new base C (0/17); total (27/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.93e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAA > REL606/1819909‑1820005 | ccATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTa < 1:5921348/51‑1 (MQ=255)ccATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTa < 1:6876608/51‑1 (MQ=255) aTTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTAcc > 1:6115766/1‑51 (MQ=255) ttCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTAccc < 1:1555209/51‑1 (MQ=255) ttCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTAc < 1:5741247/49‑1 (MQ=255) tctcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCt < 1:5713073/51‑1 (MQ=255) tctcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCt < 1:5608386/51‑1 (MQ=255) ctcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTg < 1:4974799/51‑1 (MQ=255) ctcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTg > 1:6093290/1‑51 (MQ=255) tcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTGc < 1:6547860/51‑1 (MQ=255) tcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTGc > 1:4743983/1‑51 (MQ=255) tcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTGc < 1:1584913/51‑1 (MQ=255) tcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTGc > 1:762350/1‑51 (MQ=255) aaTTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTGCCGa > 1:6847041/1‑51 (MQ=255) aTTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTGCCGAt > 1:5538820/1‑51 (MQ=255) cTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTGCCGATGGGc > 1:7193057/1‑51 (MQ=255) cGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTGCCGATGGGCgg > 1:4259464/1‑51 (MQ=255) cTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTGCCGATGGGCggtgg > 1:853850/1‑51 (MQ=255) tttGAACAAACCTTCAAAAATGCCCTGACTACCCTGCCGATGGGCggtggt > 1:5178831/1‑51 (MQ=255) tttAAACAACCCTTCAAAAATGCCCTGACTCCCCTGCCGATGGGCggtggt < 1:2173239/51‑1 (MQ=255) ttGAACAAACCTTCAAAAATGCCCTGACTACCCTGCCGATGGGCGGTGGTa > 1:3804859/1‑51 (MQ=255) ttGAACAAACCTTCAAAAATGCCCTGACTACCCTGCCGATGGGCGGTGGTa > 1:2757986/1‑51 (MQ=255) tGAACAAACCTTCAAAAATGCCCTGACTCCCCTGCCGATGGGCGGTGGTaa < 1:3305026/51‑1 (MQ=255) gAACAAACCTTCAAAAATGCCCTGACTACCCTGCCGATGGGCGGTGGTaaa < 1:6865030/51‑1 (MQ=255) aaCAAACCTTCAAAAATGCCCTGACTCCCCTGCCGATGGGCGGTGGTaaa < 1:331567/50‑1 (MQ=255) aaCAAACCTTCAAAAATGCCCTGACTCCCCTGCCGATGGGCGGTGGTa < 1:4949759/48‑1 (MQ=255) aCAAACCTTCAAAAATGCCCTGACTCCCCTGCCGATGGGCGGTGGTAAAgg < 1:4057689/51‑1 (MQ=255) cAAACCTTCAAAAATGCCCTGACTCCCCTGCCGATGGGCGGTGGTaaa < 1:6501228/48‑1 (MQ=255) cAAACCTTCAAAAATGCCCTGACTACCCTGCCGATGGGCGGTGGTAAAggc > 1:3218305/1‑51 (MQ=255) aaaCCTTCAAAAATCCCCTGACTCCCCTGCCGATGGGCGGTGGTAAAg < 1:1240395/48‑1 (MQ=255) aaCCTTCAAAAATGCCCTGACTACCCTGCCGATGGGCGGTGGTAAAggcgg > 1:2379052/1‑51 (MQ=255) aaCCTTCAAAAATGCCCTGACTACCCTGCCGATGGGCGGTGGTAAAggcgg > 1:4387969/1‑51 (MQ=255) ccTTCAAAAATGCCCTGACTCCCCTGCCGATGGGCGGTGGTAAAGGCGGCa < 1:4853436/51‑1 (MQ=255) tNAAAAATGCCCTGACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCg > 1:2275710/1‑51 (MQ=255) tCAAAAATGCCCTGACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCg > 1:6311547/1‑51 (MQ=255) cAAAAATGCCCTGACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGa > 1:7109927/1‑51 (MQ=255) aaaaaTGCCCTGACTCCCCTGCCGATGGGCGGTGGTAAAggcggc < 1:2708485/45‑1 (MQ=255) aaaaaTGCCCTGACTCCCCTGCCGATGGGCGGTGGTAAAGGCGGCa < 1:4200933/46‑1 (MQ=255) aaaaaTGCCCTGACTCCCCTGCCGATGGGCGGTGGTAAAGGCGGCa < 1:6867660/46‑1 (MQ=255) aaaaaTGCCCCAACTCCCCTGCCGATGGGCGGTGGTAAAg < 1:1604699/40‑1 (MQ=37) aaaaaTCCCCTGACTCCCCTGCCGATGGGCGGTGGTa < 1:1745598/37‑1 (MQ=38) aaaTGCCCTGACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGAttt > 1:4234912/1‑51 (MQ=255) aaaTGCCCTGACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGAttt > 1:5976942/1‑51 (MQ=255) aaTGCCCTGACTACCCTGCCGATGGGCGGTGGTAAAggcg > 1:816170/1‑40 (MQ=255) aaTGCCCTGACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTc > 1:7062560/1‑51 (MQ=255) tGCCCTGACTCCCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGa < 1:6255323/51‑1 (MQ=255) gCCCTGACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGAt > 1:27273/1‑51 (MQ=255) cccTGACTCCCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTc < 1:3284932/47‑1 (MQ=255) cccTGACTCCCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTc < 1:2539844/47‑1 (MQ=255) cccTGACTCCCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGa < 1:6052902/49‑1 (MQ=255) ccTGACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATcc > 1:292968/1‑51 (MQ=255) ccTGACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATcc > 1:1629605/1‑51 (MQ=255) cTGACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCg > 1:758876/1‑51 (MQ=255) gACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGaa > 1:4770475/1‑51 (MQ=255) | CCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAA > REL606/1819909‑1820005 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |