Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 70,068 | 0 | A | C | 50.0% | 41.4 / 13.7 | 36 | V314G (GTA→GGA) | thiP | thiamine/thiamine pyrophosphate ABC transporter permease ThiP |
Reads supporting (aligned to +/- strand): ref base A (17/1); new base C (0/18); total (17/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.19e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGA > NZ_CP009273/69987‑70157 | gccgccAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATcg > 2:905710/1‑90 (MQ=255) gccgccAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATcg > 2:162977/1‑90 (MQ=255) cgccAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATcgac > 2:791308/1‑90 (MQ=255) ccGCAACGAGGCCCACGGCCCCTGCCACAGCCCCGGTTGTGCCAGCCTTCCCGGC‑AACTGGCGATTTCCCCCATCGACGATCACCGCCAg < 1:598159/89‑1 (MQ=255) aGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTg > 2:244037/1‑90 (MQ=255) gcCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTCCCCCATCGACGATCACCg < 1:981470/67‑1 (MQ=255) gcCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCg > 2:981470/1‑67 (MQ=255) ccTGCCACAGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAcagca > 1:1210081/1‑90 (MQ=255) ccTGCCACAGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAcagca > 1:972128/1‑90 (MQ=255) tGCCACAGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAcagcagc > 2:903108/1‑90 (MQ=255) tGCCACAGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAcagcagc > 1:1072315/1‑90 (MQ=255) cccACAGCACCGGTTGTGCCAGCCCTTCCGGC‑AACTGGCGTTTTCCCCCCTCGACGATCACCGCCAGTAACGGTGGCAGCAAcagcagca < 2:1121252/89‑1 (MQ=255) caGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTAc > 2:927139/1‑42 (MQ=255) caGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTAc < 1:927139/42‑1 (MQ=255) caGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAgcgc > 2:73394/1‑90 (MQ=255) ggggTTGCGCCACCCCTTCCGGC‑ACCGGGCGTTTCCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGca < 2:1008440/88‑1 (MQ=255) cgggTTGTGCCAGCACTCCGGGC‑AACGGGCGATTTCCCCCACCGACGATCCCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGca < 2:4739/88‑1 (MQ=255) cgggTTGTGCCAGCACTCCCGGC‑ACCTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGca < 1:1136073/88‑1 (MQ=255) cgggTTGTCCCAGC‑CCTCCGGCAAACGGGCGTTTTCCCCCATCGACAATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGca < 2:186289/88‑1 (MQ=255) cgggTTGCGCCCGCTCCCCCGGC‑AACTGGCGTTTTCCCCCATCGACGACCCCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGca < 2:616669/88‑1 (MQ=255) cgggTTGCGCCAGCACTTCCGGC‑AACGGGCGTTTTCCCCCACCGACAATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGca < 2:575212/88‑1 (MQ=255) ccGGTTGTGCCAGCCCTCCCGGC‑AACTGGCGTTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGca < 1:162977/90‑1 (MQ=255) ccGGTTGTGCCAGCCCTCCCGGA‑AACGGGCTTTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGca < 1:791308/90‑1 (MQ=255) ccGGTTGGGCCAGCTCTCCCGGA‑AACTGGCGTTTCCCCCCACCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGca < 2:454048/90‑1 (MQ=255) ccGGTTGCGCCACCCCTCCCGGC‑AACTGGCGTTTTCCCCCACCGACAATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGca < 2:167055/90‑1 (MQ=255) ggTTGTGCCAGCATTCCCGGC‑AACTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGcaca < 1:220271/90‑1 (MQ=255) ggTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGcaca > 2:619372/1‑90 (MQ=255) tgCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTaa > 2:462684/1‑90 (MQ=255) tgCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTaa > 1:4739/1‑90 (MQ=255) gCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACAcct > 1:1008440/1‑89 (MQ=255) ttCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTc > 2:941050/1‑90 (MQ=255) c‑AACTGGCGTTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAAt < 1:905710/90‑1 (MQ=255) c‑AACTGGCGTTTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAAt < 2:117009/90‑1 (MQ=255) c‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAAt > 2:456609/1‑90 (MQ=255) acgggCGTTTTCCCCCATCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGc < 1:803749/87‑1 (MQ=255) ggCGATCTCCCCCACCGACAATCCCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGc < 2:329631/90‑1 (MQ=255) aCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGcc > 2:1036756/1‑51 (MQ=255) aCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGa > 1:1026486/1‑90 (MQ=255) | GCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGC‑AACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGA > NZ_CP009273/69987‑70157 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |