Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 4,531,167 | G→A | A132T (GCC→ACC) | fimB → | type 1 fimbria switch DNA invertase FimB |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,531,167 | 0 | G | A | 97.4% | 135.7 / ‑3.6 | 38 | A132T (GCC→ACC) | fimB | type 1 fimbria switch DNA invertase FimB |
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (16/21); total (17/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.47e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01 |
CAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGT > NZ_CP009273/4531032‑4531299 | cAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCg > 2:58285/1‑139 (MQ=255) aGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCgg > 2:38262/1‑139 (MQ=255) aaaaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCa < 2:447152/138‑1 (MQ=255) tACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATg < 1:38262/139‑1 (MQ=255) tACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGc > 2:239901/1‑138 (MQ=255) ccGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCAtt > 2:71304/1‑139 (MQ=255) cATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgt > 2:314450/1‑139 (MQ=255) aTGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgtg < 1:42462/139‑1 (MQ=255) cGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCt < 1:239901/139‑1 (MQ=255) aGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCttt > 2:581094/1‑139 (MQ=255) aGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCttt > 1:315047/1‑139 (MQ=255) aGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCttt > 2:291709/1‑139 (MQ=255) ggTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCg < 1:138377/139‑1 (MQ=255) tttttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAata < 1:71304/139‑1 (MQ=255) tttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAatat < 1:69547/138‑1 (MQ=255) aCGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAatat < 1:129875/132‑1 (MQ=255) aCGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAatat > 2:129875/1‑132 (MQ=255) aCGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAAt < 1:581061/138‑1 (MQ=255) tCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACt > 2:456288/1‑138 (MQ=255) tCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGAt < 2:415473/72‑1 (MQ=255) tCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGAt > 1:415473/1‑72 (MQ=255) tCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACt > 1:360724/1‑138 (MQ=255) ggCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGa < 1:60965/136‑1 (MQ=255) ggCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGa < 1:581094/136‑1 (MQ=255) ggCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGa < 1:291709/136‑1 (MQ=255) gTTTTACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTggg < 2:360724/139‑1 (MQ=255) ttttACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTggg > 1:175700/1‑138 (MQ=255) tACCATATTATCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCAt < 1:58285/138‑1 (MQ=255) tattatCTCGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCaa > 1:421151/1‑138 (MQ=255) cGACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCa > 1:457883/1‑139 (MQ=255) gACTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCa < 2:175700/138‑1 (MQ=255) aCTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAat < 2:294782/81‑1 (MQ=255) aCTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAat > 1:294782/1‑81 (MQ=255) aCTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATa < 1:314450/139‑1 (MQ=255) aCTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATa < 2:315047/139‑1 (MQ=255) aCTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATa < 1:547204/139‑1 (MQ=255) aCTTCCGGTGGTAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATa > 1:267671/1‑139 (MQ=255) gtAATACCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGt < 2:421151/138‑1 (MQ=255) | CAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGT > NZ_CP009273/4531032‑4531299 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |