Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 395,288 | 0 | A | T | 44.0% | 9.2 / 33.7 | 25 | F334I (TTC→ATC) | ampH | peptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase |
Reads supporting (aligned to +/- strand): ref base A (10/4); new base T (11/0); total (21/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.05e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCGAACGAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCA > NC_000913/395221‑395354 | ggCGAACGAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGaa > 1:484113/1‑68 (MQ=255) ggCGAACGAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGaa > 1:2043470/1‑68 (MQ=255) aaCGAGTGCCCACCACAACCGCGCCAGTATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCg < 1:2975177/68‑1 (MQ=255) aGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGccac < 1:1231947/68‑1 (MQ=255) aGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGccac > 1:3253341/1‑68 (MQ=255) tGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGcgacca > 1:1917223/1‑68 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCg > 1:1055064/1‑68 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCg > 1:2549293/1‑68 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCg > 1:2341083/1‑68 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCg > 1:868137/1‑68 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGcg > 1:3006656/1‑67 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGcc > 1:1421073/1‑68 (MQ=255) cgcCAATATTGTTCGGCGGGATCATCGCCATATGGGTAATGAAACCGCCACCACCGCCTGTCTTCTGa < 1:3789297/68‑1 (MQ=255) cgcCAATATTCTTCTGGGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGa < 1:179137/68‑1 (MQ=255) cgcCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGa < 1:3173210/68‑1 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATcc > 1:3031139/1‑68 (MQ=255) ttctGCGGGGTCAGCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCgg < 1:1481193/68‑1 (MQ=255) ttctGCGGGATTATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCgg < 1:112047/68‑1 (MQ=255) gggATCATCGCCATATAGGTAATGATACGGCGACCACCGCCAGTTTCACGAATAACATCCGgacgatc > 1:3073755/1‑66 (MQ=255) atcatcGCAATATAGGTAATGAAACCGCCGCCACCGCCTGTCTTCTGAATAATCCCCGGACGACCctc < 1:872840/68‑1 (MQ=255) ccATATAGGTAATGATACGGCGACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCggg > 1:1536821/1‑68 (MQ=255) ccATATAGGTAATGAAACGGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCggg > 1:99357/1‑68 (MQ=255) ccATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCggg > 1:1912560/1‑68 (MQ=255) ataGGTATTGAACCCGCCCCCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa < 1:3652203/68‑1 (MQ=255) ataGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGctg > 1:1446110/1‑66 (MQ=255) ataGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGCGCCa > 1:3182536/1‑68 (MQ=255) tgggTAATGAACCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:2130849/66‑1 (MQ=255) taGGTAATAAACCCGCCACCCCCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:3378978/68‑1 (MQ=255) ggggTAATAAAGCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1930003/66‑1 (MQ=255) aaGGTAATGAACCCGCCACCCCCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1004304/67‑1 (MQ=255) aGGTAATGATACGGCGACCACCGCCGGCTATACGAATACTCCCCGGACGACCCTCCTTCGGGGCCaaa > 1:3080467/1‑66 (MQ=255) aGGAAATGAACCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCata < 1:412433/68‑1 (MQ=255) ggTAATGATACGGCGACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatat > 1:3214570/1‑68 (MQ=255) ggTAATGATACGGCGACCACCGAGATTCACCTGACTAACTCCCGGACGACCCTCTTTCGGGGACatat > 1:1179298/1‑68 (MQ=255) ggTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatat > 1:3051743/1‑68 (MQ=255) gTAATGATACGGCGACCACCGCCTGCCTCCTGCATCATCCCCGGACGACCCTCTTTCGGGTCCatata > 1:3283546/1‑68 (MQ=255) aaTGAAACCGGCTCCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCAtatata < 1:3592229/68‑1 (MQ=255) aaTGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatatatc > 1:2156154/1‑67 (MQ=255) aaaCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATCTATAcctg > 1:3475512/1‑66 (MQ=255) | GGCGAACGAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCA > NC_000913/395221‑395354 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |