Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,970,032 | (G)6→8 | 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 | 100.0% | 152.1 / NA | 42 | V332V (GTC→GTG) | wzzE | Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein |
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (19/23); total (19/23) | |||||||||||
* | NC_000913 | 3,970,026 | 2 | . | G | 95.2% | 142.7 / ‑0.2 | 42 | coding (996/1047 nt) | wzzE | Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base G (18/22); total (19/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 = 3.86e-01 |
GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACT > NC_000913/3969888‑3970147 || ggCCATGTT‑AACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑g > 1:535774/1‑139 (MQ=255) tGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggggg > 2:1097/1‑139 (MQ=255) ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggggggggctgatcg > 2:980476/1‑130 (MQ=255) ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggggggggctgatcg > 1:411588/1‑130 (MQ=255) ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggggggg > 2:444694/1‑130 (MQ=255) ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggggggg < 1:444694/131‑2 (MQ=255) ggATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGgtg > 1:1190760/1‑127 (MQ=255) ggATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGgtg < 2:1190760/127‑1 (MQ=255) gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGgtgt < 2:1096561/122‑1 (MQ=255) gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGgtgt > 1:1096561/1‑122 (MQ=255) gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGtt > 1:241039/1‑139 (MQ=255) cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTg < 2:535774/139‑1 (MQ=255) aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa > 2:562841/1‑131 (MQ=255) aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa < 1:562841/131‑1 (MQ=255) aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa < 2:241039/139‑1 (MQ=255) aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa < 1:1402114/139‑1 (MQ=255) cGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca < 1:980476/139‑1 (MQ=255) tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGaaa > 2:382985/1‑124 (MQ=255) tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTAGGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGaaa < 1:382985/124‑1 (MQ=255) gCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCa > 1:294922/1‑102 (MQ=255) gCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCa < 2:294922/102‑1 (MQ=255) gCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑gggggggg > 1:1071942/1‑79 (MQ=255) gCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑gggggggg < 2:1071942/81‑3 (MQ=255) ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGgcgc > 2:436956/1‑139 (MQ=255) ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGgcgc > 2:1172229/1‑139 (MQ=255) tAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt < 2:666670/98‑1 (MQ=255) tAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt > 1:666670/1‑98 (MQ=255) aCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc > 1:1044379/1‑95 (MQ=255) aCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc < 2:1044379/95‑1 (MQ=255) cgcgATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccgcct > 1:1521363/1‑137 (MQ=255) cgcgATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt > 2:1301617/1‑130 (MQ=255) cgcgATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt < 1:1301617/130‑1 (MQ=255) aTAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCgggg > 1:114959/1‑62 (MQ=255) aTAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCgggg < 2:114959/62‑1 (MQ=255) aTAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTc > 2:194422/1‑139 (MQ=255) tAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGaaa > 2:332224/1‑93 (MQ=255) tAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGaaa < 1:332224/93‑1 (MQ=255) gtcgtGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGTTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt < 1:1097/139‑1 (MQ=255) gtcgtGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCCTATTCGAAGAGAAt < 2:1521363/139‑1 (MQ=255) gtGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa < 2:411588/139‑1 (MQ=255) tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg > 2:574112/1‑119 (MQ=255) tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg < 1:574112/119‑1 (MQ=255) tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAt > 1:1563627/1‑139 (MQ=255) gCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa < 2:800636/86‑1 (MQ=255) gCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa > 1:800636/1‑86 (MQ=255) atTATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga < 1:194422/139‑1 (MQ=255) ttATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac < 1:1172229/139‑1 (MQ=255) tATGCGGGGCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgact < 1:436956/139‑1 (MQ=255) ggggCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCTGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAaga > 2:1243445/1‑111 (MQ=255) ggggCATTGTCGGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAaga < 1:1243445/111‑1 (MQ=255) || GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACT > NC_000913/3969888‑3970147 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |