Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,782,908 | 0 | T | G | 41.7% | 36.9 / 15.3 | 51 | T252P (ACC→CCC) | gpsA | glycerol‑3‑phosphate dehydrogenase (NAD+) |
Reads supporting (aligned to +/- strand): ref base T (14/14); new base G (20/0); total (35/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.78e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCC > NC_000913/3782762‑3783056 | gAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGAATGGGTGTcgg‑taacggtaa > 2:1071577/1‑149 (MQ=255) aGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGGCTGGTTGGCGG‑TTCAGGGaa > 2:1053547/1‑147 (MQ=255) cGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑TACAGGTAAGCACCa < 2:571720/151‑1 (MQ=255) tCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGc < 2:502022/150‑1 (MQ=255) cgcgGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGccc < 1:765115/151‑1 (MQ=255) cgGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCGCCGCTGGTTTTCGG‑GACAGGGAAACACCAAAACGCCCAAcccc > 1:730402/1‑150 (MQ=255) ttcttTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGCTACGCGACTGGCTGTCTG‑GAAAGGTAAGCACCAAAACGCCCAACCCCGcc > 2:362715/1‑148 (MQ=255) cttTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGGTACCCGACTGGTTGTGGGCTAC‑GGGAAACAACAAAAAGCCCAGCCCCGCCCTGcc > 1:64365/1‑151 (MQ=255) tttCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTTGTTGTCGG‑TACAGGTGAGCACCAAATAGCCCACCCCCCCCCCGccc > 1:49802/1‑151 (MQ=255) ttCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGc < 1:618357/148‑1 (MQ=255) cGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTCCGCGCCTGGTTGTTGG‑GACCAGGAAGGACCAAAATGCCAAACCCCCCCCCGcccct > 1:1173608/1‑148 (MQ=255) gTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACCGGGTGTCGG‑GACAGGTAAACACCAAAACGCCAAGCCCCCCCCTGCCCCTa > 1:1189300/1‑150 (MQ=255) gTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTAACCGACTGGGTGTCGG‑GACAGGGAAGCCCCAGATCGCCCAACCCCCCCCTGcccc > 2:653619/1‑147 (MQ=255) gCGGTAGCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTg < 2:213343/151‑1 (MQ=255) ggTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACCGGTTGTCGG‑TACAGGGAAACACCAAAAAGCCCAACCCCCCCACGCCCCTaaaagggg > 1:1065624/1‑145 (MQ=255) gCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGGTGTCGG‑TACAGGGAAGCACCAGATTGCCCAACCCCCCCCTGCCCATAAAAGGGGCAggg > 2:59154/1‑151 (MQ=255) ccTTCCACCCCCTGACCAATCTTCTCCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAgg < 1:410082/149‑1 (MQ=255) ccTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGACTGGTTGTCGG‑GACAGGGAAGCACCAAAACGCCCAACCCCGCCATGCCCCTAAAGGGGGGAggg > 2:865057/1‑150 (MQ=255) tCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑GACAGGGAAGCACCAGATCGCCAAGCCCCCCCATGCCCATAAAGGTGGCAGGGTCgg > 2:94437/1‑151 (MQ=255) ccaccaCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGACTGGGTTTCGG‑GACAGGGAAACACCAAAACTCCAAGCCCCCCCTACCCCATAAAGGGGGCTgggcggg > 2:958844/1‑146 (MQ=255) ccAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTAAGCGGCTGGTTTTCGG‑TAAAGGAAAACACCAAAGAGCCAAGCCCCCCCCTGCCCCCTAAAGTGGCCGGGGCGGGGCCCaaagcc > 1:114149/1‑145 (MQ=255) gcgcgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTAAGCGAATGGTTGTCGG‑GACAGGGAAGCACCAAATAGCCAAGCCCCGCCATGCCCCTAAAAGGGGCAGGGGCGGCACCCCGCGCCCCACCAAGAcg > 2:1007624/1‑148 (MQ=255) cgcgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACTCCGCTGGTTGTCGG‑GACAAGTAAGCAACAGAAAGCCAAACCCCGCCCAGCCCATAAAAGGGGCAGGGTCGGCCCCCAACCCCCCCCCAAGa > 2:827755/1‑145 (MQ=255) gcgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑TAAAAGGAAGAACCAAAACGCCCAGCCCCCCCACGCCCCTAAAAGGGGGAGGGGCGGGACCCCGGGCCCCAccca > 1:886834/1‑140 (MQ=255) cgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTTCCGACTGGTTGTTGG‑TACAAGTAAGCACCACAACGCCAAACCCCGCCAAGCCCCCAAAAGAGGCAGGGGCGGCACCCCGCGCCGCACCCAAAcgc > 1:485637/1‑146 (MQ=255) aTGCCCTGACCGAGCATCATGCCAAAACGGCGGTTATGCGGCTGGGTGTTCG‑GACCGGTAAACCCCAAAAAGCCAAAGCCCCCCCTGCCCCTAAAAGGGGCAGGGGCGGGACCCCGGGCCGCCCCAAAACCCGCCATTTCAGCCAcccccc > 1:416758/1‑149 (MQ=255) gCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGGCCGGGTGTCGG‑GACAGGGAAGCACCAAAACGCCAAACCCCGCCCTGCCCATAAAAGGGGCAGGGGCGGCACCCCGCGCCCCCCCAAAACGCGGCATTTcagcccaccccccg > 1:681936/1‑142 (MQ=255) cccTGACAGAGCATCATCCCAAAACGGCGGTCACGCGACTGGTTGCCGG‑TACAGGAAAGCACCAGACCCCCACCCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCAc < 1:194571/148‑1 (MQ=255) ccTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑GGCAGGGAAGAACCAAAACGCCCAACCCCCCCCCGCCCATAAAGGGGGCAGGGGCGGCACCCCGCGCCGCACCCAGACGCGCCATTTCAGCCAGccccccggg > 1:867712/1‑145 (MQ=255) ccTGACCGAGCATCATGCCAAAACGGCGGTTACACGGCTTGTTGTCGG‑GACAGGGAAGCACCAAAGCCCCCACCCCCCCCAGCCCCATAAAAGTGGCAGGGGCGGGACCCCCCCCCCCAACAAAAAGCGACCCTTTAGCCAGCCCAccg > 2:131880/1‑147 (MQ=255) gACCGAGCATCATGCCAAAACGGCGGTTTTGCGACCGGGTGTGGG‑GACAGGGAAGCACCAAAAAGGCAAACCCCCCCCCGCCCCTTAAAGGGGCAGGGGCGGCACCCCGCGCCGCACCCAGACGCGCCATTTCACCCAGccccccgggg > 1:683462/1‑142 (MQ=255) gACCGAGCATCATGCCAAAACGGCGGTTACGCGCCTGGGTGTCGG‑GACAGGGAAGCACCACATCGCCAAACCCCCCCCTGCCCCTAAAAGTGGCAGGGTCGGCACCCCGCGCCGCACCCAGACGCGACATTTCAGCCAGCCCACGGGTg > 2:647588/1‑149 (MQ=255) gACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGa < 2:899981/150‑1 (MQ=255) ccGAGCATCATGCCAAAACGGCGGTTACGCGACTGGGTTTCGG‑GACAGGGAAGCACCAAAACGCCAAGCCCCGCCCTGCCCCTAAAAGGGGCAGGGGCGGCACCCCGCGCCGCACCCAAACGCGACATTTCAGCCAGCCCCCggggggtc > 1:1042050/1‑145 (MQ=255) cGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑TACAGGGAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAAGGGGCAGGGGCGGCACCCAGCGCCGCACCCAGACGCGACATTTCAGCCAGCCCACGGGGGATc > 2:1002030/1‑149 (MQ=255) aGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑GACAGGGAAACACCAAATCGCCAAGCCCCGCCATGCCCATAAAGGGGGCAGGGGCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGGGATCAgc > 2:905164/1‑150 (MQ=255) atcatGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGcc < 1:434085/150‑1 (MQ=255) caaaaaCGGCGGTTCCGCGACTGGTTGTCGG‑TACAGGTAACCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTAcgcg < 1:289088/149‑1 (MQ=255) aaCGGCGGTTACGCGACTGTTTGACGA‑TACAGGAAACCACCAGACCCCCAGGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTAcgc < 1:1071515/146‑1 (MQ=255) ttgcgcgACTGGTTGTCGG‑TACAGGTAACCCCCAGAGCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCAcc < 2:557679/146‑1 (MQ=255) aCGCGACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCaaa > 1:887333/1‑150 (MQ=255) aCGCGACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCaaa > 2:653998/1‑150 (MQ=255) cgcgACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGAACGCGCATTCGCACCaaa > 2:597254/1‑149 (MQ=255) gcgACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAAcc < 2:1023602/151‑1 (MQ=255) ggcgggTTGGCGG‑TGCAGGTAAGCACCAGATCGCGAGGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATGTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCg < 1:248316/146‑1 (MQ=255) gACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCg > 2:662707/1‑150 (MQ=255) gTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCg > 1:74477/1‑151 (MQ=255) gTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCg > 1:457211/1‑151 (MQ=255) cGG‑TACAGGTAAGCCCCAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGa < 2:886879/151‑1 (MQ=255) tACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACAt > 1:292725/1‑151 (MQ=255) tACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACAt > 2:500018/1‑151 (MQ=255) ggTAAGCACCAGATCCCCAACCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATcccc < 1:316978/151‑1 (MQ=255) | GAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGG‑TACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCC > NC_000913/3782762‑3783056 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |