| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,576,374 | 0 | T | C | 21.4% | 52.6 / 10.4 | 28 | T416A (ACT→GCT) | yddB | putative TonB‑dependent outer membrane receptor |
| Reads supporting (aligned to +/- strand): ref base T (12/10); new base C (3/3); total (15/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATGCGATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCT > NC_000913/1576244‑1576501 | aTGCGATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGAt < 1:53113/149‑1 (MQ=255) gATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCAc < 1:63129/131‑1 (MQ=255) gATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCAc > 2:63129/1‑131 (MQ=255) taCAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGAGCCg > 6:42605/1‑147 (MQ=255) aCAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCTCGCATCGGAATAGATGTATTCCGCGCCGAAg > 5:5402/1‑149 (MQ=255) ttGTCAATTCCCAGGCGGCCTCTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGCCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATg > 1:22140/1‑149 (MQ=255) aaTTCCCAGGCGGCCTCTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGCCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGGGaaa < 2:22140/149‑1 (MQ=255) gcCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGCCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACgg > 2:7499/1‑149 (MQ=255) tttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGaaa < 4:69909/128‑1 (MQ=255) tttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGaaa > 3:69909/1‑128 (MQ=255) tGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGCCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACgg < 1:105209/136‑1 (MQ=255) tGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGCCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACgg > 2:105209/1‑136 (MQ=255) tGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGCCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCa < 1:7499/149‑1 (MQ=255) gtatggtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCTCGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACTTTACCCACGGCGAATTTTTGCCAGTCCAGGc < 6:5402/148‑1 (MQ=255) gtatggtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGc < 5:42605/148‑1 (MQ=255) tatggtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCgt > 2:151808/1‑149 (MQ=255) gtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGttt > 2:153281/1‑149 (MQ=255) gTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCAc < 5:169882/116‑1 (MQ=255) gTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCAc > 6:169882/1‑116 (MQ=255) ttttCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAg > 6:177462/1‑149 (MQ=255) tcCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTaa < 1:151808/148‑1 (MQ=255) acaTAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTggg > 2:84344/1‑148 (MQ=255) aTTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATgtgt < 1:153281/149‑1 (MQ=255) aGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCt > 8:33748/1‑148 (MQ=255) aGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCACCGAAGTAGGGTTGATGCGAAAc < 5:76601/72‑1 (MQ=255) aGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCACCGAAGTAGGGTTGATGCGAAAc > 6:76601/1‑72 (MQ=255) gCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGaa < 7:133170/75‑1 (MQ=255) gCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGaa > 8:133170/1‑75 (MQ=255) | ATGCGATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCT > NC_000913/1576244‑1576501 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |