Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,941,183 | 0 | T | G | 55.0% | 15.5 / 24.2 | 20 | T49P (ACC→CCC) | hdfR | HTH‑type transcriptional regulator HdfR |
Reads supporting (aligned to +/- strand): ref base T (0/9); new base G (11/0); total (11/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.95e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACA‑CCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGAC > NZ_CP009273/3941098‑3941272 | gcCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGaaa < 1:556510/90‑1 (MQ=255) ccTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGaaaa < 2:595884/90‑1 (MQ=255) ccTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGaaaa < 1:617048/90‑1 (MQ=255) tGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGGGAAAAGGttacac > 2:207900/1‑86 (MQ=255) tGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCaca‑c < 2:156611/90‑1 (MQ=255) cATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGGGAAAAGGTT‑ACA‑ccccg > 2:172724/1‑88 (MQ=255) cATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGAAAAAGGTTACA‑ccccg > 1:544582/1‑88 (MQ=255) tGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACA‑CCCAGtt < 1:347948/90‑1 (MQ=255) gCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGAAAATGTTCACA‑ACCCCTTGAt > 1:371336/1‑90 (MQ=255) gTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGGAAAAGGGCACACCCCAGTTGttt > 2:586017/1‑90 (MQ=255) gCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACA‑CCCAGTTGATTTTCCAGt < 2:536978/90‑1 (MQ=255) aGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGGGAAAAGGTT‑ACA‑CCCCGTTGATTTTTCTGtttttt > 2:585903/1‑86 (MQ=255) aGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGGAAAAGGTCACA‑CCCCGGTGGTTTTTCCGTTGttt > 1:578530/1‑88 (MQ=255) aGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGGAAAAGGTCACA‑CCCCGGTGATTTTTCTGtttttt > 2:327900/1‑86 (MQ=255) aGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGAAAAGGTTCACA‑CCCCGTTGATTTTTCAGTTGttt > 1:226009/1‑88 (MQ=255) gtagtTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGGAAAAGTTCACA‑CCCAGTTGATTTTTCAGTTGTTTGAttt > 1:18652/1‑89 (MQ=255) gCGGTTAAACGGATATTGTTTCTGTGGA‑GGGTGAAAAGGTTCACA‑CCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGAc < 2:125208/90‑1 (MQ=255) ggTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACA‑CCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTg < 1:267916/90‑1 (MQ=255) ggATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACA‑CCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGAt < 2:611767/90‑1 (MQ=255) tatTGTTTCTGTGGC‑GGGGGAAAAAGGTCACA‑CCCCGTTGGTTTTTCACGTGTTTTTTTCTAAAACTCACTGCTGACTGGGGCAGATaga > 1:468306/1‑90 (MQ=255) tGAAAAGGTTCACA‑CCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGAc < 1:172724/90‑1 (MQ=255) | GCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACA‑CCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGAC > NZ_CP009273/3941098‑3941272 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |