Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,970,032 | (G)6→7 | coding (1002/1047 nt) | wzzE → | Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,970,026 | 1 | . | G | 92.7% | 132.0 / 3.0 | 41 | coding (996/1047 nt) | wzzE | Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein |
Reads supporting (aligned to +/- strand): ref base . (1/2); new base G (17/21); total (18/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01 |
GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTA > NC_000913/3969888‑3970157 | ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc > 2:514873/1‑139 (MQ=255) ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc > 2:1517433/1‑139 (MQ=255) cATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggg < 1:1498536/139‑1 (MQ=255) tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTACGTACGCCGGAAGAACCGGTAAAACGCGATAGCC‑‑CGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGAt < 2:555836/132‑1 (MQ=255) tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTACGTACGCCGGAAGAACCGGTAAAACGCGATAGCC‑‑CGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGAt > 1:555836/1‑132 (MQ=255) gAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACGTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGc < 2:1020774/139‑1 (MQ=255) ggTCCAACCCTGGATCCGCGTTTTCAGACCTAACGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTc > 2:1344310/1‑139 (MQ=255) aCCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTa > 2:1383835/1‑139 (MQ=255) ggATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACccg > 1:274422/1‑139 (MQ=255) gATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACccgc > 1:32668/1‑139 (MQ=255) tCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGGGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACccgccg < 1:660873/139‑1 (MQ=255) ttttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc < 2:274422/139‑1 (MQ=255) ttttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc < 2:163054/139‑1 (MQ=255) ttttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc < 1:1517433/139‑1 (MQ=255) cAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGaaa < 2:1162097/139‑1 (MQ=255) aCCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATaa < 1:1383835/139‑2 (MQ=255) tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAacac > 2:1293463/1‑139 (MQ=255) cGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAacac > 1:752603/1‑139 (MQ=255) tACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGAt > 1:1120/1‑82 (MQ=255) tACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGAt < 2:1120/82‑1 (MQ=255) cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGgt > 2:897096/1‑90 (MQ=255) cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gg > 1:588353/1‑70 (MQ=255) cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gg < 2:588353/70‑1 (MQ=255) cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑GGGTGGGCTGATCGGGGCTGgt < 1:897096/90‑1 (MQ=255) agaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGgc < 1:1293463/139‑1 (MQ=255) ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGc > 2:229391/1‑75 (MQ=255) ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGc < 1:229391/75‑1 (MQ=255) aaaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg < 2:752603/139‑1 (MQ=255) aaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCTTTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt < 1:514873/139‑1 (MQ=255) aTAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCg < 1:1351704/139‑1 (MQ=255) cgtcgtGCCTTCCTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCTTTAACCCGCCGTTGCTCGAAGTAACAGCACTGCTGCGGTGAGCGCAAAGGCGATCGCCGCTTATTCGAAGAGAAt < 1:34405/139‑1 (MQ=255) ccTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa > 2:64431/1‑131 (MQ=255) ccTGATGATTATGCGGGGCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa < 1:64431/131‑1 (MQ=255) tgcggggCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACccgccg > 1:499470/4‑53 (MQ=255) tgcggggCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACccgccg < 2:499470/50‑1 (MQ=255) tgcggggCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcg < 2:1003811/95‑1 (MQ=255) tgcggggCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcg > 1:1003811/4‑98 (MQ=255) tgcggggCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTac < 2:163637/128‑1 (MQ=255) tgcggggCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTac > 1:163637/4‑131 (MQ=255) ggCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga < 1:974464/129‑1 (MQ=255) ggCATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga > 2:974464/1‑129 (MQ=255) cATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt < 2:52847/99‑1 (MQ=255) cATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt > 1:52847/1‑99 (MQ=255) cATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa > 1:676028/1‑123 (MQ=255) cATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa < 2:676028/123‑1 (MQ=255) cATTGTCGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTa > 2:166871/1‑139 (MQ=255) | GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTA > NC_000913/3969888‑3970157 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |