Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,160,919 | C→A | 9.8% | C219* (TGC→TGA) ‡ | pabC → | 4‑amino‑4‑deoxychorismate lyase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,160,919 | 0 | C | A | 9.8% | 156.0 / 6.3 | 51 | C219* (TGC→TGA) ‡ | pabC | 4‑amino‑4‑deoxychorismate lyase |
Reads supporting (aligned to +/- strand): ref base C (21/25); new base A (3/2); total (24/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.56e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01 |
GGGTGATCGAAGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGTTGCCCGTACG > NC_002947/2160771‑2161063 | gggTGATCGAAGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGc < 4:157207/149‑1 (MQ=255) ggTGATCGAAGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGc > 4:182310/1‑148 (MQ=255) gTGATCGAAGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGa > 4:64891/1‑149 (MQ=255) tGATCGAAGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGAc < 5:342254/149‑1 (MQ=255) aaGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCGCAGACCTGtac < 3:313008/149‑3 (MQ=255) gTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTc > 2:134129/1‑147 (MQ=255) gTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGa < 8:340064/149‑1 (MQ=255) gTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGa < 2:365172/149‑1 (MQ=255) gTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCCGTACAGGTCGCAGACCTGtaccgcaa > 1:338552/1‑142 (MQ=255) tataCAGCAATCTGTTTCTGGTGCGCGACGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAg < 4:322644/149‑1 (MQ=255) taCAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCGCAGACCTGtac > 7:288719/1‑139 (MQ=255) taCAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCGCAGACCTGtac < 8:288719/141‑3 (MQ=255) aGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCt > 7:200717/1‑149 (MQ=255) aaTCTGTTTCTGGTGCGCGTTGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCGCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGa < 6:153660/149‑1 (MQ=255) aTCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAg > 2:227455/1‑149 (MQ=255) tttCTGGTGTGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGc < 1:347190/149‑1 (MQ=255) cTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGa < 1:224804/131‑1 (MQ=255) cTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGa > 2:224804/1‑131 (MQ=255) cTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGaa > 8:279238/1‑149 (MQ=255) cGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCTGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTAtttgttt < 5:136091/149‑1 (MQ=255) ctTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATa < 7:71117/148‑1 (MQ=255) aCCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGc < 1:134129/147‑1 (MQ=255) cgccgACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAg > 2:99265/1‑128 (MQ=255) cgccgACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAg < 1:99265/128‑1 (MQ=255) gccgACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTg > 1:401186/1‑111 (MQ=255) gccgACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTg < 2:401186/111‑1 (MQ=255) cTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCAt > 6:372115/1‑149 (MQ=255) gCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTgg > 2:134439/1‑149 (MQ=255) ccGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGc < 3:182310/149‑1 (MQ=255) tggcgtggcgGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCcgtgc < 3:64891/149‑1 (MQ=255) gcgtggcgGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCcgtgcgt > 2:93937/1‑149 (MQ=255) gtggcgGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTgg > 8:247435/1‑149 (MQ=255) ggcgTGATGCGTGCCGCACTGCTGGAACAGGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAg > 1:239249/1‑149 (MQ=255) ggcgTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGc > 5:157826/1‑149 (MQ=255) gcgTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGcc < 2:131149/149‑1 (MQ=255) aTGCGTGCCGCACTGCGGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCgcgc < 3:85949/148‑1 (MQ=255) cGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTaaa < 7:247435/149‑1 (MQ=255) ccGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTgg > 6:103681/1‑149 (MQ=255) ctgGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCg < 5:103681/147‑1 (MQ=255) ctgGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCg > 8:102098/1‑147 (MQ=255) ctgGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCg < 7:102098/147‑1 (MQ=255) tgGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTaaa > 7:50663/1‑136 (MQ=255) tgGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTaaa < 8:50663/136‑1 (MQ=255) cccGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTc > 2:117287/1‑149 (MQ=255) ctgGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTaaa < 1:227455/147‑1 (MQ=255) ctgGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACt < 6:157826/149‑1 (MQ=255) ggTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTg < 1:93937/147‑1 (MQ=255) tATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAgg < 8:200717/149‑1 (MQ=255) gCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGt > 5:55315/1‑149 (MQ=255) cGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGtt < 6:55315/149‑1 (MQ=255) cgaCAGGTCGCAGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGTTGc > 5:104825/3‑149 (MQ=255) tCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGTTGCCCGTACg > 8:44094/1‑149 (MQ=255) tCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGTTGCCCGTACg > 2:366384/1‑149 (MQ=255) | GGGTGATCGAAGGGGTATACAGCAATCTGTTTCTGGTGCGCGATGGCGTGCTGCTTACCGCCGACCTTAGCCGTTGTGGCGTGGCGGGCGTGATGCGTGCCGCACTGCTGGAACA‑GGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGTTGCCCGTACG > NC_002947/2160771‑2161063 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |