Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,899,378 | T→C | T6A (ACC→GCC) | yieK ← | putative glucosamine‑6‑phosphate deaminase YieK |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,899,378 | 0 | T | C | 95.8% | 73.5 / ‑3.7 | 24 | T6A (ACC→GCC) | yieK | putative glucosamine‑6‑phosphate deaminase YieK |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (10/13); total (11/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.58e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01 |
GTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCC > NC_000913/3899301‑3899451 | gTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAAt > 1:773828/1‑81 (MQ=255) gTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAAt < 1:773794/81‑1 (MQ=255) acacGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGAt < 1:773795/80‑1 (MQ=255) acacGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGAt < 1:773796/80‑1 (MQ=255) gacgCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAAtt < 2:773797/81‑1 (MQ=255) acgCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAAttt > 2:773823/1‑81 (MQ=255) cgcgTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTc < 1:773798/81‑1 (MQ=255) tCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCt < 1:773799/81‑1 (MQ=255) cTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCt > 1:773825/1‑80 (MQ=255) cATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGcc < 2:773800/81‑1 (MQ=255) taACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCttat < 2:773801/81‑1 (MQ=255) taACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCttat < 2:773802/81‑1 (MQ=255) cAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTtattttattt > 2:773824/1‑80 (MQ=255) gATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTaaa > 1:773822/1‑81 (MQ=255) gATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTaaa > 2:773829/1‑81 (MQ=255) cgcgACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTg > 2:773830/1‑80 (MQ=255) gACACGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTAc > 2:773851/1‑80 (MQ=255) acacGGCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACa > 1:773827/1‑80 (MQ=255) gCTCATTTCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCt < 2:773803/81‑1 (MQ=255) ttCCTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAg > 1:773840/1‑81 (MQ=255) ccTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGc > 2:773835/1‑81 (MQ=255) cTGGTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATGTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGc < 1:773804/80‑1 (MQ=255) gTAATCTTCGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGt < 1:773805/81‑1 (MQ=255) cGGCAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGcc < 2:773806/77‑1 (MQ=255) | GTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCC > NC_000913/3899301‑3899451 |
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: ‑ |