Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,526,164 | 0 | G | C | 31.1% | 102.8 / 16.8 | 45 | S147S (TCG→TCC) | hyfD | hydrogenase 4 membrane subunit |
Reads supporting (aligned to +/- strand): ref base G (5/26); new base C (14/0); total (19/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.97e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTT > REL606/2526116‑2526214 | aGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGtt < 1:6068784/51‑1 (MQ=255) aGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGtt < 1:2286923/51‑1 (MQ=255) cgcAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTgg < 1:5892074/51‑1 (MQ=255) cgcAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTgg < 1:5785756/51‑1 (MQ=255) cAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTgggg < 1:6360194/51‑1 (MQ=255) gCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGtt < 1:1910065/51‑1 (MQ=255) cGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCCCTgggg > 1:6794298/1‑48 (MQ=38) gTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCCCTgggg > 1:976145/1‑47 (MQ=38) ttCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCCTTgggg > 1:5073486/1‑46 (MQ=255) ttCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCCCTgggg > 1:3141471/1‑46 (MQ=38) tCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCCCTgggg > 1:732520/1‑45 (MQ=38) cAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATc < 1:7134381/51‑1 (MQ=255) gCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCtt < 1:3470895/51‑1 (MQ=255) gCACTAAAAGCGTTACTTATCACTCATATCGGCTCCCTgggg > 1:3451492/1‑42 (MQ=38) cACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTg < 1:3312735/51‑1 (MQ=255) aCTAAAAGCGTTACTTATCACTCATATCGGCTCCCTgggg > 1:1750216/1‑40 (MQ=38) cTAAAAGCGTTACTTATCACTCATATCGGCTCCTTGGGGGTGTATCtt > 1:5892947/1‑48 (MQ=255) tAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTgccg < 1:625988/51‑1 (MQ=255) aaaaGCGTTACTTATCACTCATATCGGCTCCTTgggg > 1:5249712/1‑37 (MQ=255) aaaaGCGTTACTTATCACTCATATCGGCTCCCTgggg > 1:843286/1‑37 (MQ=38) aaaGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTgccgcc < 1:3778384/51‑1 (MQ=255) aaGCGTTACTTATCACTCATATCGGCTCCCTgggg > 1:3070800/1‑35 (MQ=37) aGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTgccgccgc < 1:5536609/51‑1 (MQ=255) gCGTTACTTATCACTCATATCGGCTCCTTtgggg > 1:3151871/1‑33 (MQ=25) gCGTTACTTATCACTCATATCGGCTCCCTgggg > 1:235554/1‑33 (MQ=37) cGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCa < 1:6692043/51‑1 (MQ=255) cGTTACTTATCACTCATATCGGCTCCTTgggg > 1:6381741/1‑32 (MQ=255) gTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCAc < 1:2276960/51‑1 (MQ=255) tACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGc < 1:364898/51‑1 (MQ=255) aCTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCt < 1:2536735/51‑1 (MQ=255) ttATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGt < 1:61100/51‑1 (MQ=255) cACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGttttt < 1:3807853/51‑1 (MQ=255) aCTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTg < 1:749252/51‑1 (MQ=255) tCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGttttt > 1:2086087/1‑43 (MQ=255) tCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCg > 1:5538554/1‑51 (MQ=255) ggCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGa < 1:6180331/51‑1 (MQ=255) ggCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGa < 1:3676147/51‑1 (MQ=255) gCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGaa > 1:67534/1‑51 (MQ=255) gCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGaa < 1:7153726/51‑1 (MQ=255) cTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAAc > 1:1551074/1‑51 (MQ=255) cTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAAc < 1:2069577/51‑1 (MQ=255) cTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAAc < 1:908953/51‑1 (MQ=255) tCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACg < 1:11089/51‑1 (MQ=255) tCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACg < 1:1701684/51‑1 (MQ=255) cGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGt > 1:5672081/1‑51 (MQ=255) gTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGtt > 1:4192640/1‑51 (MQ=255) | AGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTT > REL606/2526116‑2526214 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |