| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,188,038 | 0 | C | . | 68.6% | 43.6 / 50.9 | 35 | coding (377/1035 nt) | yehA | putative fimbrial‑like adhesin protein |
| Reads supporting (aligned to +/- strand): ref base C (5/6); new base . (13/11); total (18/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.25e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGT > NC_000913/2187900‑2188150 | tttCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGAGGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACcgccgc < 2:417113/139‑1 (MQ=255) gTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGcccccg < 1:421406/139‑2 (MQ=255) tAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCg > 2:373393/1‑139 (MQ=255) cacaTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACAtt > 1:404914/1‑138 (MQ=255) cacaTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACAtt < 2:404914/138‑1 (MQ=255) aTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCAc > 1:398186/1‑139 (MQ=255) tGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACtt > 2:194845/1‑139 (MQ=255) tCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGc > 2:370941/1‑139 (MQ=255) cccATTACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGTCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGc < 2:129164/139‑1 (MQ=255) gCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGcc > 1:19696/1‑139 (MQ=255) gCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCg > 2:378262/1‑139 (MQ=255) aaCACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCAc < 2:398186/139‑1 (MQ=255) caGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCgccgcc > 2:559626/1‑90 (MQ=255) caGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCgccgcc < 1:559626/90‑1 (MQ=255) aGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCTACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTaa < 1:373393/139‑1 (MQ=255) tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCAc > 1:587736/1‑131 (MQ=255) tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCAc < 2:587736/131‑1 (MQ=255) tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGAGTCACCGCAAATg > 1:270950/1‑139 (MQ=255) aTTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACAtt > 2:514679/1‑96 (MQ=255) aTTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACAtt < 1:514679/96‑1 (MQ=255) tGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTgg > 1:369934/1‑139 (MQ=255) tGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTgg < 1:584520/139‑1 (MQ=255) gcgTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAAc > 2:92802/1‑133 (MQ=255) gcgTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAAc < 1:92802/133‑1 (MQ=255) cgTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGcc > 2:263143/1‑96 (MQ=255) cgTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGcc < 1:263143/96‑1 (MQ=255) tgaAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACa < 2:270950/139‑1 (MQ=255) aaGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAAc > 2:471204/1‑139 (MQ=255) tGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGAtg < 2:369934/139‑1 (MQ=255) tAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACtttt > 2:145114/1‑60 (MQ=255) tAATGTCGCAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACtttt < 1:145114/60‑1 (MQ=255) cAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTAtt < 2:19696/139‑1 (MQ=255) cAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACt > 2:145354/1‑60 (MQ=255) cAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACt < 1:145354/60‑1 (MQ=255) cAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACa > 2:125436/1‑119 (MQ=255) cAACACCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACa < 1:125436/119‑1 (MQ=255) acacCATTAATCATCACTGAACCGCCG‑CCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGt < 1:194845/139‑1 (MQ=255) caCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTa < 2:385119/64‑1 (MQ=255) aaCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTa > 1:385119/2‑64 (MQ=255) | TTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGT > NC_000913/2187900‑2188150 |
| 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: ‑ |