Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,227,550 | G→A | 8.6% | intergenic (+39/‑181) | metH → / → yjbB | homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,227,550 | 0 | G | A | 8.6% | 175.0 / 4.8 | 58 | intergenic (+39/‑181) | metH/yjbB | homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter |
Reads supporting (aligned to +/- strand): ref base G (27/26); new base A (2/3); total (29/29) | |||||||||||
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.87e-01 |
TCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAAT > NC_000913/4227412‑4227691 | tCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg < 4:8505/148‑1 (MQ=255) gcgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCt < 2:354151/149‑1 (MQ=255) aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGt < 3:98240/148‑1 (MQ=255) aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtt < 2:313289/149‑2 (MQ=255) gTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt < 1:143590/149‑1 (MQ=255) gTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt < 2:184825/149‑1 (MQ=255) tGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta < 4:258279/148‑1 (MQ=255) ttATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata > 3:162866/1‑148 (MQ=255) gTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTta < 3:230721/148‑1 (MQ=255) aaaGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTtata > 3:112112/1‑148 (MQ=255) ggTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg < 1:379909/149‑1 (MQ=255) aTGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa < 1:13168/149‑1 (MQ=255) aTGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACAGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa < 3:330453/148‑1 (MQ=255) gTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCt > 1:107767/1‑149 (MQ=255) tggcACCGAATCTGGGGTATGACGCGGACTGATTCGCAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaa < 3:276643/148‑1 (MQ=255) tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa > 1:28209/1‑149 (MQ=255) aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa < 1:78204/130‑1 (MQ=255) aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa > 2:78204/1‑130 (MQ=255) aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa > 2:265575/1‑149 (MQ=255) ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCTCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata < 3:102233/128‑1 (MQ=255) ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata > 4:102233/1‑128 (MQ=255) ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa < 3:261909/148‑1 (MQ=255) ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa > 3:10298/1‑148 (MQ=255) ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa > 3:66996/1‑148 (MQ=255) gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg > 3:27684/1‑148 (MQ=255) gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc < 4:191913/148‑1 (MQ=255) aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTa > 1:216984/1‑149 (MQ=255) gcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTaaa < 2:24693/149‑1 (MQ=255) cTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGc < 1:104589/149‑1 (MQ=255) gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGcc < 4:112112/148‑1 (MQ=255) gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGcc > 3:155788/1‑148 (MQ=255) aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAGTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAg > 1:118432/1‑149 (MQ=255) acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCa > 2:168398/1‑144 (MQ=255) acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCa < 1:168398/144‑1 (MQ=255) acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAg > 2:130581/1‑149 (MQ=255) acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGCCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCa > 4:136696/1‑148 (MQ=255) caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTc < 3:203341/146‑1 (MQ=255) caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTc > 4:203341/1‑146 (MQ=255) caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAg > 4:87027/1‑148 (MQ=255) caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAg > 4:117253/1‑148 (MQ=255) aaaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCg > 1:225542/1‑149 (MQ=255) tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTcc > 1:130744/1‑149 (MQ=255) ttttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGt < 4:405256/148‑1 (MQ=255) tttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAATCGAATGTCAAAAGGCACAGTAAAg > 3:333305/1‑124 (MQ=255) tttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAg < 4:333305/124‑1 (MQ=255) tttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTcc < 1:139756/149‑1 (MQ=255) ttACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTgg > 1:294764/1‑149 (MQ=255) tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTgg > 3:276034/1‑148 (MQ=255) tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTgg > 4:121720/1‑148 (MQ=255) aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt < 2:216984/149‑1 (MQ=255) aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt > 1:43445/1‑149 (MQ=255) acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt > 3:305425/1‑148 (MQ=255) caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTc < 2:113143/149‑1 (MQ=255) acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCaa < 4:175901/148‑1 (MQ=255) acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCaa < 4:27684/148‑1 (MQ=255) acaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAAt > 2:139273/1‑149 (MQ=255) caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTAGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTg > 4:396997/1‑138 (MQ=255) caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTAGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTg < 3:396997/138‑1 (MQ=255) | TCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAAT > NC_000913/4227412‑4227691 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |