| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,120,432 | 1 | . | C | 45.7% | 45.4 / 72.8 | 35 | coding (1123/1395 nt) | wcaJ | colanic biosynthesis UDP‑glucose lipid carrier transferase |
| Reads supporting (aligned to +/- strand): ref base . (9/10); new base C (7/9); total (16/19) | |||||||||||
| 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.88e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGATCGTTCTGCGTCGCCTGGGTCACCACTTTGTCGTTCTCCAT > NC_000913/2120297‑2120558 | cgccgcGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑c‑cc > 2:167953/1‑139 (MQ=255) gcCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGt < 2:405667/139‑1 (MQ=255) ccGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGca < 2:841241/139‑1 (MQ=255) ccGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGca < 2:1230452/139‑1 (MQ=255) cgcCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAAc > 2:944660/1‑139 (MQ=255) cAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCg < 1:15628/139‑1 (MQ=255) aGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGca > 1:126074/1‑123 (MQ=255) aGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGca < 2:126074/123‑1 (MQ=255) gCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATcc > 1:196970/1‑139 (MQ=255) tCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGt > 2:1144194/1‑139 (MQ=255) tGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGgcgc > 1:489759/1‑139 (MQ=255) tCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTa > 1:703518/1‑115 (MQ=255) tCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTa < 2:703518/115‑1 (MQ=255) cAACGAGCTGTCGATACTGTTCGTTGTGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGAACGGCGCAGAAAGTTCCCCACTTTGGt < 2:850991/139‑1 (MQ=255) aaTGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGt > 1:427297/1‑139 (MQ=255) cGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGTCCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGa < 2:37261/139‑1 (MQ=255) cGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGa < 1:167953/139‑1 (MQ=255) cTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGATCGtt < 2:983553/139‑1 (MQ=255) tCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTccc > 2:1135590/1‑111 (MQ=255) tCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTccc < 1:1135590/111‑1 (MQ=255) aTGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGATCGTTCTGCGTCGcc < 2:188954/139‑1 (MQ=255) tACTGCGTGCGGACGTGGACCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGATCGTTCTGCGTCGCCTGGGt < 1:275296/139‑1 (MQ=255) tgcgGACGTGGACCGACGATCGACAT‑CGCCCCCGGTCAGCCCATTGATAAACTGCGGCAATTAATCCAGCGCGGTACGGCGCAGCAAGTTCCCCACTTTGGTGCCCCGCGGATCGTTCTGCGTCGCCTGGGCCACCACt > 1:734118/1‑139 (MQ=255) tgcgGACGTGGACCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGATCGtt > 1:313024/1‑117 (MQ=255) tgcgGACGTGGACCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGATCGtt < 2:313024/117‑1 (MQ=255) tgcgGACGTGGACCGACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGATCGTTCTGCGTCGCCTGGGTCACCACt > 1:453624/1‑139 (MQ=255) tgcgGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGATCGTTCTGCGTCGCCTGGGTCACCACtt > 1:847351/1‑139 (MQ=255) gTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGATCGTTCTGCGTCGCCTGGGTCACCACTTTGTCGtt < 2:196970/139‑1 (MQ=255) tGGACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGAc < 1:148385/96‑1 (MQ=255) tGGACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGAc > 2:148385/1‑96 (MQ=255) gACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGATCGTTCTGCGTCGCCTGGGTCACCACTTTGTCGTtctc > 1:534543/1‑139 (MQ=255) gACAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGATCGTTCTGCGTCGCCTGGGTCACCACTTTGTCGTTCTCCAt < 2:453624/139‑1 (MQ=255) aCAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTa < 2:26473/61‑1 (MQ=255) aCAATCGACATCC‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTa > 1:26473/1‑61 (MQ=255) aCAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTAc > 2:610666/1‑61 (MQ=255) aCAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTAc < 1:610666/61‑1 (MQ=255) | CGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCGACAT‑C‑CCCCCGGTCAGCACATTGATAAACTGCGGCAATTCATCCAGCGAGGTACGGCGCAGAAAGTTCCCCACTTTGGTGACGCGCGGATCGTTCTGCGTCGCCTGGGTCACCACTTTGTCGTTCTCCAT > NC_000913/2120297‑2120558 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |