Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 788,599 | T→C | 16.3% | intergenic (+46/‑167) | ettA → / → gdhA | ADP/ATP ratio sensor/translation inhibitor/NADP‑dependent glutamate dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 788,599 | 0 | T | C | 16.3% | 72.5 / 8.7 | 37 | intergenic (+46/‑167) | ettA/gdhA | ADP/ATP ratio sensor/translation inhibitor/NADP‑dependent glutamate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (15/16); new base C (2/4); total (17/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01 |
TTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTGAAAAAACCAATAA > NC_002947/788452‑788746 | ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAACATTGGGCCACTTGGGGCCCCAttctt < 1:426234/150‑1 (MQ=255) tcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCATTCTtg < 2:1046060/150‑1 (MQ=255) gaGGTAACTGCACCGAGCACGACGCCGACCCCCA‑ACGCTGCTCGGGTGGTTCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGCTGCACAAGAATGGGGCCCCAAGTGGCCCCATTTCTGTGc < 2:456064/149‑1 (MQ=255) aGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCATTCTTGTGca < 2:86714/150‑2 (MQ=38) gTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTg > 2:217043/1‑150 (MQ=255) tAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTgg < 1:929927/150‑1 (MQ=255) acCGAGTACGAAGCCGACCGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCATTCTtgt > 2:696379/1‑137 (MQ=38) acCGAGTACGAAGCCGACCGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCATTCTtgt < 1:696379/137‑1 (MQ=38) gAAGCCGACCGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATg > 1:395254/1‑150 (MQ=255) aaGCCGACCGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATgg < 1:217043/150‑1 (MQ=255) ccGACCGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcga < 1:49745/150‑1 (MQ=255) cGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCAGAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcgaggcga < 1:447541/150‑1 (MQ=255) gCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa < 2:839766/150‑1 (MQ=255) aagaagCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAgg > 2:1023785/1‑150 (MQ=255) aagaagCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGGTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAgg > 2:999175/1‑150 (MQ=255) t‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGt < 1:1022974/150‑1 (MQ=255) tgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTAc > 2:68704/1‑150 (MQ=255) ccAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGTCTGGCGGACCGCTATAAATGGCGATGCGAAGGTGTATTTTGTATACACTTACATaaaa < 1:278077/150‑1 (MQ=255) ccAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaaa > 1:751633/1‑150 (MQ=255) ggCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTTTGCCTGGCGGACCGCTATAAATGGCCAGGCGAAGCTGTAATTTGTACACGCTTACTTTAAACGCACCAAGTAATTTCa > 2:966933/1‑150 (MQ=255) acaAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGc < 1:884154/72‑1 (MQ=255) acaAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGc > 2:884154/1‑72 (MQ=255) aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGGGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACAATTACATAAAACGCACCAAAAAATTTCaaaaaa > 2:388469/1‑150 (MQ=255) aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa > 1:765942/1‑100 (MQ=255) aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa < 2:765942/100‑1 (MQ=255) ccAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCg < 2:395254/150‑1 (MQ=255) gTAAGCCGAGTTGCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGCGGACGGTAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGccc < 1:68704/150‑1 (MQ=255) tGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGAt > 1:50090/1‑150 (MQ=255) aTGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTg > 2:1033408/1‑150 (MQ=255) ggCACAAAAATGGGGCCACTTGGGGCCCCATTCTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCgg > 2:342995/2‑141 (MQ=38) ggCACAAAAATGGGGCCACTTGGGGCCCCATTCTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCgg < 1:342995/140‑1 (MQ=38) caAGAATGGGGCCCCAAGTGGCCACATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTATGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTAGCCCTATTCCGGTGCGATTGattct > 2:554274/1‑150 (MQ=255) aGACTGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTTTAGTTTGTATACATTTACATAAAACGCCCCAAATAATTCCCAAAAACCGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTg < 1:809400/150‑1 (MQ=255) ggCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACAtt > 2:349962/1‑150 (MQ=255) gcaccATTTTTGTGCCTGGCGGACCGATATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTGAAAAAAc < 2:677214/147‑1 (MQ=255) gCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTGAAAAAAc < 1:1023785/150‑1 (MQ=255) ccccATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTa < 1:17501/110‑1 (MQ=255) ccccATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTa > 2:17501/1‑110 (MQ=255) tttttGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTATATCGTGTATACACTTACATAAAACGCACCAAATAGTTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTGAAAAAACCAATaa < 1:906543/150‑1 (MQ=255) | TTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCT‑GGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTGAAAAAACCAATAA > NC_002947/788452‑788746 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |