Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,955 | 0 | T | G | 60.0% | 21.2 / 28.8 | 30 | T495P (ACC→CCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (0/12); new base G (18/0); total (18/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.16e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCC > NZ_CP009273/822868‑823035 | ccGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGa < 1:1125920/90‑1 (MQ=255) cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑Ggggggggg > 2:205831/1‑90 (MQ=255) aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑Ggggggggggggg > 1:1568198/1‑88 (MQ=255) tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCa < 1:1048690/90‑1 (MQ=255) cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAAc < 1:937580/90‑1 (MQ=255) tgggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAAcccccgga > 1:470546/1‑86 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGGGACGCCGGGAg > 2:937580/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGGGACGCCGGGAg > 2:776819/1‑90 (MQ=255) ggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAAAGCCGGGAGGCgg > 1:610638/1‑90 (MQ=255) gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTtc < 1:742796/90‑1 (MQ=255) gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTtc < 1:79760/90‑1 (MQ=255) cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAACGCCGGAAGGCGGgtcgtc > 1:891919/1‑90 (MQ=255) gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGaaaa > 1:1143546/1‑89 (MQ=255) ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc < 2:35153/90‑1 (MQ=255) tACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc < 1:476915/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGTGGGGGGGGACGCCGGAGGGCGGGTTGGCGGGAAAAAAAAAGTCCCGGTCAttc > 2:1276375/1‑88 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGGAACGCCCGAAGGCGGGTTGTCCAGAAAAAGAAAGTGCGGGTGATgg > 2:235850/1‑89 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAACGCCGGGAGTGGGGTCGTCGAGAAAAAGAATGTGCGGGGCATGc > 2:1027209/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAACGCCGGGAGGCGGGTTGTCGAGGAAAAAAAAATCCGGGTTATGc > 1:148693/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGAGGGGGGGGACCCCGGGAAGCGGGTCGGCGAGGAAAAGAAAGTTCGGGTCATGc > 2:60360/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGc > 1:32911/1‑90 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg < 1:60360/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGGAACGCCGGGAGGCGGGTCGTCGAGAAAAAAAATGTCCGGTTCATGCATccg > 2:1490583/1‑88 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAACAGAAAGTCCGGGTCATGCCTCAg > 1:683374/1‑90 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAAAGCCGGGAGGCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCCTGCATCAg > 1:874603/1‑90 (MQ=255) tgtgCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:1597129/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 2:514497/90‑1 (MQ=255) aGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGGCACGCCCGAAGGCGGGTCGTCGGGAAAAAAAAAGGCCGGGGCATGCATCAGCGAACAGgcc > 2:436753/1‑90 (MQ=255) aGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAAAGCCGGAAGGCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAAgcc > 1:166591/1‑90 (MQ=255) aGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAAAAATTCCGGGTTATGCATCAGCGAAAAGgcc > 1:402090/1‑90 (MQ=255) cACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGcc < 1:1276375/90‑1 (MQ=255) cACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGcc < 1:235850/90‑1 (MQ=255) | CCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCC > NZ_CP009273/822868‑823035 |
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: ‑ |