Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,086,499 | 0 | C | A | 41.3% | 65.3 / 22.9 | 46 | V100V (GTG→GTT) | rhaM | L‑rhamnose mutarotase |
Reads supporting (aligned to +/- strand): ref base C (6/21); new base A (19/0); total (25/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.07e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAA > NZ_CP009273/4086410‑4086583 | agaTGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAac < 2:1467089/90‑1 (MQ=255) atatGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGGAAAAAAattctt > 2:589435/1‑90 (MQ=255) gTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAg < 2:607868/90‑1 (MQ=255) gTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAAAACCTCTTTCAg > 1:910000/1‑90 (MQ=255) tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGGAAAAAACTTCTTTCAGc > 2:1252138/1‑90 (MQ=255) tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGGAAAAAAATTCTTGCCGc > 1:1159093/1‑90 (MQ=255) aTTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTccc < 2:1292318/90‑3 (MQ=255) aaaTTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCgggat < 2:1315400/90‑3 (MQ=255) aTTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGt < 2:619879/90‑1 (MQ=255) tAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAacac < 1:1281273/58‑1 (MQ=255) tAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTAt < 1:675882/90‑1 (MQ=255) tAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGGAAAaaac > 2:1281273/1‑58 (MQ=255) tAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGGAAAAAACTTCTTTTAGCTCCCTACT‑CACCCGGGTGtttt > 2:675882/1‑88 (MQ=255) aaaTATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATcc < 2:512162/90‑1 (MQ=255) atTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGt < 2:1044671/90‑1 (MQ=255) tgtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGGAAAAAACCTTTTTCTGCTCGCGCCT‑CTCCCGGGGGTTATCCGGGGTa > 1:806021/1‑90 (MQ=255) tgtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGGAAAAAACCTTTTTCAGCTCGCCACT‑CACCGGGCTGTTTTTCGGGGTa > 1:57378/1‑90 (MQ=255) gttccGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAg < 1:1111399/87‑1 (MQ=255) gtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGTAAAACACTTCTTGCAGCTCGCTACT‑CAc > 2:851313/1‑71 (MQ=255) gtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGTAAAACACTTCTTGCAGCTCGCTACT‑CAc < 1:851313/71‑1 (MQ=255) gtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGTAAAACACTTCTTGCAGCTCGCTACT‑CAc < 1:1100286/71‑1 (MQ=255) gtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGTAAAACACTTCTTGCAGCTCGCTACT‑CAc > 2:1100286/1‑71 (MQ=255) gtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGGAAAAAACTTCTTTCAGCCCGCTACT‑CACCCGGCTGTTATCCggggttg > 2:1329998/1‑86 (MQ=255) tgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGc < 2:1042269/90‑1 (MQ=255) gCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCg < 2:1319939/90‑1 (MQ=255) ggAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGGAAAAAACATCTTGCCGCTCGCTTTT‑TCCCGGGCAGTTTTTCGGGGTTGCGGGc > 1:1042269/1‑90 (MQ=255) ggAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGGAAAAAAACTCTTTCAGCCCGCCACT‑CCCCCGGCTGTTTTTCGGGGTTGCgggg > 1:698408/1‑89 (MQ=255) aGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGGGAAAAACTTCTTGCAGCTCGCTACT‑CCCCCGGCTGTTATCCGGGTTAGCGGGcct > 1:825478/1‑88 (MQ=255) aGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGGAAAAAACTTCTTTCAGCCCGCTACT‑CACCGGGCGGTTATTCGGGTTATCGGGcct > 1:1152837/1‑88 (MQ=255) cATCTCTGGCACATTGGGCAATTACGGCA‑GGGAAAACACTTCTTTCCGCTCGCTGCT‑CCCCCGGCTGTTTTTCGGGGTAACGGGGATaaa > 1:1433078/1‑89 (MQ=255) tctcTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACAt < 2:101633/90‑1 (MQ=255) tctcTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACAt < 2:1433078/90‑1 (MQ=255) tGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGt < 2:271880/90‑1 (MQ=255) tGGCACATTGGGCAATTACGGCA‑GGGAAAAAACTTCTTTCAGCCCGCTCCTCCCCCGGGCGGTTTTCCGGGTT‑GCGGGGATAAAATCggg > 1:739966/1‑89 (MQ=255) ggCACATTGGGCAATTACGGCA‑GGGAAAAAACTTCTTTCAGCTCGCTACT‑CACCCGGCTGTTTTCCGGGTTAGCGGGCATAACATCGGTc > 2:799255/1‑90 (MQ=255) aTTGGGCAATTACGGCA‑GGGGAAAAAATTCTTTCAGCGCGCTACT‑CCCCCGGGGGTTTTTCGGGGTAGGGGGCATTAAATCTGTCatata > 1:1186607/1‑90 (MQ=255) ttGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCt > 2:1174400/1‑53 (MQ=255) ttGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCt < 1:1174400/53‑1 (MQ=255) ttGGGCAATTACGGCA‑GGTAAAAAACTTCTTGCAGCTCCCTACT‑CACCCGGCTGTTATCCCGGGTAGGGGGCATAACATCGGGCatatat > 2:23069/1‑90 (MQ=255) tGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATAtt < 2:1159093/90‑1 (MQ=255) gggCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCg > 2:581646/1‑67 (MQ=255) gggCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCg < 1:581646/67‑1 (MQ=255) gggCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATAttt < 1:1329998/90‑1 (MQ=255) gggCAATTACGGCA‑GGGAAAACACTTCTTTCAGCTCGCTACT‑CACCCGGCTGGTATCCGGGTTAGCGGGCATAAAAACGGGCATATAttt > 2:470531/1‑90 (MQ=255) gggCAATTACGGCA‑GGGAAAAAACTTCTTTCAGCCCCCCACT‑CACCCGGGTGTTTTTCCGGGTAGGGGGGATAAAAATGGTCATATAttt > 1:242054/1‑90 (MQ=255) gCAATTACGGCA‑GGTAAAAAACTTCTTTCAGCTCCCTACT‑CACCCGGCTGTTTTCCGGGGTTGCGGGCATAAAAACGGTCATATAttttc > 2:928984/1‑88 (MQ=255) gCAATTACGGCA‑GGGAAAAAACTTTTTTCAGCTCGCTACT‑CACCCGGCTTTTTTCCGGGGTAGCGGGGATAAAATCGGTCATATAttttc > 1:1073678/1‑88 (MQ=255) gCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGtt < 2:242054/90‑1 (MQ=255) gggaaaaaaCTTCTTTCAGCTCGCTACT‑CCCCCGGCTGTTATCCGGGGTAGAGGGCAAAACATCGGTCCTATATTTCCCCCCCCGTTTGc > 1:1233416/4‑90 (MQ=255) tAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCaa < 1:1348123/90‑1 (MQ=255) | AGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCA‑GGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAA > NZ_CP009273/4086410‑4086583 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |