Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,543,600 | A→G | 6.4% | A162A (GCA→GCG) | fimA → | major type 1 subunit fimbrin (pilin) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,543,600 | 0 | A | G | 6.4% | 186.7 / 2.9 | 62 | A162A (GCA→GCG) | fimA | major type 1 subunit fimbrin (pilin) |
Reads supporting (aligned to +/- strand): ref base A (32/26); new base G (2/2); total (34/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01 |
GCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTA > NC_000913/4543457‑4543745 | gCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCgg < 8:25268/149‑1 (MQ=255) cGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGTCTGAACGGGTGGTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCTGGCGCGTTATTTTGCAACCgtgg < 8:95117/147‑1 (MQ=255) ggTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAAc > 6:30637/1‑140 (MQ=255) ggTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAAc < 5:30637/140‑1 (MQ=255) aacaaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCGACCGGGGCCGCAAcccat < 7:49242‑M2/148‑6 (MQ=255) acaaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa > 8:42‑M2/1‑142 (MQ=255) caaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag < 8:353000‑M2/149‑9 (MQ=255) caaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag > 3:101357‑M2/1‑141 (MQ=255) caaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag < 8:83955‑M2/149‑9 (MQ=255) gTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACCGCCCTGAATAACGGAACCACTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCACcccataagcgct < 8:151335‑M2/148‑13 (MQ=255) tgCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCGACCGGGGCCGCAAcccataagcgctaacttaag < 5:12693‑M2/149‑21 (MQ=255) aGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt < 2:130642‑M2/149‑24 (MQ=255) cTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaa > 8:259667‑M2/1‑121 (MQ=255) tGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCCTTATTTTGCAACCGGGGCCCCAAcccataatcgctaacttaagggatgaac < 2:197016‑M2/148‑29 (MQ=255) gACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaacca > 8:32897‑M2/1‑118 (MQ=255) aCAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatc > 1:303033‑M2/1‑117 (MQ=255) aCGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAATCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagggctaacttatgggttgaaccatctgaa > 1:197016‑M2/1‑112 (MQ=255) ggTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCATAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaat > 4:9822‑M2/1‑109 (MQ=255) gTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatg > 2:13707‑M2/1‑108 (MQ=255) gTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccacctgaagaatg > 2:233330‑M2/1‑108 (MQ=255) ctgcGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcga < 8:88276‑M2/149‑45 (MQ=255) tgcGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgac < 7:212683‑M2/149‑46 (MQ=255) cgcTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgc > 4:107486‑M2/1‑102 (MQ=255) cTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcct > 1:378656‑M2/1‑100 (MQ=255) gCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCGTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtg < 7:318445‑M2/149‑55 (MQ=255) gATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcct > 3:86546‑M2/1‑91 (MQ=255) tGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgt > 7:31864‑M2/1‑89 (MQ=255) tGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgtt < 7:46037‑M2/147‑62 (MQ=255) gACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaacaatctgaagaatgcgacgcctcggtgcctcgttaagac < 7:153604‑M2/149‑67 (MQ=255) aTTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgat < 1:209270‑M2/149‑70 (MQ=255) tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatg < 6:361133‑M2/149‑71 (MQ=255) tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatg < 4:6416‑M2/149‑71 (MQ=255) ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcg < 8:23289‑M2/149‑76 (MQ=255) aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGTCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttc > 1:284359‑M2/1‑69 (MQ=255) aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttc < 6:199447‑M2/149‑81 (MQ=255) aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatc < 4:174112‑M2/101‑33 (MQ=255) aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttc < 2:145621‑M2/148‑84 (MQ=255) aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcg > 7:83955‑M2/1‑58 (MQ=255) aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgtt < 5:339140‑M2/149‑93 (MQ=255) ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCGACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgta > 6:12693‑M2/1‑52 (MQ=255) ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgta > 8:46037‑M2/1‑52 (MQ=255) ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttg < 7:205748‑M2/147‑96 (MQ=255) ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttg < 2:280947‑M2/147‑96 (MQ=255) gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCGACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtag > 8:49242‑M2/1‑51 (MQ=255) aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgt > 7:23289‑M2/1‑46 (MQ=255) ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcag > 6:339140‑M2/1‑42 (MQ=255) cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaat > 4:200903‑M2/1‑41 (MQ=255) cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaat < 3:200903‑M2/127‑87 (MQ=255) ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgtttt < 6:56547‑M2/133‑95 (MQ=255) ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgtttt > 5:56547‑M2/1‑39 (MQ=255) ccTTTCCAGTCGCGTTATTTTGCAACCGGGGCCGCAAcccataagggctatgttaagggtttaaccatctgaagaatgagacggctgggtgggtcgttaagacgatgcatcgagttcttcaattgcgttttgtaggctgtcagggata > 8:318445‑M2/1‑37 (MQ=255) ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggata > 2:209270‑M2/1‑37 (MQ=255) ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtc > 7:353000‑M2/1‑32 (MQ=255) ggCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccac > 7:88276‑M2/1‑29 (MQ=255) gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtttaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcatttgcgttttgtaggctgtcagggatactgtcccac > 4:147192‑M2/1‑28 (MQ=255) cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacg < 5:57090‑M2/148‑122 (MQ=255) cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgcagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaat > 6:180784‑M2/1‑25 (MQ=255) cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttctccaattgcgttttgtaggctgtcagggatactgtcccacgaat > 3:6416‑M2/1‑25 (MQ=255) gTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatg > 1:43950‑M2/1‑24 (MQ=255) ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatgg > 1:145621‑M2/1‑23 (MQ=255) aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggc > 8:17943‑M2/1‑21 (MQ=255) gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacct > 1:130642‑M2/1‑16 (MQ=255) | GCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTA > NC_000913/4543457‑4543745 |
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: ‑ |