Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 745,091 | T→A | 15.8% | C283* (TGT→TGA) | ybgK → | putative allophanate hydrolase, subunit 2 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 745,091 | 0 | T | A | 15.8% | 186.3 / 20.5 | 74 | C283* (TGT→TGA) | ybgK | putative allophanate hydrolase, subunit 2 |
Reads supporting (aligned to +/- strand): ref base T (40/22); new base A (7/5); total (47/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.49e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATC > NC_000913/744992‑745189 | gACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAgtgt > 1:739061/1‑100 (MQ=255) gACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAgtgt < 1:747432/100‑1 (MQ=255) gACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAgtgt > 1:957729/1‑100 (MQ=255) ccaccGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTc > 1:2573581/1‑100 (MQ=255) accGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCAc > 1:1399961/1‑100 (MQ=255) ccGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACt > 2:1831/1‑100 (MQ=255) cGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTg < 1:2238835/100‑1 (MQ=255) cGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTg < 2:2263026/100‑1 (MQ=255) ggtggtTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTgg < 2:2423878/100‑1 (MQ=255) ggtggtTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTgg < 2:2564184/100‑1 (MQ=255) gtggtTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGa < 1:1721678/100‑1 (MQ=255) gtTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGAACACTGGacaa < 2:1930893/100‑4 (MQ=255) ttACCCGCGTATTGCTTGTATCATTGAGGCTGATATATACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATATTGTCCAGTGAACACTGGacaaa < 1:536226/100‑5 (MQ=255) aCCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGc < 2:1176801/100‑1 (MQ=255) cccGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTTCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCa > 1:2153780/1‑100 (MQ=255) cccGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGc < 2:2348056/99‑1 (MQ=255) cccGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCa > 1:647951/1‑100 (MQ=255) cccGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCa > 2:460781/1‑100 (MQ=255) cccGCGGAGTGCCCGTATCATTGAGGCTGATATCTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGAACACTGGacaaaatg < 1:229852/100‑8 (MQ=255) ccGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGAACACTGGacaaaatga < 2:1379863/100‑9 (MQ=255) cgcgTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACt < 2:1350686/100‑1 (MQ=255) gcgTATTGCCTTTCTCATTGATGCTGATATGTTCCATCTTGCGTAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGAACTCTGGataaaatga < 1:1539126/98‑9 (MQ=255) aTTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAg < 2:1074374/100‑1 (MQ=255) ttGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGc < 1:2426617/100‑1 (MQ=255) gCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCAc < 2:1604496/100‑1 (MQ=255) gCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCAc > 2:972454/1‑100 (MQ=255) gTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCa > 1:1139653/1‑100 (MQ=255) tATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCaa > 2:709967/1‑100 (MQ=255) aTCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCaa < 2:2573600/99‑1 (MQ=255) tCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGa < 2:1002644/100‑1 (MQ=255) cATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGAt < 2:2150660/100‑1 (MQ=255) ggCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCaa > 1:1040805/1‑100 (MQ=255) cTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACg < 2:555696/100‑1 (MQ=255) cTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACg < 1:2130959/100‑1 (MQ=255) tGATCTGTACCACCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGATCACTGGACGCGGCACTAAAAGCACGGCAATATCAGCAAc > 1:992198/1‑98 (MQ=255) atGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAtt > 1:73518/1‑100 (MQ=255) gTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAtttt > 1:791124/1‑100 (MQ=255) gTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAtttt > 1:340320/1‑100 (MQ=255) gTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAtttt > 1:1655329/1‑100 (MQ=255) gTACCATCTGGCGCAAATTCCGCTCGGTCAACCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAtttt > 2:1489962/1‑100 (MQ=255) aCCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGa < 2:1682314/100‑1 (MQ=255) ngcgcAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTa > 1:427995/2‑100 (MQ=255) ggCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTa > 2:2213341/1‑100 (MQ=255) ggCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTa > 2:2290675/1‑100 (MQ=255) gcgcAAATTCCGCTCGGTAAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCCCTAAAAGCACGGCAAGATCAGCAACGTCATTTTGAACAAtt < 1:820745/98‑1 (MQ=255) cgcAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGc > 1:434209/1‑100 (MQ=255) cgcAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGc > 2:2211027/1‑100 (MQ=255) cgcAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGc > 2:813603/1‑100 (MQ=255) aTTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTggc > 2:260797/1‑100 (MQ=255) aTTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTggc < 1:799244/100‑1 (MQ=255) tCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTggcgg < 2:1130928/100‑1 (MQ=255) ccGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTggcggc > 2:797601/1‑100 (MQ=255) cGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCt > 2:677822/1‑100 (MQ=255) cGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCt > 2:2120594/1‑100 (MQ=255) cGCTCCTTCACCCGATTCAATTTGACCAGGGTTCACTAGCAGAGGAACTAAAAGCACGGCAAGATCAACTACGTTATATTGacaaatcacc > 1:2139850/1‑82 (MQ=255) gCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTg > 2:152877/1‑100 (MQ=255) cTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGc < 2:953159/100‑1 (MQ=255) cGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGcac > 2:2378426/1‑100 (MQ=255) cGGTCAGCCGATTCATTTTGTCCAGTGATCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGcac > 1:300109/1‑100 (MQ=255) ggTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGcaca > 1:472559/1‑100 (MQ=255) ccGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGaaa > 1:182930/1‑100 (MQ=255) gATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAAt > 1:930914/1‑100 (MQ=255) gATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAAt > 2:91118/1‑100 (MQ=255) gATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAAt > 2:91214/1‑100 (MQ=255) gATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAAt > 2:2238723/1‑100 (MQ=255) ttCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTg < 1:1151452/100‑1 (MQ=255) gcctcttCCAGTGAACACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAcc > 2:2225283/7‑100 (MQ=255) gcctcttCCAGTGAACACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAcc > 2:1008529/7‑100 (MQ=255) aTTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAcc > 2:1297670/1‑100 (MQ=255) tttGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGTAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTg > 1:707828/1‑100 (MQ=255) tcttCCAGTGAACACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGa > 2:1583952/4‑100 (MQ=255) cttCCAGTGAACACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGaa > 2:313711/3‑100 (MQ=255) gTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAAc > 1:672248/1‑100 (MQ=255) tCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACg > 1:272132/1‑100 (MQ=255) tCCAGTGAACACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACg > 2:2014156/1‑100 (MQ=255) ccAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGc > 1:2249827/1‑100 (MQ=255) ccAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGc > 2:155146/1‑100 (MQ=255) gtgtTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGa > 1:1844903/1‑100 (MQ=255) tgtTCACTGGAAGAGTCACTAACAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGAt < 1:1807922/100‑1 (MQ=255) gtTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATc > 2:1083112/1‑100 (MQ=255) gtTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATc > 2:64577/1‑100 (MQ=255) | GACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATC > NC_000913/744992‑745189 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |