| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,250,830 | 0 | T | . | 66.7% | 14.0 / 33.6 | 21 | coding (9/1071 nt) | dhaK | dihydroxyacetone kinase subunit K |
| Reads supporting (aligned to +/- strand): ref base T (4/3); new base . (8/6); total (12/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGAT > NC_000913/1250756‑1250907 | aGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAAtttttt > 2:226485/1‑80 (MQ=255) aGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCattat < 2:226450/80‑1 (MQ=255) aGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAA‑TTTTTCATTAttt > 1:226499/1‑81 (MQ=255) ggatgCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAA‑TTTTTCATTATTTTGCTc > 2:226472/1‑80 (MQ=255) gatgCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAA‑TTTTTCATTATTTTGCTcc < 2:226451/80‑1 (MQ=255) cTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAAt < 2:226452/81‑1 (MQ=255) tCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGTTCAA‑TTTTTCATTATTTTGCTCCAGCAATTACg > 1:226477/1‑81 (MQ=255) gCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAA‑TTTTTCATTATTTTGCTCCAGCAATTACgg > 2:226482/1‑80 (MQ=255) aGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAA‑TTTTTCATTATTTTGCTCCAGCAATTACGGTAgg > 1:226471/1‑81 (MQ=255) tCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAA‑TTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGc < 2:226453/81‑1 (MQ=255) cAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGga > 1:226476/1‑81 (MQ=255) ttgttCGTCCAGTACGTCTTGCACATCATTGTTCAA‑TTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGgatgat > 1:444689/1‑81 (MQ=255) ttgttCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGgatga > 2:226486/1‑81 (MQ=255) ttCGTCCAGTACGTCTTACACATCATTGATCAA‑TTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGt < 1:226454/80‑1 (MQ=255) tCCAGTACGTCTTGCACATCATTGATCAA‑TTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAAc < 2:226455/81‑1 (MQ=255) gTACGTCTTGCACATCATTGATCAA‑TTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGaca < 2:226456/81‑1 (MQ=255) gCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGAc > 2:226483/1‑80 (MQ=255) cacaTCATTGATCAA‑TTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCac > 1:226487/1‑81 (MQ=255) acaTCATTGATCAA‑TTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCact > 2:226481/1‑81 (MQ=255) atTGATCAA‑TTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCactgact < 1:226457/80‑1 (MQ=255) aTCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGcc < 2:226458/80‑1 (MQ=255) aaTTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGAt > 2:226492/1‑80 (MQ=255) | AGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGAT > NC_000913/1250756‑1250907 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |