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 | 97.9% | 162.6 / ‑3.3 | 47 | coding (996/1047 nt) | wzzE | Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (22/24); total (23/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.89e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01 |
GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑G‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTA > NC_000913/3969888‑3970157 | ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc > 2:1416976/1‑139 (MQ=255) ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCgg > 1:100274/1‑139 (MQ=255) tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCgggg > 1:249950/1‑139 (MQ=255) gAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGc < 1:873616/139‑1 (MQ=255) ggATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAAcc > 2:241798/1‑137 (MQ=255) ggATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAAcc < 1:241798/137‑1 (MQ=255) gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGgtg > 1:1401423/1‑120 (MQ=255) gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGgtg < 2:1401423/120‑1 (MQ=255) aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATaacaac > 2:1133442/1‑139 (MQ=255) gaagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCaaa > 1:1126110/1‑139 (MQ=255) aaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt < 1:1133442/139‑1 (MQ=255) gATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTc < 2:100274/139‑1 (MQ=255) aCGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATaac < 2:220294/90‑3 (MQ=255) aCGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATaac > 1:220294/1‑88 (MQ=255) aCGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGaa > 1:363347/1‑139 (MQ=255) cgtcgtGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAAGAACAACACTGCTGCg < 2:956820/101‑1 (MQ=255) cgtcgtGCCTTCCTGATGATTATGCGGGGCATTGTC‑GCGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTCCg > 1:956820/1‑101 (MQ=255) tGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg < 2:519292/139‑1 (MQ=255) gCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCAt < 2:614991/59‑1 (MQ=255) gCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCAt > 1:614991/1‑59 (MQ=255) gCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc > 1:694639/1‑75 (MQ=255) gCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc < 2:694639/75‑1 (MQ=255) gCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgt > 1:1010364/1‑139 (MQ=255) tgatgaTTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt < 1:450849/112‑1 (MQ=255) tgatgaTTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt > 2:450849/1‑112 (MQ=255) atgatTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCAt < 1:261586/50‑1 (MQ=255) atgatTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCAt > 2:261586/1‑50 (MQ=255) atgatTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTact < 1:78545/139‑1 (MQ=255) ttATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgact < 2:1126110/139‑1 (MQ=255) tATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa < 1:863538/63‑1 (MQ=255) tATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa > 2:863538/1‑63 (MQ=255) tATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt > 2:249723/1‑105 (MQ=255) tATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt < 1:249723/105‑1 (MQ=255) tATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac > 2:1032646/1‑137 (MQ=255) tATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac < 1:1032646/137‑1 (MQ=255) tgcggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAtt > 2:369240/4‑107 (MQ=255) tgcggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAtt < 1:369240/104‑1 (MQ=255) tgcggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga < 2:654244/131‑1 (MQ=255) tgcggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga > 1:654244/4‑134 (MQ=255) ggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt > 1:1398107/1‑55 (MQ=255) ggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt < 2:1398107/55‑1 (MQ=255) gCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGaa > 1:1358606/1‑120 (MQ=255) gCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGaa < 2:1358606/120‑1 (MQ=255) cATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt < 1:414046/89‑1 (MQ=255) cATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt > 2:414046/1‑89 (MQ=255) cATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt < 2:1429048/110‑1 (MQ=255) cATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt > 1:1429048/1‑110 (MQ=255) cATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTa < 1:1416976/139‑1 (MQ=255) | GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑G‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTA > NC_000913/3969888‑3970157 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |