Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,476,217 | 0 | C | A | 73.3% | 45.4 / 23.8 | 30 | G75V (GGT→GTT) | yheO | putative DNA‑binding transcriptional regulator YheO |
Reads supporting (aligned to +/- strand): ref base C (4/4); new base A (6/16); total (10/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.84e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATT > NC_000913/3476142‑3476297 | tAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATc < 2:610062/81‑1 (MQ=255) gCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCtt < 2:610064/81‑1 (MQ=255) gCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCtt < 2:610063/81‑1 (MQ=255) tttAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCg > 1:610098/1‑81 (MQ=255) tttAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCg < 2:610065/81‑1 (MQ=255) tAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCgg < 2:610066/80‑1 (MQ=255) aaCGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGt < 1:610067/81‑1 (MQ=255) ctgctATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGt > 1:610111/1‑81 (MQ=255) ctgctATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGt > 1:610120/1‑81 (MQ=255) ctATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCa > 2:610099/1‑81 (MQ=255) tCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCAcc < 2:610068/80‑1 (MQ=255) cGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAAt > 2:610123/1‑81 (MQ=255) atCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAAt > 2:610124/1‑81 (MQ=255) tGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCt < 1:610069/81‑1 (MQ=255) tGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCt < 1:610070/81‑1 (MQ=255) gcagcaTACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCt > 1:610113/1‑80 (MQ=255) cagcaTACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTg < 2:610071/80‑1 (MQ=255) gcaTACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAAc < 1:610072/81‑1 (MQ=255) cGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACAtttt < 2:610073/81‑1 (MQ=255) tAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAg < 2:610074/81‑1 (MQ=255) aGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACAttt > 1:610114/1‑77 (MQ=255) cgcCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATc < 2:610075/81‑1 (MQ=255) aGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGc < 1:610076/81‑1 (MQ=255) aGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGc < 2:610077/81‑1 (MQ=255) ggTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCa < 2:610078/81‑1 (MQ=255) gTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCaa < 1:610079/81‑1 (MQ=255) gTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCa > 2:610117/1‑80 (MQ=255) cAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCa > 1:610122/1‑78 (MQ=255) aTTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGc < 1:610080/81‑1 (MQ=255) gCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCaaaa < 2:610081/81‑1 (MQ=255) ccAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGAtt > 2:610129/1‑81 (MQ=255) | TAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATT > NC_000913/3476142‑3476297 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |