Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,963 | 0 | T | G | 52.9% | 10.4 / 37.4 | 34 | D492A (GAC→GCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (3/13); new base G (18/0); total (21/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.03e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT > NZ_CP009273/822876‑823052 | 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:79760/90‑1 (MQ=255) gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTtc < 1:742796/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) gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa < 2:683374/90‑1 (MQ=255) gggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCCGTTTAAAACCt > 1:1378237/1‑90 (MQ=255) gTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa > 1:494847/1‑83 (MQ=255) gTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa < 2:494847/83‑1 (MQ=255) gTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTaa > 2:1336975/1‑90 (MQ=255) tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAt > 2:1072616/1‑90 (MQ=255) tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAt > 1:1083845/1‑90 (MQ=255) | TCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT > NZ_CP009273/822876‑823052 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |