| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 1,369,793 | G→A | 58.9% | E35E (GAG→GAA) | pspE → | thiosulfate:cyanide sulfurtransferase (rhodanese) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,369,793 | 0 | G | A | 58.9% | 32.8 / 91.8 | 78 | E35E (GAG→GAA) | pspE | thiosulfate:cyanide sulfurtransferase (rhodanese) |
| Reads supporting (aligned to +/- strand): ref base G (21/11); new base A (27/19); total (48/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.38e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01 | |||||||||||
AAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAAT > NC_000913/1369696‑1369892 | aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCa > 2:1985140/1‑100 (MQ=255) aGGCTTACTTGCTCTGGCACTGGTTTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGt < 1:3175970/99‑1 (MQ=255) cTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTCTCAGCAAGAACACGTTCAg > 1:890833/1‑100 (MQ=255) cTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAg > 2:483538/1‑100 (MQ=255) ttacttacTCTGGCACAGGTGTTTTCAATGCACGTTTACGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTACAgg > 1:301746/1‑100 (MQ=255) aCTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAgggg < 2:992495/100‑1 (MQ=255) tGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGcc > 1:341823/1‑100 (MQ=255) ctctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCAt > 2:186116/1‑100 (MQ=255) tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCaa < 1:205337/100‑1 (MQ=255) tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGATCAGGGGGCCATCaa > 2:3346252/1‑100 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAata < 2:2235841/100‑1 (MQ=255) cACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAata > 2:15257/1‑99 (MQ=255) cACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAatat > 1:2206695/1‑100 (MQ=255) gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcccc < 2:3025211/100‑1 (MQ=255) gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcccc < 1:3231133/100‑1 (MQ=255) tgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCt < 2:1629825/100‑1 (MQ=255) ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTga < 2:2297604/100‑1 (MQ=255) ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaag < 1:1555774/100‑1 (MQ=255) ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaag < 2:1610454/100‑1 (MQ=255) cACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaagga > 1:2830708/1‑100 (MQ=255) cTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaaggaag < 2:2365937/100‑1 (MQ=255) cTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaaggaag < 2:2206710/100‑1 (MQ=255) gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg > 2:305962/1‑100 (MQ=255) gCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTg > 2:2211661/1‑100 (MQ=255) cccGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGa > 1:2450171/1‑100 (MQ=255) cctttttCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGa < 1:951908/96‑1 (MQ=255) ccGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGaa < 1:1247412/100‑1 (MQ=255) cGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGaaa > 1:3379353/1‑100 (MQ=255) cGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGaaa > 1:3379616/1‑100 (MQ=255) cGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGaaa > 1:471141/1‑100 (MQ=255) tttCGGCGCTGAACACTGGATCGAGGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAaga < 1:1941615/99‑1 (MQ=255) tttCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAagag < 1:2003246/100‑1 (MQ=255) gccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCa < 1:483536/100‑1 (MQ=255) gccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCa > 2:1816112/1‑100 (MQ=255) ccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAt > 2:1456840/1‑100 (MQ=255) ccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAt > 2:2188932/1‑100 (MQ=255) ccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAt > 1:1605970/1‑100 (MQ=255) cgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAtt > 1:1317809/1‑100 (MQ=255) gcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTg > 2:1666274/1‑100 (MQ=255) cTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGc > 2:1029240/1‑100 (MQ=255) tGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGcc < 2:1298748/100‑1 (MQ=255) gAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCa > 2:1485684/1‑100 (MQ=255) acacTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAcc > 2:1631447/1‑100 (MQ=255) acTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgc > 1:2670074/1‑100 (MQ=255) cTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgcc > 2:72329/1‑100 (MQ=255) ggATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGt > 2:1440995/1‑100 (MQ=255) ggATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGt > 2:2446333/1‑100 (MQ=255) ggATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGt > 2:1081037/1‑100 (MQ=255) ggAACGATGTGAGTGTTCTAGAGCAGTATCAGCAAGAACACGTCCAGCGGACCATCCTTATCCCCCTGAAGGAAGTGAAAGAGCGCACTGCCAccgcg < 1:2556431/98‑2 (MQ=255) gATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGtt > 2:1158159/1‑100 (MQ=255) gATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGtt > 2:902785/1‑100 (MQ=255) aTCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTc < 2:2966485/100‑1 (MQ=255) cGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCg < 2:2128259/100‑1 (MQ=255) cGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCg < 2:258452/100‑1 (MQ=255) gATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCTACCGCCGTTCCgg > 1:3344894/1‑100 (MQ=255) tGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt > 2:1103605/1‑100 (MQ=255) tGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt > 1:2232612/1‑100 (MQ=255) ttCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaa > 2:1835949/1‑100 (MQ=255) tCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaaa > 1:1025560/1‑100 (MQ=255) tCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTGTCGGAt > 2:1340110/1‑97 (MQ=255) cGTGTTCCAGAGCAGTATCAGCAAGAAAACGTTCAGGGGGCCATCAATATTCCACTGAAGGAAGTGAAAGAGCGCATTGCCACCACCGTTCCGGCTaaa > 1:522181/1‑99 (MQ=255) gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg > 2:3340423/1‑100 (MQ=255) ccAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca > 2:3146963/1‑100 (MQ=255) cAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGacac < 2:1207195/100‑1 (MQ=255) cAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGacac < 1:1924855/100‑1 (MQ=255) agagCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACAcc < 1:848892/100‑1 (MQ=255) gagCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAAAAAAACGACACCg > 1:1605532/1‑100 (MQ=255) agCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGt < 2:998199/100‑1 (MQ=255) aGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaa < 2:3379382/100‑1 (MQ=255) aGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaa < 2:3379645/100‑1 (MQ=255) gTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaaa > 1:585159/1‑100 (MQ=255) aTCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgt > 2:778972/1‑100 (MQ=255) aTCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAg > 1:1015222/1‑99 (MQ=255) tCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgtg < 1:863141/100‑1 (MQ=255) cAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgtc > 1:2052585/1‑98 (MQ=255) cAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgtgt > 2:3384111/1‑100 (MQ=255) aGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGAAACCGAGCAAGTGTa > 1:1536329/1‑100 (MQ=255) gCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTAt < 2:2315035/100‑1 (MQ=255) agaaCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGc > 2:595019/1‑100 (MQ=255) gCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAAt > 2:28480/1‑100 (MQ=255) | AAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAAT > NC_000913/1369696‑1369892 |
| 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: ‑ |