| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 324,033 | 0 | T | G | 56.5% | 17.5 / 22.2 | 23 | T138P (ACC→CCC) | betB | betaine‑aldehyde dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base T (3/7); new base G (13/0); total (16/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.89e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGC > NZ_CP009273/323945‑324117 | aGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGta < 2:93013/90‑1 (MQ=255) tttCCACAGGCCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGt < 2:158879/90‑1 (MQ=255) acaGGGCAATCGGGATCGGGTAGTTCCACGCCCAAACCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCAc < 1:589418/90‑1 (MQ=255) ggATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCt < 2:709781/90‑1 (MQ=255) gATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGGATACACAAAGGACGTTTTACGCAACGGGAGCTg > 2:361525/1‑90 (MQ=255) tCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGtat > 2:326996/1‑57 (MQ=255) tCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGtat < 1:326996/57‑1 (MQ=255) tCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGc < 2:543339/90‑1 (MQ=255) gggTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGGTCTGGCTg > 2:717348/1‑90 (MQ=255) gggTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGGATACACAAAAGACGTTTCACGCCACGGGATCTGGCTg > 1:500015/1‑90 (MQ=255) ggTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGGATACACAAAAGACGTTTCCCGCAACCGGAGCTGGCTGc > 1:210255/1‑90 (MQ=255) gTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGGATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTc > 2:730171/1‑90 (MQ=255) aTGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGGATAAACAAAAGACGTTTCCCCCAACCGGGTCTGGCTGCCTTCCcgcg > 1:523094/1‑90 (MQ=255) aTGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGGATAAACAAAAGACGGTTCACCCCACGGGATCTTGCTGCCTTCCcgcg > 1:60234/1‑90 (MQ=255) accccTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATc < 2:275571/88‑1 (MQ=255) tGCCACTACGCCCAGCGGTTCGCGGCGGGGATACACAAAGGACGTTTTACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGGTCAGccc > 2:211175/1‑90 (MQ=255) tGCCACTACGCCCAGCGGTTCGCGGCGGGGATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCACCGCCGGGATCAGccc > 2:82999/1‑90 (MQ=255) tGCCACTACGCCCAGCGGTTCGCGGCGGGGATACACAAAGGACGGTTTACGCAACGGGAGCTGGCTGCCTTCCCGCGCCCGGATCAGccc > 2:170708/1‑90 (MQ=255) aCTACGCCCAGCGGTTCGCGGCGGGGATACACAAAAGACGTTTCACGCAACGGGATCTGGCTGCCTTCCCGCGCCGGGGTCAGCCCcgcg > 2:8400/1‑90 (MQ=255) tACGCCCAGCGGTTCGCGGCGGGTATACACAAAAGACGTTTCACGCAACGGGATCTGGCTGCCTTCCCGCGCCGGGGTCAGCCCGGCgga > 2:74041/1‑88 (MQ=255) tACGCCCAGCGGTTCGCGGCGGGGATACACAAAAGACGTTTCACGCAAAGGGATCTCGCTGCCTTCCCGCCCCGGGGGCCGCCCGGCgga > 1:163914/1‑88 (MQ=255) aCGCCCAGCGGTTCGCGGCGGGGATACACAAAAGACGTTTCACGCAACGGGGTCTTGCTGCCTTCCAGCGCCGGGGTCAGCCCCGCggag > 2:580725/1‑90 (MQ=255) cGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCgtagt < 1:397607/90‑1 (MQ=255) gCCCAGCGGTTCGCGGCGGGTATACACAAAAGGCGTTTTCCGCCACGGGGTCTGGCTGCCTTCCCGCGCCCGGGTCAGCCCGGCggagga > 1:343925/1‑88 (MQ=255) ccAGCGGTTCGCGGCGGGGATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGGGGAGTAcc > 2:712486/1‑89 (MQ=255) gcggGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGc < 2:60234/90‑1 (MQ=255) | AGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGC > NZ_CP009273/323945‑324117 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |