| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 353,591 | 0 | T | C | 52.9% | 11.5 / 21.4 | 17 | K271E (AAA→GAA) | cynR | transcriptional regulator CynR |
| Reads supporting (aligned to +/- strand): ref base T (6/2); new base C (0/9); total (6/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.26e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGC > NZ_CP009273/353507‑353675 | gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc < 2:371734/90‑1 (MQ=255) gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc > 2:213498/1‑90 (MQ=255) aacaccCGCCCTTTTACCCACCCCCATGTGAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGC‑AACAAAA‑CCGCCGt < 2:33287/90‑1 (MQ=21) cacccGCGCATTTACCCAACGCCTTGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGC‑AAAAAAA‑CCGCCGTtc < 2:708693/87‑1 (MQ=21) cacccGCCCTTTTACCCACCCCCATGTGAAAAAATCCCTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGC‑AAAAAAA‑CCGCCGTtc < 2:555807/87‑1 (MQ=255) aTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTgg > 2:242751/1‑90 (MQ=255) aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg > 1:157431/1‑90 (MQ=255) aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AAC‑AAA‑CCGCCGTTCTCTCCAGTAGTGGCggggg > 2:1068423/1‑89 (MQ=255) caaaaaaTGCCTTCGCGGCGGCTTTCGGCCAGTTTTTTTCCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa < 2:439323/88‑1 (MQ=255) aaaaaacTCCCTTCGCGGCGGCTTTCGGCCA‑CTTTTTTTCCCCCCCAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa < 2:1079678/88‑1 (MQ=255) aaaaaaTCCCTCCGGGGGGGTTGTCTGCCAGTTTTTTTCCCCCCGC‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt < 2:564061/89‑1 (MQ=255) gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGTtc > 1:307782/1‑55 (MQ=255) gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAACtttt > 1:117770/1‑89 (MQ=255) gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGC‑AACAAAA‑CCGCCGTtc < 2:307782/55‑1 (MQ=255) cgtcTGTCTCCCACTTTTTTTCCCCCCCC‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAt < 2:531474/87‑1 (MQ=255) gcTGTCTCCCATTTTTTTTCCCCCCGA‑AAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt < 2:1200112/90‑1 (MQ=255) gcTGTCGGCCAGCTTTTTTCCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt < 2:483126/90‑1 (MQ=255) tttttcCCCCCGA‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 2:485568/89‑1 (MQ=255) tttttCCCCCCA‑AAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 1:186758/90‑1 (MQ=255) | GTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGC > NZ_CP009273/353507‑353675 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |