Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,127,250 | C→T | 7.9% | intergenic (+26/+35) | rpmE → / ← yiiX | 50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,127,250 | 0 | C | T | 7.9% | 127.9 / 3.7 | 52 | intergenic (+26/+35) | rpmE/yiiX | 50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein |
Reads supporting (aligned to +/- strand): ref base C (25/21); new base T (2/2); total (28/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01 |
CTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTC > NC_000913/4127105‑4127398 | cTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCCAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCACCAAGCGTTTCAACATCCCGGGCGCCAAATAAGACCGGATTTCCGAAAAAAGCcacgcacg > 2:139475/1‑149 (MQ=255) aCCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑gcgc < 2:403086/149‑1 (MQ=255) ccTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATACCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCt > 1:8237/1‑149 (MQ=255) cGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGTGCG‑GCGCttttttc < 1:726794/149‑2 (MQ=255) cGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGTGCG‑GCGCttttttc < 2:450798/149‑2 (MQ=255) gtgtGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTtg < 2:802071/149‑1 (MQ=255) gtgtGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTtg < 1:274772/149‑1 (MQ=255) gtgtGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGATACCGGTGGCCGTTTTGTCCGCTTCCACAAGTGTTTCAACATGCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTtg < 2:523646/149‑1 (MQ=255) gtgCAGCAAGTGCCACCCGTTCTTAACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCTACAAGCGTTTCAACATGCCGGGCAGCAAAGAAGACCGGATTTTCGAAAAAAGCGGCGCACG‑GCGCTTTTTTtgtt > 1:590763/1‑148 (MQ=255) cagcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCCCG‑CCGCCTTTTTTGTGccc > 1:337706/1‑148 (MQ=255) agcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCTGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTg > 1:75021/1‑149 (MQ=255) caAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGCTGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACCATCGTTTCACCATCCCTGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGaa < 1:789321/149‑1 (MQ=255) ccACCCGTTCTTCACGGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTgg < 1:173771/149‑1 (MQ=255) ttcttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGCTTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTAtt > 2:38099/1‑149 (MQ=255) aaCAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAg < 1:481682/149‑1 (MQ=255) aaCAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTAAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTCTGTTTATGCCTGAATTTTGGTCGTATTACCACGCAAACAg > 1:618532/1‑149 (MQ=255) aGCGTGATGTTGCTGCCGGTGGCCGAGTTGACCTCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTACGAAATAAGCGCCACACG‑GCGCTTTTTTTGTGCCTGATTTTTGGTCGTATTACCACGAAAACAGCgg < 1:488806/149‑1 (MQ=255) gTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCtt > 1:406498/1‑149 (MQ=255) gCTACCGGTGGCCGTGTTGCCCGCTTCAACAAGCTCCTCGACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTgg < 1:149434/149‑1 (MQ=255) gCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACGGCGGCCATTCTTTgg > 1:138511/1‑149 (MQ=255) cGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATCAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGAcc > 1:354042/1‑149 (MQ=255) gTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgt > 2:49286/1‑149 (MQ=255) tGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGAACG‑TCTTTTGTTGTGGGCATCCCCTTTAGGCTTATTACCACGAAAACAGCGTCCATTCTTTGGCGACCgtg < 1:190523/149‑1 (MQ=255) gtTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGtt < 1:7209/149‑1 (MQ=255) ccGCTTCAACAAGCGTTTCAACATCCCGGGCAACAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTATGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGgc < 2:447410/149‑1 (MQ=255) cgtttCAACAAGCGTTTCAACATCCAGGCCAGCAAATCAGACCGTACGGTTAAAAAAAATG‑CTAACGTGCGCTTAAGTTGTGCCTGAACTTTGTTCGTATAACCCCGGCAACAGCGGCCATTCTTTGGCGGCCGGGGTGGGCGGTGgcg > 2:369421/4‑149 (MQ=255) caaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTCGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATcccga > 2:596389/1‑145 (MQ=255) tttCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCt > 1:525150/1‑149 (MQ=255) ttCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCtt > 2:181225/1‑149 (MQ=255) aaCATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTGGTGGCGCATCGAAAACGGCCtgtg > 2:675994/1‑149 (MQ=255) aCATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCATGAAGTTTGGTAGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGGGGTGTGGCGCATCGAAAACGGCTTGTgg > 1:789137/1‑149 (MQ=255) tCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCgc > 1:164767/1‑148 (MQ=255) tCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTGGTGGCga > 1:760486/1‑149 (MQ=255) gggCAGCAAATAAGACCGGATTTACGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACGACGCAAACAGCGGCCATTCTTTGGGGACCGTGGTGAGTTGTGGCTCATCGAAAACGGCTTGTGGCGCGacc > 2:2200/1‑149 (MQ=255) ccagccAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGaccac < 2:590763/148‑1 (MQ=255) agcaAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCg < 1:287837/149‑1 (MQ=255) caAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGtt < 2:774509/149‑1 (MQ=255) aactAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGttt < 1:129456/146‑1 (MQ=255) aaaTAAGTCCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGttt < 2:381981/149‑1 (MQ=255) aaaTAAGCCCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGttt < 1:264665/149‑1 (MQ=255) aTAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGGGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTGTcc > 1:463255/1‑149 (MQ=255) tAAGACCGGAGTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTCGCGCCTGAATTTTGGTCGTCTTACCACGAAACCAGCGCCCATTCTTCGGCGACCGTGGTGAGTTGTGGCGCATCGAAATCGTCTTGTGGCGAGACCACCGTTTCCt < 2:388464/149‑1 (MQ=255) gtcgggATGTCCGAAGAAAGCGCCGCACG‑GCGCTTTTTTGGTGCCTGAATGTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTCGTGAGTTGTGGCGCATCGAAAACGGCCTGTGGCGAGATCACCGTGTCCTCaa < 2:16669/145‑1 (MQ=255) ccGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGGTGTGGCGCATCGAAAACGGCTTGTGGCGAGAACACCGGTTCCTCAAGc > 2:13696/1‑149 (MQ=255) tttCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGCGTTGTGGCGCATCGAAAACGGCTTCTGGCGAGACCACCGCTTCCTCCAGCGGAat > 2:15265/1‑149 (MQ=255) ccGAAAACAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCTTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATAtt < 1:480675/149‑1 (MQ=255) ccGAAAAAAGCGCCGCACG‑GCGCTTGTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCGCCGTTTCCTCAAGCGGAATAtt > 1:401131/1‑149 (MQ=255) gAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCt > 2:152404/1‑149 (MQ=255) gAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCt > 2:654618/1‑149 (MQ=255) gAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGAACACGGTTTCTTCAAGCGGAATATTCt > 2:152718/1‑149 (MQ=255) aaaaaaGCGCCGTGCG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCtt > 1:243324/1‑149 (MQ=255) aaaaGCGCCGTGCG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGACCCCTCAAGAGGAATATTCTTGc > 1:733110/1‑149 (MQ=255) aaaaGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTGTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGc < 1:750088/149‑1 (MQ=255) cACG‑GCGCTTTTTTTTTTTTTTTATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGACTATTCTTGCCGTGACGttt > 2:63948/1‑148 (MQ=255) | CTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACG‑GCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTC > NC_000913/4127105‑4127398 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |