Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 760,187 | 0 | G | C | 48.8% | 38.6 / 24.1 | 41 | P184P (CCG→CCC) | sucD | succinate‑‑CoA ligase subunit alpha |
Reads supporting (aligned to +/- strand): ref base G (6/15); new base C (20/0); total (26/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.10e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGC > NZ_CP009273/760098‑760276 | ggTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCg < 1:488974/90‑1 (MQ=255) cTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCATCCCCGGgtct > 2:689264/1‑90 (MQ=255) cTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCATCCCCGGctct > 2:533265/1‑90 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGGCCCCGTCCCCGGgtctta > 2:826848/1‑88 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCGATCCCCGGCTCTaa > 2:723804/1‑90 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCGTCCCCGGctctta > 2:2086857/1‑88 (MQ=255) aCCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAAc < 1:357045/90‑1 (MQ=255) gCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGAc < 1:2003103/90‑1 (MQ=255) gCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGAc < 2:1870867/90‑1 (MQ=255) gCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCCATCCCCGGCTCTTACTTTTTTGAc > 2:428799/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGGCCCCCTCCCCGGGCCTTACTTTTTTGACAt > 2:1885681/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCCATCCCCGGCTCTAACTTTTTTGACAt > 2:491929/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCCATCCCCGGCTCTAACTTTTTTGACAt > 2:512336/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGGCCGCTCCCCGGCGCTAACTTTTTCGAcct > 1:923774/1‑88 (MQ=255) aGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGaaa < 1:826848/90‑1 (MQ=255) ccACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGt > 2:486667/1‑90 (MQ=255) ccACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCCATCCCCGGCTCTAACTTTTTCGACATTCTCCAAATGt > 2:1450979/1‑90 (MQ=255) ccACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCAACCCCGGGCCTAACTTTTTTTACATTCTCCAAAAGt > 1:348442/1‑90 (MQ=255) ccACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGAACCGATCCCCGGCGCTAACTTTTTTTGCATTTTCCAAATgg > 1:815156/1‑89 (MQ=255) ccACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGAACCCCTCCCCGGGGCTAACTTTTTTGGAAATCTCGaaaaag > 2:1747775/1‑87 (MQ=255) aCGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTc < 1:1196417/90‑1 (MQ=255) aCGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCATCCCCGGCCCTAACTTTTTCGACATTCTCTAAAAGttt > 2:1222426/1‑89 (MQ=255) gATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGaa < 2:1170199/90‑1 (MQ=255) gATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCATCCCGGGCTCTAACTTTTTCGACATTCTCGGAATGTTTGaa > 2:695434/1‑90 (MQ=255) gTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCATCCCCGGCTCCTACTTTTTTTACAATTTCGGAAAGTTTGaaaaaaaaa > 1:1581930/1‑86 (MQ=255) gTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCATCCCCGGCTCTAACTTTTTTTACATTCTCGAAATGTTTGaaaaaaaaa > 1:375653/1‑86 (MQ=255) ggTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAg < 2:1783386/90‑1 (MQ=255) cAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGAcc < 2:671276/90‑1 (MQ=255) cAGTCGACCTGTGTCGGTATCGGCGGTGACCCCCTCCCCGGGTCTAACTTTTTTGACATTTTCGAAATGTTTTAAAAAAAACCGCAGAcc > 1:677919/1‑90 (MQ=255) gTCGACCTGTGTCGGTATCGGCGGTGACCCCATCCCCGGCTCTAACTTTTTTGACATTTTCGGAATGTTTGGAAAAAAACCGCAGACCGa > 1:455376/1‑90 (MQ=255) gTCGACCTGTGTCGGTATCGGCGGTGACCCCATCCCCGGCTCTAACTTTTTCGACATTTTCGAAATGTTTGAAAAAGAACCGCCGAccca > 2:160217/1‑88 (MQ=255) cTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGAt < 1:1822274/90‑1 (MQ=255) gtgtCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCg < 1:869000/90‑1 (MQ=255) ggTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTTTTGACATTCTCGAAAAGTTTGAAAAAAAACCCCAGGCCGAAGCGAACgtggtg > 2:258726/1‑90 (MQ=255) gTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGtgatga < 2:429140/90‑1 (MQ=255) cggcggTGACCCGATCCCGGGCTCTAACTTTTTCGACATTTTCGAAATGTTCGAAAAAAATCCGCAGACCGAAGCGATCGTGATGATCgg > 2:1304700/1‑90 (MQ=255) cggcggTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCgg < 2:1625878/90‑1 (MQ=255) cggTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTga < 2:992760/90‑1 (MQ=255) aCCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATcg < 1:1885681/90‑1 (MQ=255) gATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAg > 1:2000516/1‑55 (MQ=255) gATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAg < 2:2000516/55‑1 (MQ=255) aTCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGt > 1:2062475/1‑90 (MQ=255) gggCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCg < 2:1399497/58‑1 (MQ=255) gggCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCg > 1:1399497/1‑58 (MQ=255) gggCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAgcgc < 2:2062475/90‑1 (MQ=255) gggCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAgcgc < 1:732022/90‑1 (MQ=255) gggCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAgcgc > 1:1944487/1‑90 (MQ=255) | GGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGC > NZ_CP009273/760098‑760276 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |