Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,681,932 | T→C | P33P (CCA→CCG) | ydgC ← | GlpM family inner membrane protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,681,932 | 0 | T | C | 96.3% | 274.8 / 0.2 | 81 | P33P (CCA→CCG) | ydgC | GlpM family inner membrane protein |
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (44/34); total (45/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.82e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01 |
GACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAG > NC_000913/1681837‑1682024 | gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGaa > 1:609303/1‑100 (MQ=255) gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa < 2:620533/100‑1 (MQ=255) ccACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATc > 1:1673465/1‑100 (MQ=255) aCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCg > 1:916657/1‑100 (MQ=255) cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg > 1:225984/1‑100 (MQ=255) cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg > 2:64778/1‑100 (MQ=255) cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg > 1:153766/1‑100 (MQ=255) tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc > 2:342183/1‑100 (MQ=255) tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc < 1:1404235/100‑1 (MQ=255) tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc > 2:1108161/1‑100 (MQ=255) aaaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCg > 2:35027/1‑100 (MQ=255) aaaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCg > 1:117740/1‑100 (MQ=255) ataaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGGAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCg < 1:1842305/97‑1 (MQ=255) aaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGa < 1:50406/100‑1 (MQ=255) aaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGa > 1:838309/1‑100 (MQ=255) aaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGata < 1:796632/100‑1 (MQ=255) atgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGatat < 2:759520/100‑1 (MQ=255) gatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATataa > 2:1173432/1‑100 (MQ=255) gcgcGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAtttt > 2:1442957/1‑100 (MQ=255) cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt > 1:1901708/1‑100 (MQ=255) cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt < 2:1673452/100‑1 (MQ=255) gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc > 1:1617342/1‑100 (MQ=255) gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc < 1:1905583/100‑1 (MQ=255) cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCg < 2:153765/100‑1 (MQ=255) gTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTAGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGt < 2:726105/100‑1 (MQ=255) tAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGtt > 1:76641/1‑100 (MQ=255) aaGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCACGGCGATATAATAATTTTTCGtt > 2:1921389/1‑99 (MQ=255) aGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGtttt > 2:107983/1‑100 (MQ=255) aGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGtttt < 1:723637/100‑1 (MQ=255) ggCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttttt > 1:1977060/1‑100 (MQ=255) ggCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttttt < 1:1682896/100‑1 (MQ=255) gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg < 2:354245/100‑1 (MQ=255) gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg > 1:901054/1‑100 (MQ=255) catcaacctggcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc < 2:1897094/90‑1 (MQ=255) cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc > 2:805650/1‑100 (MQ=255) cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc < 1:1173438/100‑1 (MQ=255) ttCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCt < 2:1665989/100‑1 (MQ=255) cAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaa < 2:1845759/100‑1 (MQ=255) aTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaaa > 1:766587/1‑100 (MQ=255) tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac > 1:75152/1‑100 (MQ=255) gCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac < 2:1078649/99‑1 (MQ=255) ccGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAacac > 1:710384/1‑100 (MQ=255) cgatctCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAcc < 2:1944069/95‑1 (MQ=255) gcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCa > 2:1599057/1‑100 (MQ=255) ttCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATc < 1:923601/100‑1 (MQ=255) tCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCa < 2:1995109/100‑1 (MQ=255) cGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGACATAATAATTTTTCGTTTTTGCTAAAACACCagatcg > 1:1944087/1‑95 (MQ=255) cTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa < 1:1722170/100‑1 (MQ=255) cTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa > 1:287135/1‑100 (MQ=255) ggCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGc < 1:1093082/100‑1 (MQ=255) ggCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGc < 1:1877479/100‑1 (MQ=255) ggCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGc > 1:265463/1‑100 (MQ=255) gCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGc < 1:805653/99‑1 (MQ=255) caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc > 1:1205148/1‑100 (MQ=255) caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc > 2:1158162/1‑100 (MQ=255) caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc > 1:235159/1‑100 (MQ=255) acaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTa < 1:64778/100‑1 (MQ=255) acaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGAAATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTa > 1:142733/1‑100 (MQ=255) aTGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCa > 2:1052792/1‑100 (MQ=255) tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc > 2:279982/1‑100 (MQ=255) cgcgATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCg > 2:1333516/1‑100 (MQ=255) cgATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAg > 1:1120545/1‑100 (MQ=255) aaGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGcgcg < 2:142732/100‑1 (MQ=255) gcgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGc > 1:1987666/1‑100 (MQ=255) gcgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGc < 1:1108166/100‑1 (MQ=255) cgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCt > 2:1053355/1‑100 (MQ=255) gcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCtt < 1:35027/100‑1 (MQ=255) cAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCttt < 2:1981440/100‑1 (MQ=255) aaaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTg < 2:997172/100‑1 (MQ=255) aaaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTg < 2:766584/100‑1 (MQ=255) aaaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGCTTTTGCTAAAACACCAAGCAGCTGCACTACCAGCGCACCGAGCGCGGCTTTg < 1:139664/100‑1 (MQ=255) aaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGa > 2:441930/1‑100 (MQ=255) aGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTTTTTATACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAt > 2:1763721/1‑100 (MQ=255) aGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAt < 2:648557/100‑1 (MQ=255) ggTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAtt > 1:1320388/1‑100 (MQ=255) tCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAc < 2:225983/100‑1 (MQ=255) tCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAc > 2:1204332/1‑100 (MQ=255) cGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTGTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAcc < 1:1514036/100‑1 (MQ=255) gggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCa > 1:1828588/1‑100 (MQ=255) gggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCa > 1:1104578/1‑100 (MQ=255) ggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAg > 1:156913/1‑100 (MQ=255) | GACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAG > NC_000913/1681837‑1682024 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |