Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,125,627 | T→C | A40A (GCA→GCG) | wcaI ← | putative colanic biosynthesis glycosyl transferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,125,627 | 0 | T | C | 97.9% | 163.8 / ‑4.6 | 48 | A40A (GCA→GCG) | wcaI | putative colanic biosynthesis glycosyl transferase |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (21/26); total (22/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.58e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01 |
ACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTAC > NC_000913/2125494‑2125757 | aCAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTa < 2:902489/140‑1 (MQ=255) aGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAAt > 2:718601/1‑140 (MQ=255) gCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGa > 1:665696/1‑140 (MQ=255) cGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGgcggcgcg > 1:1468858/1‑132 (MQ=255) cGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGgcggcgcg < 2:1468858/132‑1 (MQ=255) cGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGAc > 1:910136/1‑140 (MQ=255) tCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGc < 2:1527013/140‑1 (MQ=255) tCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGc < 2:469853/140‑1 (MQ=255) ttttGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTtgtg < 1:1505571/140‑1 (MQ=255) ttGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCTCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTgcc < 2:1422959/140‑1 (MQ=255) tGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTgccg < 2:163297/140‑1 (MQ=255) gCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTgccgcc < 2:1050076/140‑1 (MQ=255) caTACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCgccagc < 1:727061/140‑1 (MQ=255) caTACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCgccagc < 2:147795/140‑1 (MQ=255) caTACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCgccagc < 2:227247/140‑1 (MQ=255) tGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTAt > 2:185603/1‑140 (MQ=255) tGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTAt > 1:843328/1‑140 (MQ=255) tGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTAt > 2:382703/1‑140 (MQ=255) gcggcCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTAttt > 1:573158/1‑140 (MQ=255) ccccTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGa > 2:275413/1‑140 (MQ=255) cccTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGAt < 1:370337/140‑1 (MQ=255) ccTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATg < 2:665696/140‑1 (MQ=255) ctcGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGc < 1:590400/129‑1 (MQ=255) ctcGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGc > 2:590400/1‑129 (MQ=255) gTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCg < 1:275413/140‑1 (MQ=255) tagagcCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCgg > 2:367158/5‑140 (MQ=255) tAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATg > 1:554862/1‑126 (MQ=255) tAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATg > 1:915215/1‑126 (MQ=255) tAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATg < 2:554862/126‑1 (MQ=255) tAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATg < 2:915215/126‑1 (MQ=255) tAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCg > 1:452655/1‑129 (MQ=255) tAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCg < 2:1269791/129‑1 (MQ=255) tAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCg > 1:1269791/1‑129 (MQ=255) tAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCg < 2:452655/129‑1 (MQ=255) tAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCgg > 2:1282724/1‑140 (MQ=255) cgcCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGa < 1:1282724/140‑1 (MQ=255) ccAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGagt > 1:316832/1‑140 (MQ=255) aGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGagtag > 1:439750/1‑140 (MQ=255) ggAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTaa < 2:129662/140‑1 (MQ=255) ggAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTaa < 1:367158/140‑1 (MQ=255) aaTAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATg > 1:175675/1‑140 (MQ=255) aGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCg < 2:439750/140‑1 (MQ=255) gCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCa > 2:710337/1‑140 (MQ=255) gCCCACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCa > 1:1351566/1‑140 (MQ=255) ccACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGt < 2:316832/140‑1 (MQ=255) cACCTGCCATTGCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTa < 2:843328/140‑1 (MQ=255) gCGGGTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATc > 1:185241/1‑140 (MQ=255) ggTAGTAAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCgg < 1:710337/140‑1 (MQ=255) tagtaAGGCGGCGCGGTAATGACCAGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGt > 2:591724/1‑140 (MQ=255) gtaAGGCGGCGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTAc < 1:167412/140‑1 (MQ=255) | ACAGGCGTTTCAGGGTGCTCGGCTGTTTTGGCACATACAGCGGGCAGCGCCACACCGTGGCGGCCCCCTCTTCTCGTTTGTAGCGCCAGGCGGAATAGTTCTCGCCCACCTGCCATTGCGGGTAGTAAGGCGGTGCGGTAATGACCCGCACCTCATGACCTTGTGCCGCCAGCCATTCCACCATCTCGCCGGTGTATTTGCCGATGCCGGTTAACTCCGGCGAGTAGTTAATGCCGTAGACCAGTATTTTCATAATCCGGGTAC > NC_000913/2125494‑2125757 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |