Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,955 | 0 | T | G | 51.2% | 24.7 / 37.0 | 41 | T495P (ACC→CCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (3/17); new base G (21/0); total (24/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.52e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAA > NZ_CP009273/822866‑823044 | ttCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGAGt > 1:1510419/1‑90 (MQ=255) cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGgggggggg > 2:816500/1‑90 (MQ=255) cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAggg < 1:1017514/90‑1 (MQ=255) cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAggg < 1:997703/90‑1 (MQ=255) aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGgggggggggggg > 2:1283168/1‑88 (MQ=255) aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGgggggggggggg > 1:1896513/1‑88 (MQ=255) aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAggggggg > 1:1139225/1‑88 (MQ=255) aGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGACgg < 1:791953/90‑1 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCCACGACGGAAg > 2:750589/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCGACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAg > 2:1239557/1‑90 (MQ=255) gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAGGCggggc > 2:1401034/1‑87 (MQ=255) tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGGAAGCCGGGAGGCGGGgcgtcgg > 1:534664/1‑89 (MQ=255) ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCGGGTCGTCGAGaa > 2:320666/1‑90 (MQ=255) gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAg < 1:836225/90‑1 (MQ=255) gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCCGGAGGCGGGTCGTCGAGaaaaaa > 1:159796/1‑87 (MQ=255) cTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGGAAAAGAAAGt > 2:124702/1‑90 (MQ=255) ctctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGAGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCgg < 2:1510419/90‑1 (MQ=255) cATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCGAAGCCCGAAGGCGGGTCGTCGGGAAAAAAAAAGTACCGGTTctg > 1:228042/1‑88 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGc < 2:228042/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGc < 1:2135301/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGc < 1:649455/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGAAGTCGGGTCGTCGAGAAAAAAAAAGT‑CCCGGTCATGc > 2:107278/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGGCGGGTCGTCGAGAAAAAGAAAGT‑CCGGGTCATgg > 1:1977453/1‑89 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGAGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGc > 1:1002589/1‑90 (MQ=255) tGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGGTCGTCGAGAAAAAGAAAGT‑CCGGGTCATGCa > 2:88372/1‑90 (MQ=255) gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGGAGGCGGGTCGTCGGGAAAAAGAAAGT‑TCGGGTTATgggt > 1:1559549/1‑87 (MQ=255) cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAATCGGGTTGTCGGGGAAAAGAAAGT‑TCGGGTGATGCaac > 2:1719171/1‑88 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAGGCGGGTCGTCGGGAAAAAAAAAGT‑CCGGGTCATGCATCAg > 2:791548/1‑90 (MQ=255) tGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGCATCAGCg < 1:534357/90‑1 (MQ=255) tGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGT‑CCGGGTCTTGCATCAGCg > 2:135991/1‑90 (MQ=255) tgttCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGCATCAGCGAAc < 2:128428/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGCATCAGCGAAc < 2:1828000/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGCATCAGCGAAc < 2:1978776/90‑1 (MQ=255) gCAGCCAAAATTCACGTCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGCATCAGCGAACAgg > 1:1209667/1‑90 (MQ=255) gCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTGGGGTCTTCGTGAAAAAAAAAGG‑CCGGCTCCTGCCACAGCGGACaag > 1:1828000/1‑88 (MQ=255) gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGGCGGGAAAAAGAAAGG‑CCGGGTCCTGCATCAGCGAACaag > 1:701814/1‑88 (MQ=255) gCAGCCAAAATTCACGGCGGGGGAGGGGGGGAACGCCGGAAGTCGGGTCGTCGAGAAACAAAAAGT‑CCGGGTCATGCATCAGCGGACaag > 1:426024/1‑88 (MQ=255) aGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGGTCGTCGAGAAAAAAAAAGT‑CCGGGTCATGCATCAGCGGACAgggc > 1:548701/1‑88 (MQ=255) aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAGGCGGGTCGTCGAGAAAAAGAAAGG‑TCGGGTCATGCCTCCGCGGACAGgcc > 1:1182850/1‑90 (MQ=255) cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 1:169800/90‑1 (MQ=255) cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 1:1283168/90‑1 (MQ=255) cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 1:471646/90‑1 (MQ=255) cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGATCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 2:534664/90‑1 (MQ=255) cggGGGGGGGGGGCAACCCCGGGAGTCGGGTTGTCGAGGAAAAGAAAGT‑CCGGGTCATGCATCAGCGAACAGGCCACCGCCCGCCGCCGt > 2:1018213/1‑90 (MQ=255) ggTGAGGGGGTCAACGACGGAAGTAGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGttt < 1:750589/90‑1 (MQ=255) tGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa < 2:1948968/90‑1 (MQ=255) | TTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGT‑CCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAA > NZ_CP009273/822866‑823044 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |