Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB122 | 1,791,678 | 0 | T | G | 32.6% | 32.0 / 11.8 | 47 | L182R (CTG→CGG) | lpxB | Lipid‑A‑disaccharide synthase |
Reads supporting (aligned to +/- strand): ref base T (23/8); new base G (5/10); total (29/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.13e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACCTGATGCTGACGCTGCTGCCGTTCGAGGCGCGGTTCTACGAAGAGCAGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCCAACAAGTGCCGGGTGTGCGCTTCGTGTTG > ACB122/1791545‑1791811 | gtccTGATGCTGACGCTGCTGCCGTTCGAGGCGCGGTTCTACGAAGAGCAGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGctgg < 7:245089/134‑1 (MQ=255) ctgctgCCGTTCGAGGCGCGGTTCTACGAAGAGCAGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGc > 4:83177/1‑136 (MQ=255) gAGGCGCGGTTCTACGAAGAGCAGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGCCACCATCCCGCTTGAAGCCGACCGCCCGG‑CCGCGCGGGCCGCGCTGGGCCTGGGCGCGGGGCCGGtggtgggg > 1:191759/1‑133 (MQ=255) cgcgGTTCTACGAAGAGCAGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGAt > 4:120113/1‑136 (MQ=255) cgcgGTTCTACGAAGAGCAGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGAt > 4:114583/1‑136 (MQ=255) ggTTCTACGAAGAGCAGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGGTCCGCGCGGGGCCTGGGCGAGGGGCCGGTGGTGGCGGTggttc > 2:243433/1‑131 (MQ=255) ttCTACGAAGAGCAGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCgg > 8:302476/1‑136 (MQ=255) ttCTACGAAGAGCAGGGTGTGCCAGTCCGTTTTGTCGGCCACCCGCTGCCAGACACCATACCGCTTGAAGCCTACCGCCCGGCCCCCGCCTCCCCCCCGCGCCCCGGGCGCGGGCCCGTTGGGGCCGTTGATGccc > 8:296220/1‑135 (MQ=255) cTACGAAGAGCAGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCGGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGc < 3:146355/136‑1 (MQ=255) agagCAGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCgg > 7:172559/1‑136 (MQ=255) aGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CTGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTg > 2:64101/1‑135 (MQ=255) aGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGCTGCCGGGCAGGCGGGGTgg > 8:154036/1‑136 (MQ=255) aGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTgg > 7:151811/1‑136 (MQ=255) aGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCGGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTg > 2:237454/1‑135 (MQ=255) aGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCCCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCCTTGATGCCGGGCCGCCGGGGTgg > 3:236957/1‑136 (MQ=255) gtgtGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGGTGATGCCGGGCAGCCGGGGGGGCg > 7:13735/1‑135 (MQ=255) aGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGGCCCGCTGGGCCGGGGCGCGGGGCCGGTGGGGGCGGTGCTGCCGGGCAGCCGGGGTGGGGAAGTAggg > 6:41170/1‑136 (MQ=255) aGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAgg > 1:3985/1‑135 (MQ=255) aGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAgg > 1:272460/1‑135 (MQ=255) aGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCAGGGGCTGGGCGAGGGGGCGGTGGTGGGGTTGCTGCCGGGCAGCCGGGGTGGGGAAGTAgg > 2:301574/1‑135 (MQ=255) aGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGGGCCGCGGGGGGGCGGGGCGGGGGGGCGGGGGGGGCGGTGTGGGCGGGGCGCCCGGGGGGGGGAGTTggg > 5:184214/1‑136 (MQ=255) aGTGCGTTATGTCGGCCACCCGCTGGCAGACACCATACCGCAAGAAGCCGACCGCCCCG‑CCGCGCCTGCCGCGCTGGGCCTGGGCGAGCGGCCGGAGGAGGCGTTGATGCACGGCAGGCGGGGTGGCGAAGTAggg > 6:152829/1‑136 (MQ=255) gTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGGCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGCTGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCt > 2:224704/1‑135 (MQ=255) gTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGTTGCCGGGCCGCCGGGGGGGCGGAGTAGGGCGCtt > 6:278934/1‑136 (MQ=255) ttGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGGGCCGCGCTGGGCCGGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTgg < 6:94154/135‑1 (MQ=255) tGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGgc > 5:171442/1‑136 (MQ=255) gtggcccACCCGCGGCCACAAACCAACCCGCTTCAACCCGCCCGCCGGC‑CCGGGCGGCCCGCGCGGGGCCGGGGCGCGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGgc < 7:167264/130‑1 (MQ=255) cccGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCGGGGCGCGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGtt < 2:3985/135‑1 (MQ=255) cccGCTGGAACACACCATACCGCTTGACGCCGCCCGCCCGG‑CCGCGCGGGCCGCGCGGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGtt < 2:272460/135‑1 (MQ=255) gCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGC‑CCGCGGGTGCCGCGCGGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCt < 2:9646/135‑1 (MQ=255) aCCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCAGGGCGCGGGGCCGGTGGTGGCGTTGATGCCGGGCCGCCGGGGGGGCGAAGTAGGGCGCTTGGGCGCGCTGGTCCTCGCCGCCGCCGAgcgc > 7:129337/1‑136 (MQ=255) gCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCt > 4:275271/1‑135 (MQ=255) gAAGCCGACCGCCCGG‑CCGCGCGTGGCGCGCGGGGCCAGGGCGCGGGGCCGGTGGTGGCGTTGATGCCGGGCCGGCGGGGTGGCGGAGTAGGGCGCTTGGGCGCGCGGGTACGCGCCGCCGCCGCGCGccgccgc > 1:310990/1‑130 (MQ=255) gAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGGGTTGATGCCGGGCCGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGCGCGCCTCAGc > 2:110768/1‑135 (MQ=255) acgccgccCGCCGGG‑CGGGGCGGGCCGCGCGGGGCCGGGGGGCGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCca < 5:85161/134‑1 (MQ=255) aaGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGcc > 8:234453/1‑135 (MQ=255) gcggcccGCCCGG‑CGGCGGGGGCGGCGCTGGGCGGGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCcaac < 6:270697/131‑1 (MQ=255) gcccGCCCGG‑CCGCGCGGGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCCAACAAg < 5:62018/134‑1 (MQ=255) ccccgccgGG‑CGGCGGGGGCCGGGGGGGGCGGGGGGGCGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCCAACAAg < 5:238241/134‑1 (MQ=255) ggg‑CCGGGCGGGCCGCGCGGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCCAACAAGTGCCggg < 7:141221/135‑1 (MQ=255) tgcgcgTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCCAACAAGTGCCGGgtgt < 1:64101/134‑1 (MQ=255) cgggggTGCCGGGCGGGGCCGGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCCAACAAGTGCCGGgtgt < 1:243433/130‑1 (MQ=255) gcgcgGGCGGCGCTGGGCCGGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCCAACAAGTGCCGGGTGTgc < 7:46233/136‑1 (MQ=255) gcgcgGGCCGCGCTGGGCCGGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCCAACAAGTGCCGGGTGTgc < 7:46235/136‑1 (MQ=255) gcgcgGGCCGCGCGGGGCCGGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCCAACAAGTGCCGGGTGTgc < 3:83177/136‑1 (MQ=255) cgggggCCGCGGGGGGCGTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCCAACAAGTGCCGGGTGTgc < 1:114557/130‑1 (MQ=255) cgcgTGCCGCGCTGGGCCTGGGGGCGGGGGCGGTGGTGGGGTTGTTTCCGGGCCGCCGGGGGGGCGCAGTAGGGCGCTTGGGCGCGCGGTTGCCCGCCGCCGCCGCGCGCCGCCGCCACCCAGTGGCggggtgggg > 6:3252/1‑130 (MQ=255) gcTGGGCCGGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCCAACAAGTGCCGGGTGTGCGCTTCGTGTTg < 6:171442/136‑1 (MQ=255) tgggcctgggcGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGaa > 6:59302/1‑56 (MQ=255) tgggcctgggcGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGaa < 5:59302/56‑1 (MQ=255) | GACCTGATGCTGACGCTGCTGCCGTTCGAGGCGCGGTTCTACGAAGAGCAGGGTGTGCCAGTGCGTTTTGTCGGCCACCCGCTGGCAGACACCATACCGCTTGAAGCCGACCGCCCGG‑CCGCGCGTGCCGCGCTGGGCCTGGGCGAGGGGCCGGTGGTGGCGTTGATGCCGGGCAGCCGGGGTGGCGAAGTAGGGCGCTTGGGCGCGCTGTTCCTCGACGCCGCCGAGCGCCTCAGCCAACAAGTGCCGGGTGTGCGCTTCGTGTTG > ACB122/1791545‑1791811 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |