| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 69,900 | 0 | A | C | 55.6% | 31.2 / 42.3 | 36 | V370G (GTG→GGG) | thiP | thiamine/thiamine pyrophosphate ABC transporter permease ThiP |
| Reads supporting (aligned to +/- strand): ref base A (12/4); new base C (0/20); total (12/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.45e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATAC > NZ_CP009273/69815‑69982 | gATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTg > 2:40865/1‑90 (MQ=255) tAGTGTTGTGGAGAAGAAAAAAGAAGCCGGTAGCCAGCACAACCCCCGGCAGGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGAc < 2:174734/90‑1 (MQ=255) ttttgtggAGAAGAAAAAAGAAGCCGGTAGCCAGCAAAACCCCCGGCATGGCGAGGATCAACATGCC‑CTCATCCCCAGCCCCTGACCCGc < 2:92192/88‑1 (MQ=255) gttgttGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGc > 1:210741/1‑90 (MQ=255) tgAGCAAAAAAAAAAAGCCGGTACCCAGCACAACCCCCGGCATGGGGAGGATAAACATGCCGCTCATCTCCAGCCCCTGCCCCGCCAGCa < 2:237360/90‑1 (MQ=255) gCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCAttt > 2:221330/1‑90 (MQ=255) aagaagCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCcgc > 1:588778/1‑90 (MQ=255) gCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACGTGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCg > 2:138053/1‑90 (MQ=255) gCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCg > 1:366463/1‑90 (MQ=255) gCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCg > 1:386941/1‑90 (MQ=255) ccGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGc < 1:249316/90‑1 (MQ=255) ccGGTAGCCACCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGc < 1:5314/90‑1 (MQ=255) ccGGTACCCAGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGc < 2:207151/90‑1 (MQ=255) aGCCAGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTc < 1:175004/90‑1 (MQ=255) cacaACCCCCGGCATGGCGAGGACCAACATGCCCCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt < 1:102855/90‑1 (MQ=255) cacaACCCCCGGCAGGGAGAGGTCCAACATCCCGCTCAGCCCCAGCACCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACt < 2:231843/90‑1 (MQ=255) caATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGc > 1:237714/1‑90 (MQ=255) caATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGc > 2:216246/1‑90 (MQ=255) aTCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTc > 2:384445/1‑90 (MQ=255) cccccGGCATGGCGAGGATCAACTTGCCCCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTcc < 1:186193/89‑1 (MQ=255) ggCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagc < 1:221330/90‑1 (MQ=255) ggCAGGGCGGGGACCAACATGCCGCTCATCCCCACCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagc < 2:293283/90‑1 (MQ=255) ggagaATCAACATGCCGCTCACCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 2:497531/89‑1 (MQ=255) cGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 1:216246/90‑1 (MQ=255) aTCAACATGCCCCTCATCTCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTc < 1:21043/90‑1 (MQ=255) cAACATCCCCCTCCCCCCCACCCCCTCCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAg < 1:424977/90‑1 (MQ=255) cacgCCCCTCTTCCCCCGCCCCTCCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCAc < 2:210741/87‑1 (MQ=255) cccTCACCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 2:571096/88‑1 (MQ=255) cGCTCATCTCCAGCACCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 1:386266/90‑1 (MQ=255) cGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 2:216966/90‑1 (MQ=255) cGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 2:580052/90‑1 (MQ=255) cGCTCACCCCCACCCCCTCCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 1:466299/90‑1 (MQ=255) cGCCCATCTCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 1:287785/90‑1 (MQ=255) cATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACaa > 1:72497/1‑90 (MQ=255) ctcCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATac > 1:293949/1‑90 (MQ=255) ccccAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATac < 2:72497/88‑1 (MQ=255) | GATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATAC > NZ_CP009273/69815‑69982 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |