Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,691,173 | A→T | 10.4% | noncoding (168/184 nt) | glmY ← | sRNA activator of glmS mRNA, glmZ processing antagonist |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,691,173 | 0 | A | T | 10.4% | 190.9 / 7.8 | 59 | noncoding (168/184 nt) | glmY | sRNA activator of glmS mRNA, glmZ processing antagonist |
Reads supporting (aligned to +/- strand): ref base A (27/24); new base T (4/2); total (31/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.78e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATA > NC_000913/2691029‑2691321 | gcagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCg > 4:557651/1‑149 (MQ=255) aaTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATT‑CCCGGCtt < 2:337275/149‑1 (MQ=255) aaTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCAGGCtt > 4:557942/1‑149 (MQ=255) gCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTa > 3:349659/1‑149 (MQ=255) cAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATT‑CCCGGCTTTGttgt < 4:411570/149‑3 (MQ=255) cAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTAt < 4:207968/149‑1 (MQ=255) ccAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAg < 1:215638/149‑1 (MQ=255) aTTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGt < 3:557651/149‑1 (MQ=255) aTTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGt < 3:126264/149‑1 (MQ=255) gTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTAcc > 3:453869/1‑117 (MQ=255) gTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTAcc < 4:453869/117‑1 (MQ=255) tGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTc > 3:71623/1‑149 (MQ=255) gTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACg > 4:430921/1‑149 (MQ=255) tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGt < 3:302901/149‑1 (MQ=255) gCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGcc < 4:264407/149‑1 (MQ=255) ggggAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGGCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCg < 3:557942/149‑1 (MQ=255) aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt > 1:94514/1‑149 (MQ=255) aaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTg > 3:292066/1‑149 (MQ=255) ggCCAGCGTTTCAAGGTGTTACTCTCGACAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGatat > 4:311095/1‑149 (MQ=255) gCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTAccc > 3:423707/1‑91 (MQ=255) gCCAGCGTTTCAAGGTGTTACTCTCGTAAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTAccc < 4:423707/91‑1 (MQ=255) cAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGc > 3:135662/1‑149 (MQ=255) gCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGata > 1:362776/1‑143 (MQ=255) gCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGata < 2:362776/143‑1 (MQ=255) cAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCTTTg > 2:405039/1‑88 (MQ=38) cAAGGTGTTACTCTCGTCAGACGCGAATAGCCTAATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCTTTg < 1:405039/88‑1 (MQ=38) ggTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCtt < 3:325236/83‑1 (MQ=39) ggTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCtt > 4:325236/1‑83 (MQ=38) tACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGaa < 3:311095/149‑1 (MQ=255) tACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGaa < 1:177478/149‑1 (MQ=255) tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATg < 3:294539/82‑1 (MQ=255) tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATg > 4:294539/1‑82 (MQ=255) gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGAc > 4:499278/1‑149 (MQ=255) tAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGAt < 3:192892/102‑1 (MQ=255) tAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGAt > 4:192892/1‑102 (MQ=255) tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt > 2:151615/1‑149 (MQ=255) gATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCgg > 2:470521/1‑50 (MQ=25) aTGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGc < 1:391441/149‑1 (MQ=255) aaCCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGt < 2:543153/120‑1 (MQ=255) aaCCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGt > 1:543153/1‑120 (MQ=255) ggAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAttt < 3:499278/149‑1 (MQ=255) gTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATg < 3:430921/149‑1 (MQ=255) gATACAACCAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt < 3:456919/77‑1 (MQ=255) gATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt > 4:456919/1‑77 (MQ=255) aTACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCAc > 4:61345/1‑149 (MQ=255) aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGAtata > 3:485480/1‑78 (MQ=255) aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGAtata < 4:485480/78‑1 (MQ=255) aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTc < 4:487871/129‑1 (MQ=255) aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTc > 3:487871/1‑129 (MQ=255) aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTcc < 4:71623/149‑1 (MQ=255) acaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca > 3:144327/1‑106 (MQ=255) acaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca < 4:144327/106‑1 (MQ=255) acaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTccc > 4:44297/1‑149 (MQ=255) aaGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa > 3:105739/1‑149 (MQ=255) aGCCGGG‑TAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg > 2:566578/1‑149 (MQ=255) ccGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAggg > 3:203305/1‑149 (MQ=255) ccGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAggg > 1:346353/1‑149 (MQ=255) ggg‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGc > 2:76443/1‑149 (MQ=255) ggg‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGc < 4:135662/149‑1 (MQ=255) aCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGa > 3:621524/1‑72 (MQ=255) aCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGa < 4:621524/72‑1 (MQ=255) aCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATa < 2:94514/149‑1 (MQ=255) | GCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATA > NC_000913/2691029‑2691321 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |