Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,145,089 | 0 | G | A | 43.5% | 19.6 / 10.1 | 23 | E142K (GAA→AAA) | yegJ | DUF2314 domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base G (4/9); new base A (8/2); total (12/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.61e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTA‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTAT > NZ_CP009273/2144963‑2145226 | aCAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑AACAATTATTCTTTATTTGTTACCTA‑TCCAAAAAAAGAAACcgac > 2:696671/1‑138 (MQ=255) gAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGTGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTA‑TCCAAAAGAAGAAGCCGACGAGCTAagag < 2:576048/139‑1 (MQ=255) gAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTA‑TCCAAAAGAAGAAGCCGACGAGCTAagag < 1:450401/139‑1 (MQ=255) ggAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑AACATTTATCCTTTTTTTTGTACCTATTCCAAAAAAAAAAG‑CAACCAGCTCAAAGCAAAAATAGtt > 1:870883/1‑139 (MQ=255) ggAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTTATTCTTTTTTTGTCACCCA‑TCCAAAAAAAGAAACCGACGAGCCAAGAGCAAAAATAAtt > 1:699684/1‑139 (MQ=255) ggAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTATAACCATT‑ATTCTTTTTTTGTGTCCCA‑ACCAAAAAAAAAAGCCGGCGAGGTTAGAGAAAAAATAAtt > 1:523424/1‑139 (MQ=255) gAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTA‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTaa < 1:47590/139‑1 (MQ=255) gAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑AACATTTTTTCTTTTTTTGTCACCTT‑TCCAAAAAAAGAAGCCGACCGGCTAAGAGCAAAAATAATTTGCtta > 2:70864/1‑137 (MQ=255) aTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTA‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAg < 1:696671/139‑1 (MQ=255) agaaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTA‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGt < 2:222353/139‑1 (MQ=255) aTGTATACGCGCGGGGATAAAATTTACGGCGGTTTT‑TCCATTTTTCCTTTTTTTGTGACCCA‑TCCAAAAAAAGAAGCCGACCAGGTACAAAGAAAATTAATTTGCTTATAAACCGTTTTATTTTTTTTTTTAACTGGCg > 2:759164/1‑139 (MQ=255) tGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTA‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCgg < 1:880505/139‑1 (MQ=255) cgcgGGGATAAAATTTACGGCGGTTAT‑ACCATTTATCCTTTTTTTGTTACCTA‑TTCCAAAAAAAAAACCGACCAGCTAAAAGGAAAAATAATTTGCTTAATAGCCCGTTTATTTTTTTTTTTAAATGGGGGGTgggttt > 1:670466/1‑134 (MQ=255) cgGGGATAAAATTTACGGCGGTTTT‑ACCATTTTTTCTTTTTTTGTCACCCA‑TCCAAAAAAAAAAACCGGCGAGCTAAGAGCAAAAATAGTTTGCTAATTAGCCGTTATTTTTTTTTTTTTAAATTGCGGGTGGAttttt > 1:201755/1‑139 (MQ=255) gcgtTTAT‑CCAATTGACCCTTTATTGGCCCCCTA‑CCCAAAAGAAGAAGCCGACGAGCAAAGAGCAAAATTAGTCCGCTAACCAGCCTTTATATTTATTGTTAAAACGGGCGGTTAATTTGTTGTCTTCAACCGCCATga < 2:728168/139‑1 (MQ=255) tat‑cCCATTGATCCTTTATTGTTCACCTA‑TCCAAAAAAAGAACCCGACGATCAAAGAGCAAATTTGGTTCTCTAACCAGCCTTTATATTTATTGTTTAACCTGGGGTTTAAATTGTGGTCTTCAACCGCCATGAGaaa < 1:495916/138‑1 (MQ=255) ttGATCCTTTATTGGTCACCTA‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGGATTGTTGTCTTCAACCGCCATGGGAAACTCActct > 2:213825/1‑139 (MQ=255) tGATCCTTTATTGGTCACCTA‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCCTGAGAAACTCACTCTg > 1:894596/1‑139 (MQ=255) tGATCCTTTATTGGTCACCTA‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCACGAGAAACTCCCTCTg > 1:849211/1‑139 (MQ=255) ggTCCCCTA‑CCAAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTATTCCGCTAACCAGCCGTTATATTTATTGTTAAACCTGGGGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGcc < 2:670466/138‑1 (MQ=255) aCCTA‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGc > 1:632327/1‑139 (MQ=255) cttA‑TCAAAAAAAAAAAGCCTCCGACTAAATATCAAATTTAGTGCGCTATCCGCCCTTTATTTTTTTTGTTAAACCGGGGGATTGATTTTTTGTTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGcc < 2:849211/137‑1 (MQ=255) aaaaaGAAGAAGCCGACGAGCTAAGAGCAAAATTGTTTCGCAAATCACCGTTTATATTTATTGTTTAAACTGGCGTTTGAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGtt < 1:386524/138‑1 (MQ=255) gaagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTCAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTat > 1:695961/1‑138 (MQ=255) | ACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAT‑ACCATTGATCCTTTATTGGTCACCTA‑TCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTAT > NZ_CP009273/2144963‑2145226 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |