Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 353,601 | 0 | C | A | 59.3% | ‑4.9 / 22.0 | 27 | L267F (TTG→TTT) | cynR | transcriptional regulator CynR |
Reads supporting (aligned to +/- strand): ref base C (9/2); new base A (0/16); total (9/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.17e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG > NZ_CP009273/353514‑353677 | aTTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑‑CCGCCGcaacaa > 1:117104/1‑90 (MQ=255) tCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑‑CCGCCGCAACAAAA‑c > 2:515193/1‑90 (MQ=255) acaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTCCCAGCTATTTTT‑‑CCGCCGCAACAAAA‑CCGCCGtt > 1:481757/1‑90 (MQ=255) cacccGCGCTTTTCCCCAACCCCTTGTAAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTC‑‑CCCCCGAAAAAAAA‑CCGCCGTtc < 1:309093/87‑1 (MQ=12) cacccGCGCTTTTACCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTC‑‑CCGCCGCAACAAAA‑CCGCCGTtc < 2:154200/87‑1 (MQ=25) cacccGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGCTTTTTTC‑‑CCGCCGCAAAAAAA‑CCGCCGTtc < 2:146748/87‑1 (MQ=21) cgcgCTTTTACCCACCCCCATGTCAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGCTTTTTCC‑‑CCGCCGCAAAAAAA‑CCGCCGTtctctc < 2:181815/90‑1 (MQ=19) ccgcTTTTACCCACCCCCATGTGAAAAAATCCTTCCGGGGGGGTTGTCTCCCAGTTTTTTTC‑‑CCCCCGAAAAAAAA‑CCGCCGTTCTCTcc < 2:434435/89‑1 (MQ=255) cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCa > 2:409890/1‑90 (MQ=255) tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTg > 2:374280/1‑90 (MQ=255) tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTg > 2:50568/1‑90 (MQ=255) aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg > 1:462313/1‑90 (MQ=255) aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg > 1:282907/1‑90 (MQ=255) aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg > 2:535079/1‑90 (MQ=255) aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgctg > 1:244944/1‑87 (MQ=255) gCCTTGTGAAAAAATGCCTGCGCGGCGGCTGTCTGCCAGTTTTTTTT‑‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCa < 2:448543/90‑1 (MQ=37) cccATGTGAAAAAATGCCTTCGCGGGGGCTGTCGGCCAGCTTTTTTC‑‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCa < 2:217010/89‑1 (MQ=255) caggtgAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGCTTTTTTC‑‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAag < 2:38254/87‑1 (MQ=37) aaaaaaTCCCTCCGGGGGGGCTTTCCCCCAGTTTTTTTC‑‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt < 2:370102/89‑1 (MQ=255) aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTC‑‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa < 2:68689/90‑1 (MQ=38) tccgcgGCGGCT‑TCATCCAG‑TATTTTTCCCCCCCGCAAC‑AAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 1:390488/88‑1 (MQ=255) ggggcTGTCTGCCAGTTTTTTTC‑‑CCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa < 2:586366/88‑1 (MQ=255) gcggtTGTCGGCCATTTTTTTTC‑‑CCGCCAAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa < 1:139378/90‑1 (MQ=255) gcTGTCCCCCCGTTATTTTC‑‑CCCCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt < 2:367800/90‑1 (MQ=255) tttttc‑‑CCCCCACAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 2:351433/89‑1 (MQ=255) ttttt‑‑CCCCCCAAAAAAAACCCGCCGCTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 1:75147/90‑1 (MQ=255) tttc‑‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 2:46727/90‑1 (MQ=255) tttc‑‑CCCCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 2:608272/90‑1 (MQ=255) | ATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG > NZ_CP009273/353514‑353677 |
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: ‑ |