Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,062,357 | 0 | T | A | 24.0% | 38.7 / 10.7 | 25 | intergenic (+22/+34) | asnV/ldtA | tRNA‑Asn/L,D‑transpeptidase linking Lpp to murein |
Reads supporting (aligned to +/- strand): ref base T (11/8); new base A (2/4); total (13/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.78e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGC > NC_000913/2062219‑2062498 | tCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCttt < 3:95373/149‑1 (MQ=255) aCCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCtttt < 4:114898/148‑1 (MQ=255) ccATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCtttttt < 1:73128/149‑2 (MQ=255) ccATGCACTGCAATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCTGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGATCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCtttttt < 1:70678/149‑2 (MQ=17) cATGCACTGCCATTACTATGCCCCCCTTTCACAAGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGGCTGTTAATCCGGATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGc > 4:106902/1‑149 (MQ=18) cATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGc > 3:17974/1‑149 (MQ=255) cACTGCCATTAATATGCGCCCCGTTCACACGAGTCCTCTGTAGTTCAGTCGGTAGAAGGTCGGACTGGTAATCCGTATGTCACTGGGTCGAGTCCAGTCAGAGCAGCCAAATGCAAAAAAGCCTGCTATCTAGCAGGCTTTTTGCttt < 3:106902/148‑1 (MQ=18) gATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCa > 4:72778/1‑105 (MQ=255) tAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGaa < 4:17974/149‑1 (MQ=255) gTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAAc > 1:51490/1‑148 (MQ=255) cAGTCGGTAGAACGGCTGACTGTTAATCCTTATGTCACTGGTTCGAGTACAGTCAGAGTAGCCAAATTCAAAAAAGCCTGATTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGATCTTAAAACATCTGTCTTTAACAAGAACTAAt > 8:78220/1‑149 (MQ=255) gAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAgg > 6:14581/1‑149 (MQ=255) aaggtttCGAGTCCAGTCAGATGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgcagc < 7:78220/144‑1 (MQ=255) aCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTACTTACCAACGCTCTTAACACATCTGTCTAGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgcag > 4:74508/1‑148 (MQ=255) aCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgcag > 6:157583/1‑148 (MQ=255) aGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTa > 3:125078/1‑149 (MQ=255) aGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTa > 4:96351/1‑149 (MQ=255) agagGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACAt < 5:157583/148‑1 (MQ=255) ggagCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAAc < 5:14581/149‑1 (MQ=255) agcaaaaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACa < 4:105406/131‑1 (MQ=255) agcaaaaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACa > 3:105406/4‑134 (MQ=255) aaaaaaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGtt < 4:125078/149‑1 (MQ=255) aaaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTAc < 7:72533/129‑1 (MQ=255) aaaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTAc > 8:72533/1‑129 (MQ=255) tGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGc > 1:74037/1‑149 (MQ=255) | TCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGC > NC_000913/2062219‑2062498 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |