Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,220,009 | 0 | G | T | 35.3% | 47.8 / 27.1 | 34 | V472F (GTT→TTT) | aceK | isocitrate dehydrogenase kinase/phosphatase |
Reads supporting (aligned to +/- strand): ref base G (13/9); new base T (5/7); total (18/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.75e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTCGCCACTG > NC_000913/4219862‑4220157 | aaTATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGt > 8:126667/1‑149 (MQ=255) atCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGtttt > 2:85595/1‑149 (MQ=255) gctggAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGTTTTTTTa < 5:68268/147‑1 (MQ=255) gctggAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGTTTTTTTATg < 5:128065/149‑1 (MQ=255) ctggAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGa < 4:81372/149‑1 (MQ=255) gTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTAGGATGAAATTTGc < 1:85595/149‑1 (MQ=255) aaTACGGTAACGCTATTCGCCAGCGTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGTTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGAc < 4:119608/149‑1 (MQ=255) tACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACAt > 1:165790/1‑149 (MQ=255) aCGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCccgccg < 2:165790/149‑1 (MQ=255) cGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCccgccgc > 1:57954/1‑149 (MQ=255) gCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGTTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTAt > 2:102830/1‑149 (MQ=255) cGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTATGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGAc > 6:88978/1‑148 (MQ=255) aaCATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGa > 6:113266/1‑77 (MQ=255) aaCATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGa < 5:113266/77‑1 (MQ=255) ttttCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGTTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGc > 5:170141/1‑149 (MQ=255) ggCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTg > 5:153164/1‑149 (MQ=255) gCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGTTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTAGCCAGCGAACCGAgg > 8:7869/1‑149 (MQ=255) gCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTg < 6:153164/148‑1 (MQ=255) gACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTa > 2:174346/1‑149 (MQ=255) gACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTa > 5:162706/1‑149 (MQ=255) aTGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGt > 5:123557/1‑142 (MQ=255) tGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGc < 7:126667/149‑1 (MQ=255) gTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGTTTTTTTATGATTACGATGAAATTTGCTACAt < 2:133974/69‑1 (MQ=255) gTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGTTTTTTTATGATTACGATGAAATTTGCTACAt > 1:133974/1‑69 (MQ=255) tttAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGTTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAACGtc < 6:170141/147‑1 (MQ=255) tAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAAc < 2:31175/119‑1 (MQ=255) tAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTAAGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAAc > 1:31175/1‑119 (MQ=255) aaaaaCTTCGGTGTCACCCGTCACGGGCGTGTGTTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATc < 2:152529/106‑1 (MQ=255) aaaaaCTTCGGTGTCACCCGTCACGGGCGTGTGTTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATc > 1:152529/1‑106 (MQ=255) cccGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTgg > 6:168297/1‑119 (MQ=255) ggCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGttt < 6:162706/149‑1 (MQ=255) gCGTGTGTTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTc < 1:102830/149‑1 (MQ=255) gtgtGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTCGc < 1:174346/149‑1 (MQ=255) gtgGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTCGCCa > 2:43800/1‑149 (MQ=255) gTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTCGCCACTg > 6:1225/1‑149 (MQ=255) | AATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTCGCCACTG > NC_000913/4219862‑4220157 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |