Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 276,935 | C→A | 6.1% | T74K (ACG→AAG) | mmuM → | CP4‑6 prophage; S‑methylmethionine:homocysteine methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 276,935 | 0 | C | A | 6.1% | 232.1 / 5.3 | 65 | T74K (ACG→AAG) | mmuM | CP4‑6 prophage; S‑methylmethionine:homocysteine methyltransferase |
Reads supporting (aligned to +/- strand): ref base C (30/31); new base A (2/2); total (32/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01 |
GGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCG > NC_000913/276795‑277053 | ggAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCcggcg > 1:301041/1‑148 (MQ=255) gAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCcggcgg > 2:709077/1‑148 (MQ=255) gtgtAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCcgcg < 1:709077/148‑1 (MQ=255) tAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGa < 4:998581/148‑1 (MQ=255) cTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAg > 1:1029835/1‑146 (MQ=255) cTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAg < 2:1029835/146‑1 (MQ=255) aGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCcggcgg > 1:1052793/1‑128 (MQ=255) aGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCcggcgg < 2:1052793/128‑1 (MQ=255) aGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTg < 2:301041/148‑1 (MQ=255) tgtgGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCcg < 2:507219/113‑1 (MQ=255) tgtgGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCcg > 1:507219/1‑113 (MQ=255) tgtgGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTg < 3:787123/137‑1 (MQ=255) ggtgGTCAGCCAAAGTGCGGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATGTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTg > 4:787123/2‑137 (MQ=255) cAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCcgg > 2:803198/1‑108 (MQ=255) cAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCcgg < 1:803198/108‑1 (MQ=255) cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTcg > 2:1031643/1‑113 (MQ=255) cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTcg < 1:1031643/113‑1 (MQ=255) cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGa > 1:821025/1‑121 (MQ=255) cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGa < 2:821025/121‑1 (MQ=255) cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAgcgc < 1:470104/137‑1 (MQ=255) cAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAgcgc > 2:470104/1‑137 (MQ=255) aaGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCcggcg > 1:329983/1‑104 (MQ=255) aaGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCcggcg < 2:329983/104‑1 (MQ=255) ggTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGc > 1:23232/1‑148 (MQ=255) gTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCa > 3:446218/1‑148 (MQ=255) tAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGa < 1:1053379/118‑1 (MQ=255) tAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGa > 2:1053379/1‑118 (MQ=255) aGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCaaa > 4:489693/1‑148 (MQ=255) aaaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGAAGCCCGCCGGCGTcgcctgata < 1:298112/105‑6 (MQ=38) aaaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGAAGCCCGCCGGCGTcgcctgata > 2:298112/1‑100 (MQ=38) aaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGa > 2:337483/1‑116 (MQ=255) aaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGa < 1:337483/116‑1 (MQ=255) aaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGAAGCCCGCCGGCGTcgcct > 3:700704/1‑99 (MQ=255) aaCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCAATCAGGCGAAGCCCGCCGGCGTcgcct < 4:700704/100‑2 (MQ=255) aCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCcgcgcg < 3:756988/103‑1 (MQ=255) aCCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCcgcgcg > 4:756988/1‑103 (MQ=255) ccGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTgg < 2:23232/148‑1 (MQ=255) tCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGcgcg > 1:600714/1‑148 (MQ=255) cgAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGaa > 4:582228/1‑111 (MQ=255) cgAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGaa < 3:582228/111‑1 (MQ=255) tCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGa > 2:670438/1‑79 (MQ=255) tCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGa < 1:670438/79‑1 (MQ=255) tCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGaa > 3:995996/1‑148 (MQ=255) tGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTgg < 2:930843/83‑1 (MQ=255) tGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTgg > 1:930843/1‑83 (MQ=255) tGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCTTTCGactgcg < 2:591986/103‑6 (MQ=255) tGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCTTTCGactgcg > 1:591986/1‑98 (MQ=255) aTTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCg < 3:719117/117‑1 (MQ=255) aTTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCg > 4:719117/1‑117 (MQ=255) ttACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGAtt > 3:75104/1‑106 (MQ=255) ttACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGAtt < 4:75104/106‑1 (MQ=255) ttACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGt > 3:768182/1‑117 (MQ=255) ttACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGt < 4:768182/117‑1 (MQ=255) ttACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTAt < 4:995996/148‑1 (MQ=255) tactacCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGa < 1:664550/140‑1 (MQ=255) tactacCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGa > 2:664550/1‑140 (MQ=255) actacCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCg > 3:27400/1‑114 (MQ=255) actacCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCg < 4:27400/114‑1 (MQ=255) acCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTaa > 2:369079/1‑127 (MQ=255) acCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTaa < 1:369079/127‑1 (MQ=255) gggggCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAAc < 2:600714/148‑1 (MQ=255) ggggCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAAcc < 3:489693/148‑1 (MQ=255) aTGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCg < 4:446218/148‑1 (MQ=255) cTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGt < 1:403942/85‑1 (MQ=255) cTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGt > 2:403942/1‑85 (MQ=255) | GGAAGCGCGAGGGTGTAACTTAGCCGACAGCCTGTGGTCAGCCAAAGTGCTGGTAGAAAACCCGGAGCTTATCCGCGAAGTGCATCTTGATTACTACCGGGCGGGGGCGCAATGCGCGATCACTGCCAGCTATCAGGCGACGCCGGCGGGCTTCGCCGCGCGAGGTCTGGATGAAGCGCAGTCGAAAGCGCTGATTGGCAAAAGCGTGGAGCTGGCGCGTAAAGCCCGTGAAGCGTATCTGGCAGAGAACCCGCAGGCG > NC_000913/276795‑277053 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |