Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,965,782 | 0 | T | G | 74.2% | 63.6 / 19.6 | 31 | E389D (GAA→GAC) | tar | methyl‑accepting chemotaxis protein II |
Reads supporting (aligned to +/- strand): ref base T (5/3); new base G (6/17); total (11/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.47e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAA > NZ_CP009273/1965696‑1965864 | aGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCa > 2:123168/1‑90 (MQ=255) tttGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGc < 1:184007/90‑1 (MQ=255) ctctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCAcgcgc < 1:909089/90‑1 (MQ=255) tctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCAcgcgcg > 2:580992/1‑90 (MQ=255) tctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCAcgcgcg > 1:770518/1‑90 (MQ=255) gccgccTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCa > 2:710600/1‑90 (MQ=255) gccgccTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCa > 1:28005/1‑90 (MQ=255) cgccTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAAc < 1:1034880/90‑1 (MQ=255) cgccTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAAc < 1:433523/90‑1 (MQ=255) gccTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAAcc > 2:966789/1‑90 (MQ=255) cTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACcgc > 2:28813/1‑90 (MQ=255) gggCGCTGCGACTGGCAAGATTCCGCCCTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCgg < 2:528043/90‑1 (MQ=25) cTGCGACTGGAAAGATTCCGCCCTTCCCCCGCCACCACGGCAAAACCCCGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTc < 2:1054294/90‑1 (MQ=255) gccgggCAAGATCCCGCCCTCCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAgcg < 2:599431/86‑1 (MQ=25) cTGGCAAGATTACGCCCTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg < 2:1109424/90‑1 (MQ=255) cTGGCAAGATTACGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg < 1:544966/90‑1 (MQ=255) ggCAAGATTCCGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAg < 2:896320/90‑1 (MQ=37) ggCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAg < 1:710600/90‑1 (MQ=255) ggCAAGATTACCCCCTCCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAg < 2:129512/90‑1 (MQ=37) gATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGatat < 1:1065471/90‑1 (MQ=255) ccTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTgg < 1:1141910/89‑1 (MQ=255) ttCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCCCGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGCTATTAGTCTGGaa > 2:1091075/1‑90 (MQ=255) ttCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGaa > 1:317263/1‑90 (MQ=255) tCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAg > 1:1086759/1‑90 (MQ=255) ccGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCaa < 1:123168/90‑1 (MQ=255) gCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATa < 2:807765/90‑1 (MQ=255) gCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATa < 1:600170/90‑1 (MQ=255) gCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCAAGGATATTAGTCTGGAAGGCAATACCGTCGATa < 2:24412/90‑1 (MQ=255) aaaaCCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAc < 2:533114/90‑1 (MQ=255) aaaCCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACg < 2:1086759/90‑1 (MQ=255) ggCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATaa > 1:848798/1‑90 (MQ=255) | AGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAA > NZ_CP009273/1965696‑1965864 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |