Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,086,207 | 0 | T | G | 28.0% | 38.7 / 14.4 | 25 | G7G (GGT→GGG) | yqgC | protein YqgC |
Reads supporting (aligned to +/- strand): ref base T (1/17); new base G (7/0); total (8/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.66e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCAGTAACTGTAACCAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGAACGAATCTTCACCAGAACGGTGAGACAGGTT > NC_000913/3086145‑3086267 | tCAGTAACTGTAACCAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGTAt < 1:1182395/65‑1 (MQ=255) tCAGTAACTGTAACCAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGTAt < 1:101448/65‑1 (MQ=255) aGTAACTGTAACCAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGTATGc < 1:2159760/65‑1 (MQ=255) aCTGTAACCAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGGATGCAAAg > 1:1252991/1‑65 (MQ=255) tGTAACCAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGTATGCAAAGCa < 1:614447/65‑1 (MQ=255) tGTAACCAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGGGTGCAAAGCa > 1:762442/1‑65 (MQ=255) tGTAACCAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGGATGCCAAGCa > 1:772650/1‑65 (MQ=255) tGTAACCAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGGATGCAAAGCa > 1:709156/1‑65 (MQ=255) aCCAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGGATGCCAAGCagaag > 1:1513390/1‑63 (MQ=255) ccAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGGATGGAAAGCAGAGAt > 1:324579/1‑65 (MQ=255) ccAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGGATGCAAAGGAGAGAt > 1:2022386/1‑65 (MQ=255) ccAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGGATGCAAAGCCGAGAt > 1:56005/1‑65 (MQ=255) ccAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGGATGCAAAGCAGAGAt > 1:2073769/1‑65 (MQ=255) ccAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGGATGCAAAGCAGAGAt > 1:327085/1‑65 (MQ=255) ccAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGGAAGCAAAGCagaggt > 1:1230287/1‑63 (MQ=255) tCCTGAGAAGCGCCGATATGCGTATAACATCGGCAGGTATGCAAAGCAGAGATGCAGAGTGCggg < 1:2947910/65‑1 (MQ=255) ccTGAGAAGTGCCGAGATGGGTATAACATCGGCAGGGGTGCAAAGCAGAGATGCAGAGGGCgggg > 1:810816/1‑65 (MQ=255) ccTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGTATGCAAAGCAGAGATGCAGAGTGCgggg > 1:1236049/1‑65 (MQ=255) ccTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGGATGGAAAGCAGAGATGCAGAGGGCgggg > 1:901888/1‑65 (MQ=255) cTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGa < 1:954250/65‑1 (MQ=255) gaAGCGCCGAGATGGGTATAACATCGGCAGGTATGCAAAGCAAAGACGCCGAGCGGGGGGAACGa > 1:2941598/1‑65 (MQ=255) gaAGCGCCGAGATGGGTATAACATCGGCAGGGATGCAAAGCAGAGATGCAGAGGGCGGGGGACGa > 1:3000901/1‑65 (MQ=255) gaAGCGCCGAGATGGGTATAACATCGGCAGGGATGCAAAGCAGAGATGCAGAGGGCGGGGAAAGa > 1:2738119/1‑65 (MQ=255) aGCGCCGAGATGGGTATAACATCGGCAGGTATTCAAAGCAGAGATGCAGAGTGCGGGGAACGAAt < 1:247759/65‑1 (MQ=255) aGCGCCGAGATGGGTATAACATCGGCAGGGATGCAAAACAGAGATGCAGAGGGCGGGGAACGAAt > 1:729768/1‑65 (MQ=255) agaTGGGTATAACATCGGCAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGAACGAATCTTCAcc < 1:1709913/65‑1 (MQ=255) gaTGGGTATAACATCGGCAGGTATGCAAAGCAGAGAAGCAGAGGGCGGGGAACGCATCTTCCCCa > 1:1985045/1‑65 (MQ=255) aTGGGTATAACATCGGCAGGGATGCAAAACAGAGATGCAGAGTGCGGGGAACGAATCTTCACCAg > 1:1464136/1‑65 (MQ=255) tGGGTATAACATCGGCAGGGAAGCAAAGCAGGGGTGCAGAGAGCGGGGGACGAAAATTCACCAGa > 1:2013064/1‑65 (MQ=255) gTATAACATCGGCAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGAACGAATCTTCACCAGAACg < 1:3048465/65‑1 (MQ=255) gTATAACATCGGCAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGAACGAATCTTCACCAGAACg < 1:2069435/65‑1 (MQ=255) tataACATCGGCAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGAACGAATCTTCACCAGAACgg < 1:3203705/65‑1 (MQ=255) taACATCGGCAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGAACGAATCTTCACCAGAACGGTg < 1:2402053/65‑1 (MQ=255) cATCGGCAGGGATGCAAAGCAGAGATGCAGAGTGGGGGGAACGAATCTTCACCAGAACGGGgaga > 1:317252/1‑65 (MQ=255) aTCGGCAGGGATGCAAAGCAGAGATACAGAGTGCGGGGGACGAATCTTCACCAGAACGGGGAGAc > 1:212578/1‑65 (MQ=255) tCGGCAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGAACGAATCTTCACCAGAACGGTGAGACa < 1:2912830/65‑1 (MQ=255) cGGCAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGAACGAATCTTCACCAGAACGGTGAGACaa < 1:1241760/65‑2 (MQ=255) cGGCAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGAACGAATCTTCACCAGAACGGTGAGACAg < 1:1886753/65‑1 (MQ=255) gCAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGAACGAATCTTCACCAGAACGGTGAGACAGGt < 1:1302715/65‑1 (MQ=255) cAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGAACGAATCTTCACCAGAACGGTGAGACAGGtt < 1:384428/65‑1 (MQ=255) | TCAGTAACTGTAACCAGCTCTTGAATCCTGAGAAGCGCCGAGATGGGTATAACATCGGCAGGTATGCAAAGCAGAGATGCAGAGTGCGGGGAACGAATCTTCACCAGAACGGTGAGACAGGTT > NC_000913/3086145‑3086267 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |