Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,743,058 | T→A | 7.2% | V109E (GTG→GAG) | yiaK → | 2,3‑diketo‑L‑gulonate reductase, NADH‑dependent |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,743,058 | 0 | T | A | 7.2% | 159.2 / 3.9 | 56 | V109E (GTG→GAG) | yiaK | 2,3‑diketo‑L‑gulonate reductase, NADH‑dependent |
Reads supporting (aligned to +/- strand): ref base T (24/26); new base A (2/2); total (27/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 = 4.68e-01 |
GATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCC > NC_000913/3742910‑3743205 | gATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCATTTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGTCAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCtggt < 5:202377/149‑1 (MQ=255) gATATCATTCCTGATGACAAACCCAAACGTATAACAAGCCTCTGCGCAATTGAACAGTGGGACGCCCAGCTTTCGAACGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCAATTGACATGGCTGCCTATCACTGTAATGGTCtggt > 4:114639/1‑149 (MQ=255) aTTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCAc > 1:191622/1‑148 (MQ=255) aTTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCAc < 5:181385/148‑1 (MQ=255) ccTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTAc < 1:150530/148‑1 (MQ=255) ccTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGAACGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTAc < 5:141543/148‑1 (MQ=255) ccaaaccaaaCGAATAACAAGCATATGCGCAATTCAACAGACCGACGCACAGCGATCGATCGGTAACCTGACAGCCAAAAACATGATGGATCGCGCCATAGACCTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGcc > 5:166305/1‑148 (MQ=255) caacccaaACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCaa > 4:225699/1‑149 (MQ=255) ccaaACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCAc > 3:199129/1‑149 (MQ=255) aaCGTATAACCAGCCGCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTgg < 8:211828/149‑1 (MQ=255) gTATAACCCGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATg > 6:89934/1‑149 (MQ=255) gTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATg > 6:149726/1‑149 (MQ=255) taACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATgcg < 6:159597/148‑1 (MQ=255) aaCCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGcgcg > 5:170555/1‑149 (MQ=255) ccAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCgcggc < 4:43312/149‑1 (MQ=255) aGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGcggcgg < 6:114932/149‑1 (MQ=255) gCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGcggcgg > 6:175686/1‑148 (MQ=255) ccTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCa < 2:99529/149‑1 (MQ=255) aaTTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCt > 4:112543/1‑149 (MQ=255) gACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGACTGGc < 7:154569/139‑1 (MQ=255) gACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGACTGGc > 8:154569/1‑139 (MQ=255) gTTCGATCGGTAATTTGACAGCGTAATAGATGATGGATCGCGCTATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTTCGTAATGCCAACCTCTGGATGCGCGGCGGCAGTTTCGGCTGGCAGGCGGCGGaaata < 4:190228/143‑3 (MQ=255) gTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGaaaaa > 3:190228/1‑143 (MQ=255) gTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCta > 6:87411/1‑148 (MQ=255) gATCGGTAACCTGACAGCGAAAAAGATGATGGATGGGGACATTGAACTGGCTGCCGATCACGGTATTTGTCTGGTGGCACTACGTAATGCCATCCACTGGATGCGCGGCGGCTGCTACGGGTGGCAGGCGGCGGAAAAATGCTATAtt < 6:166305/148‑1 (MQ=255) gATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTg > 1:102980/1‑149 (MQ=255) gATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTg < 6:210353/149‑1 (MQ=255) agcGGTAACCTGACAGCGAAAACGATGATGGAGCGCGCCATTGCACGGGCTGCCGATCACGGTATTGGTCTGGGGGCACTACGTAAAGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTgg < 3:161393/147‑1 (MQ=255) tCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGAGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTgg > 7:34493/1‑148 (MQ=255) tCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGAGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTgg > 7:73320/1‑148 (MQ=255) gTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAggcgg > 8:219775/1‑126 (MQ=255) aCCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGc > 6:168125/1‑149 (MQ=255) cTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTg < 5:168125/149‑1 (MQ=255) cTGACAGCGAAAAAGATGATGGATAGAGCAATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTTGCATTTGCTg < 4:199129/149‑1 (MQ=255) tGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTATGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTg > 7:143948/1‑148 (MQ=255) aCAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGa < 5:175686/148‑1 (MQ=255) aCAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGa < 5:199703/148‑1 (MQ=255) aCAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGAGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTgg < 8:73320/147‑1 (MQ=255) aCAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGAGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTgg < 8:34493/147‑1 (MQ=255) gAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACt > 1:184922/1‑149 (MQ=255) aaaaaGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTc < 6:55830/149‑1 (MQ=255) aaaaGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTcc < 6:170555/149‑1 (MQ=255) atgGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCg > 1:108209/1‑148 (MQ=255) ggATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACt > 6:102628/1‑68 (MQ=255) ggATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACt < 5:102628/68‑1 (MQ=255) aTCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGCGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCAGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATg > 2:202979/1‑149 (MQ=255) gAACTGGCTGCCGATCACGGTATTTGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTggg < 3:114639/148‑1 (MQ=255) aaCTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGGATTGCCGCCGTGGGgc > 6:168537/1‑149 (MQ=255) ggCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGc < 1:163555/103‑1 (MQ=255) ggCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGc > 2:163555/1‑103 (MQ=255) tGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCaa > 3:64354/1‑146 (MQ=255) tCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTc < 5:149726/149‑1 (MQ=255) gTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGc < 1:169406/95‑1 (MQ=255) gTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGc > 2:169406/1‑95 (MQ=255) gTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATa < 5:89934/149‑1 (MQ=255) tATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAACGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATgc < 8:149017/119‑1 (MQ=255) tATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAACGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATgc > 7:149017/1‑119 (MQ=255) ttGGGCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTTATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGc < 3:62478/149‑1 (MQ=255) gtggCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAAccc > 5:4709/1‑149 (MQ=255) | GATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCC > NC_000913/3742910‑3743205 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |