Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,017,405 | 0 | A | C | 42.1% | 17.2 / 10.9 | 19 | D449A (GAC→GCC) | ygfK | putative oxidoreductase, Fe‑S subunit |
Reads supporting (aligned to +/- strand): ref base A (11/0); new base C (0/8); total (11/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.32e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATT > NC_000913/3017345‑3017469 | tAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGGGGCAGAAGACCTg > 1:3262877/1‑65 (MQ=255) aCCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGcc > 1:502/1‑65 (MQ=255) ccATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCg > 1:276058/1‑65 (MQ=255) ccATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCGGCCg > 1:611886/1‑65 (MQ=255) ccATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGAACTGCCg > 1:453686/1‑65 (MQ=255) cATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGc > 1:3226340/1‑65 (MQ=255) ggAATACACCCAGAAACGCTGAAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGCCCTGCCGCTGa < 1:1113541/65‑1 (MQ=255) gAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGat > 1:1215494/1‑64 (MQ=255) aTACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCg > 1:3877922/1‑65 (MQ=255) tACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGa > 1:3325364/1‑65 (MQ=255) acacCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGAc > 1:1165544/1‑65 (MQ=255) cacCCAGAAACACTGGAAGCCAGAAGAGCGTATGGAAGGGGCAGAAGACCTGCCGCTGACCGACt < 1:1303682/65‑1 (MQ=255) cacCAAAAAACCCTGAAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGCCCTGCCGCTGACCGACt < 1:3855455/65‑1 (MQ=255) agcccAGAAACACTGGAAGCCAGAAGAGCGTATTAAAGTGGCAGAAGCCCTGCCGCTGACCGACt < 1:1187223/63‑1 (MQ=255) acCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTg > 1:1113272/1‑65 (MQ=255) acCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTg > 1:18759/1‑65 (MQ=255) acCCAGAAACACTGGAAGCCAGAAGAGCGGATTGAAGTGGCAGAAGAACTGCCGCTGACCGACTg > 1:1460683/1‑65 (MQ=255) aaaaaCACGGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGCCCTGCCGCTGACCGACTGCTAc < 1:2548026/63‑1 (MQ=255) aaaaACTGGAAGCCAGAAGAGCGTTTTGAGGTGGCAGAAGCCCTGCCGCTGACCGACTGCTACGt < 1:3483701/65‑1 (MQ=255) aaaCACGGGAAGCAAGAAGAGGGTATGGAGGTGGCAGAGGCCCTGCCGCTGACCGACTGCTACGt < 1:3018447/65‑1 (MQ=255) acacTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTg > 1:2614272/1‑65 (MQ=255) gAAGCCAGAAGAGCGTATTGAAGGGGCAGAAGCCCTGCCGCTGACCGACTGCTACGTTGCCCCCt < 1:2889379/65‑1 (MQ=255) aaaGCGAGCAGAGCTTATTGAAGGGGAGAAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCt < 1:1677960/64‑1 (MQ=255) aaaGCCAGAAGAGCGTATTGAAGTGGCATAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCt < 1:1648056/64‑1 (MQ=255) aGCCAGAAGGGCGTATTGAACTGGCAGAAGCCCTGCCGCTGACCGACTGCTACGTTGCCCCCtgt < 1:3288796/65‑1 (MQ=255) aGCCAGAAGAGGTTATTGAAGTGGCAGAAGCCCTGCCGCTGACCGACTGCTACGTTGCCCCCtgt < 1:2431713/65‑1 (MQ=255) aGCCAGAAGAGCGTATTGAGGTGGCAGAAGCCCTGCCGCTGACCGACTGCTACGTTGCCCCCtgt < 1:1372423/65‑1 (MQ=255) aGCCAGAAGAGCGTATTGAAGTGGCAGAAGCCCTGCCGCTGACCGACTGCTACGTTGCCCCCtgt < 1:974493/65‑1 (MQ=255) aGCCAGAAGAGCGTATTGAAGTGGCAGAAGCCCTGCCGCTGACCGACTGCTACGTTGCCCCCtgt < 1:2154906/65‑1 (MQ=255) aGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCtgt > 1:4030085/1‑65 (MQ=255) aGCCAAAAGAGCGTATTGAAGGGGCAAAAACCCTGCCGCTGACCGACTGCTACGTTGCCCCCtgt < 1:2329740/65‑1 (MQ=255) agagCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTActg > 1:3980604/1‑65 (MQ=255) agagCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTActg > 1:1509435/1‑65 (MQ=255) agagCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTActg > 1:1781113/1‑65 (MQ=255) agCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTActgcc > 1:259322/1‑65 (MQ=255) aTTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTgcgc > 1:3521836/1‑65 (MQ=255) gAAGTGGCAGAAGCCCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTAt < 1:2013175/65‑1 (MQ=255) gAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTAt > 1:3553720/1‑65 (MQ=255) gAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGCTGCCCCCTCTGTTATTGCCTGCGCTAt > 1:1443098/1‑65 (MQ=255) gTGGCAGAAGGCCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATcaa < 1:3723844/65‑1 (MQ=255) aagcCCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGat < 1:395426/65‑1 (MQ=255) aCCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATAtt > 1:161961/1‑65 (MQ=255) aCCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATAtt > 1:50100/1‑65 (MQ=255) aCCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATAtt > 1:770542/1‑65 (MQ=255) aCCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATAtt > 1:882462/1‑65 (MQ=255) | TAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATT > NC_000913/3017345‑3017469 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |