Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 760,181 | 0 | G | C | 54.5% | 4.3 / 13.7 | 11 | P182P (CCG→CCC) | sucD | succinate‑‑CoA ligase subunit alpha |
Reads supporting (aligned to +/- strand): ref base G (0/5); new base C (6/0); total (6/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.16e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGC‑GGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGG > NZ_CP009273/760098‑760270 | ggTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGC‑GGTGACCCGATCCCg < 1:222298/90‑1 (MQ=255) aTGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCCACCCCGGGCCCTAAtttt > 1:485432/1‑90 (MQ=255) aTGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGC‑GGGGAACCCCACCCCGGCGCTAACtttt > 1:278338/1‑89 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGC‑GGTGACCCCATCCCCGGCTCTTACTTTTTTGACAt > 2:119912/1‑90 (MQ=255) aGACCACGGATTACGGTTTCGGTCAGTCGACCTGGGTCGGTATCGGC‑GGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGaaa < 2:451994/90‑1 (MQ=255) gACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGC‑GGTGACCCCATCCCCGGGTCTTACTTTTTTGGCATTTTCCAAAt > 2:118811/1‑90 (MQ=255) ccACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGC‑GGTGACCCCATCCCCGGCTCTAACTTTTTCGACATTCTCGAAATGt > 2:102636/1‑90 (MQ=255) cAGTCGACCTGTGTCGGTATCGGC‑GGTGACCCCATCCCCGGCCCTAAATTTTTCCACATTTTGGAAAAGTTTGAAAAAAAACCGCAGAcc > 2:409466/1‑90 (MQ=255) gACCTGTGTCGGTATCGGC‑GGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGc < 1:102636/90‑1 (MQ=255) ccTGTGTCGGTATCGGC‑GGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGa < 2:150393/90‑1 (MQ=255) ggc‑ggTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGt < 2:172780/90‑1 (MQ=255) gATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATcggcgg > 1:145881/1‑90 (MQ=255) gATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATcggcgg > 1:231927/1‑90 (MQ=255) | GGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGC‑GGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGG > NZ_CP009273/760098‑760270 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |