Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 353,591 | 0 | T | C | 58.3% | 11.2 / 27.7 | 24 | K271E (AAA→GAA) | cynR | transcriptional regulator CynR |
Reads supporting (aligned to +/- strand): ref base T (8/2); new base C (0/14); total (8/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.12e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG > NZ_CP009273/353502‑353677 | cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 2:113666/1‑90 (MQ=255) cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 1:268515/1‑90 (MQ=255) attcattcCCGCCAACACCCGCGCTTTTATCCAACCCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGc‑a < 2:36411/90‑1 (MQ=18) ttCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGc‑aa > 1:351240/1‑90 (MQ=255) ccGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑cc > 1:410622/1‑90 (MQ=255) gCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑ccgc > 1:276781/1‑90 (MQ=255) aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGt > 2:339051/1‑90 (MQ=255) cacccGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGC‑AAAAAAA‑CCGCCGTtc < 1:132548/87‑1 (MQ=21) cacccGCGCATTTATCCAACGCAAGGTAAAAAAATGCTTTCGGGGGGGTTGTTTCCCAGTTTTTTTCCCCCCGC‑AACAAAA‑CCGCCGTtc < 1:286722/87‑1 (MQ=255) cgcgCATTTACCCAACCCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCATTTTTTTTCCCCCCGC‑AACAAAA‑CCGCCGTtctctc < 1:289763/90‑1 (MQ=21) aTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCagta > 2:50191/1‑90 (MQ=255) aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg > 1:72594/1‑90 (MQ=255) aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg > 1:83702/1‑90 (MQ=255) aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg > 2:91028/1‑90 (MQ=255) cAAAAATGCCTTCGCGGCGGCTGTCGCCCAGTTTTTTTCCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt < 1:182049/90‑1 (MQ=255) aaaaaaTGCCTTCGCGGGGGCTGTCTGCCAGTTTTTTTCCCCCCGA‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt < 1:445538/89‑1 (MQ=37) aaaaaTTCCTTCGGGGGCGCTTTCCCCCAGTTTTTTTCCCCCCCA‑AAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt < 2:344760/90‑1 (MQ=255) ctccgcgGGGGCTGTCTGCCATTTTTTTTCCCCCCGC‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg < 2:83702/87‑1 (MQ=255) ctccgcgGCGGCTGTCTGCCAGTTTTTTTCCCGCCGC‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg < 1:172524/87‑1 (MQ=255) tgggggcggtTGTCTGCCAGCTTTTTTCCCCCCGA‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 1:432829/86‑1 (MQ=255) tCGGGGGGGTTGTCTGCCAGTTTTTTTCCCGCCGA‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 1:163098/90‑1 (MQ=255) tCGCGGCGGCTTTCGGCCAGTTTTTTTCCCCCCGA‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 1:113666/90‑1 (MQ=255) tCGCGGCGGCTCTC‑CCCAACTTTTTTTCCCCCCCAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 1:488280/90‑1 (MQ=255) ccggggcggcTGTCCCCCATTTTTTTTCCCCCCAC‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 1:273476/86‑1 (MQ=255) ccggggcgccTGCCCGCCAGTTTTTTCCCCCCCCA‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 2:228850/86‑1 (MQ=255) ccgcgGCGGCTTTCTGCCACTTTTTTTCCCCCCAC‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 1:461498/89‑1 (MQ=255) ttttttCCGCCGC‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 2:470499/89‑1 (MQ=255) tttttcCCGCCGC‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 1:108675/89‑1 (MQ=255) tttttcCCCCCGC‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 1:326693/89‑1 (MQ=255) tttcCCGCCGA‑AAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 1:371829/90‑1 (MQ=255) | CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGC‑AACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG > NZ_CP009273/353502‑353677 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |