| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 417,378 | A→G | 19.9% | Q362R (CAG→CGG) | waaG → | lipopolysaccharide glucosyltransferase I |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 417,378 | 0 | A | G | 19.9% | 55.5 / 9.3 | 32 | Q362R (CAG→CGG) | waaG | lipopolysaccharide glucosyltransferase I |
| Reads supporting (aligned to +/- strand): ref base A (14/11); new base G (3/3); total (17/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01 | |||||||||||
GACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCA > NC_002947/417230‑417526 | gACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCTTGCTATGGTCTGGCATGTGCTGATAGCGCCGATCTGTACAGCCTGCCgcag > 2:312346/1‑150 (MQ=255) gCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTAAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCAtgctg < 1:290260/150‑1 (MQ=255) tGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgtac < 2:575037/150‑4 (MQ=255) gCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGt > 1:115115/1‑150 (MQ=255) ttCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgta > 2:750282/1‑130 (MQ=255) ttCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgta < 1:750282/132‑3 (MQ=255) ttCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCaa > 2:332571/1‑150 (MQ=255) aaCAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTggg < 1:564929/142‑1 (MQ=255) aaCAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTggg > 2:564929/1‑142 (MQ=255) aaGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCg < 2:621237/150‑1 (MQ=255) gcagcTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATg > 1:209930/1‑150 (MQ=255) gcgcATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAAccc > 1:605847/1‑150 (MQ=255) gcgcATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGACCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAAccc > 1:809173/1‑150 (MQ=255) aTGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCa > 2:628191/1‑150 (MQ=255) cTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGACCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAAc < 2:809173/150‑1 (MQ=255) tCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACg > 1:205826/1‑150 (MQ=255) gAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCCATCTGTACAGCATGCCGAAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGc‑c < 1:700069/150‑1 (MQ=255) aaGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑Ct > 2:744196/1‑150 (MQ=255) cGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgta > 2:638964/1‑81 (MQ=39) cGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgta < 1:638964/83‑3 (MQ=39) ccAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGc < 1:628191/150‑1 (MQ=255) cAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGcc > 2:191208/1‑150 (MQ=255) ggCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTg < 1:332571/150‑1 (MQ=255) cTGATAGCGCCGATCTGTAAAGTCTGCCGCTGCATGCAGCCGCTGTGATCCCAGCCCACCCCACCGCTCGGGGCTGATTGTTGCCGATGTTTTTAAATGC‑TGGTTGGGCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGa < 2:736500/150‑1 (MQ=255) tGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg > 1:736500/1‑150 (MQ=255) tGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg < 1:91453/150‑1 (MQ=255) cgcCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGCTCTCGACCCATCCCTT‑AAACGCTGTTTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTAt < 1:744196/150‑1 (MQ=255) aTCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGa > 1:494590/1‑150 (MQ=255) cggcAGCATGCTGCGGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTg > 1:485748/4‑150 (MQ=255) tACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGCCTTCGa > 2:4541/1‑109 (MQ=255) tACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGCCTTCGa < 1:4541/109‑1 (MQ=255) ccaagCCGAAGCATGCTGCCGACGTGATCCTTGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGAGACGGGCGGGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAgg < 2:209930/146‑1 (MQ=255) ccGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGcgcc > 2:532571/1‑150 (MQ=255) cagcaTGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCa < 2:605847/150‑1 (MQ=255) | GACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGC‑CTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCA > NC_002947/417230‑417526 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |