Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,475,282 | 0 | T | G | 55.0% | 12.7 / 19.1 | 20 | T163P (ACC→CCC) | pepA | leucyl aminopeptidase |
Reads supporting (aligned to +/- strand): ref base T (2/7); new base G (11/0); total (13/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.64e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGATCTTTTGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGC > NZ_CP009273/4475199‑4475343 | gagaTCTTTTGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCAc > 2:434243/1‑90 (MQ=255) tgctgcTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAacac < 1:434243/90‑1 (MQ=255) tgctgcTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAacac < 2:334770/90‑1 (MQ=255) gcggcAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGga < 1:521169/90‑1 (MQ=255) aTCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGAcgc < 1:116391/90‑1 (MQ=255) aGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGTGGCTGCACCTTGAACACCATNTTACGCAACGGACGACGCgggttg > 2:179646/1‑86 (MQ=255) cGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCCCGGTGAAAACCCCCCTTCCCCGCGGACGACGCGGGTCGccct > 1:419059/1‑90 (MQ=255) gTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCa < 2:421003/47‑1 (MQ=255) gTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCa > 1:421003/1‑47 (MQ=255) gTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGAACACCATCTTTCGCAGCGGGCGACGCGGGTCGCTCtt > 2:181700/1‑90 (MQ=255) gTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGAAAACCAACTTACGCAGCGGGCGACGCGGGGCGCTCtt > 2:295807/1‑90 (MQ=255) cgcgcgCTCACCGCTGGTCAGTTCACGGCGGGGCGGGCCGGTGTACACCCCCCTTCCCCGCGGGGGACGCGGGGCGCTCCTGTTCGtttt > 1:344254/1‑90 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAacac < 1:139036/47‑1 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAacac > 2:139036/1‑47 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGGTCGCTCTTGTTCGttttt > 2:159379/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGGACGTTGAACACCATCTTTCGCAGCGGGCGACGCGGGTTGCTCTTTTTTGttttt > 2:419540/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCCCGTTGGAAACCCCCTTTCCCCGGGGGCGGCGCGGGTGGCCCTTTTTTGttttt > 1:371773/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGAACACCATCTTTCCCAGCGGACGACGCGGGTTGCTCTTTTTTGttttt > 2:318984/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACCGTGAAAACCCTCTTACGCAGCGGGCGAGGCGGGTTGCTCTTTTTTGttttt > 2:258155/1‑89 (MQ=255) cgcTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGc < 2:59888/90‑1 (MQ=255) | GAGATCTTTTGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGC > NZ_CP009273/4475199‑4475343 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |