Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 908,794 | 0 | A | G | 74.5% | 87.5 / 29.0 | 51 | G167G (GGT→GGC) | ltaE | low‑specificity L‑threonine aldolase |
Reads supporting (aligned to +/- strand): ref base A (6/7); new base G (20/18); total (26/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.55e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCG > NC_000913/908656‑908926 | gTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCAc < 2:437103/140‑1 (MQ=255) gTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCAc < 1:942376/140‑1 (MQ=255) aTCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGt < 1:842315/140‑1 (MQ=255) tCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTc > 1:1006246/1‑140 (MQ=255) tCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTc > 1:274831/1‑140 (MQ=255) ttACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATg < 2:696674/140‑1 (MQ=255) ttACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATg < 2:274831/140‑1 (MQ=255) ttACCGAAGAGAAATGCACCGACAGGCGTGCCAAGACCTTTCGACAGGGAAATGGTGAACGAATCACAATATTGTGTGATCTCTTTAAGTTCGCAGCCGTAAGCCGCCACGGCATTAAAGATGCGGGCGCCATCAACATg < 2:753511/140‑1 (MQ=255) aaCCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTgcgc > 1:286764/1‑140 (MQ=255) aCTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTcgcg > 2:804763/1‑140 (MQ=255) aCTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTcgcg > 2:557284/1‑140 (MQ=255) aCTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTcgcg > 2:339703/1‑140 (MQ=255) aCTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTcgcg > 1:709342/1‑140 (MQ=255) aCTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTcgcg > 1:802434/1‑140 (MQ=255) cAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTcc > 2:509611/1‑140 (MQ=255) cTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGct < 1:1118393/140‑1 (MQ=255) tttCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGctt < 2:802434/140‑1 (MQ=255) ttCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGcttc > 1:1005949/1‑140 (MQ=255) tCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGcttct < 2:1090808/140‑1 (MQ=255) cGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGcttctt < 1:509611/140‑1 (MQ=255) ggCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAgg < 1:557284/140‑1 (MQ=255) gCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGt > 1:688970/1‑140 (MQ=255) tGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTcc < 2:1244197/140‑1 (MQ=255) gAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCcgcg > 1:281783/1‑140 (MQ=255) aCGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGc < 1:886932/140‑1 (MQ=255) aCGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGc < 2:524557/140‑1 (MQ=255) aTCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAaca > 1:1150115/1‑140 (MQ=255) ttGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGt > 1:664229/1‑140 (MQ=255) tGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGtt > 1:248829/1‑140 (MQ=255) tGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGtt > 1:281755/1‑140 (MQ=255) tGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGtt > 1:331163/1‑140 (MQ=255) cGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTtgt > 1:142376/1‑140 (MQ=255) aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTGGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACCCTTTGCCGTTGTGGgt < 1:804763/140‑1 (MQ=255) aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt < 2:1005949/140‑1 (MQ=255) aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt < 1:1150662/140‑1 (MQ=255) aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt < 2:248829/140‑1 (MQ=255) aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt < 2:142376/140‑1 (MQ=255) aTCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt < 2:1006246/140‑1 (MQ=255) tctTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGtt > 1:155411/1‑140 (MQ=255) tttCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGtttt > 1:1393356/1‑140 (MQ=255) tCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt > 2:389534/1‑133 (MQ=255) tCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGgt < 1:389534/133‑1 (MQ=255) gTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCAGCAACACTTTGCCGTTGTGGGTGTTTTCCAGa < 2:281783/140‑1 (MQ=255) ttCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGAc > 1:1145981/1‑140 (MQ=255) tCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACGCTTTGCCGTTGTGGGTGTTTTCCAGACt < 2:274174/140‑1 (MQ=255) gCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGa < 1:416665/140‑1 (MQ=255) cAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAg > 1:1227268/1‑140 (MQ=255) cAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAg > 2:533615/1‑140 (MQ=255) ccacGGCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCgg < 2:709342/140‑1 (MQ=255) gCATTAAAGATGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGaa < 1:533615/140‑1 (MQ=255) aTGCGCGCGCCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATAtc > 1:755182/1‑140 (MQ=255) tGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATAtcg > 1:372132/1‑140 (MQ=255) | GTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCG > NC_000913/908656‑908926 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |