Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,550,663 | 0 | A | C | 28.9% | 88.0 / 15.5 | 45 | N68T (AAC→ACC) | fdnH | formate dehydrogenase N subunit beta |
Reads supporting (aligned to +/- strand): ref base A (30/2); new base C (0/13); total (30/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.44e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTT‑CCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTC > NC_000913/1550530‑1550801 | cGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTT‑CCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgc > 2:694466/1‑140 (MQ=255) tGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTT‑CCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgccg > 1:974344/1‑140 (MQ=255) aTTACAAAGCAGAAGTCGCAAAACTTATCGACGTTT‑CCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGAt < 1:213822/140‑1 (MQ=255) aCAAAGCAGAAGTCGCAAAACTTATCGACGTTT‑CCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTg > 1:216583/1‑140 (MQ=255) aGCAGAAGTCGCAAAACTTATCGACGTTT‑CCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAgcg > 1:820629/1‑140 (MQ=255) aaGAAGTCCCAAAGTTTATCG‑CGTTTCCCACGTGGTTGGGCTGAAAAGCTTGGCGGGTGGCGTGTTGGGAGGGAAACGACATCCGGGATGAGGGGGGGACCGGCGGCGGGGTTTAGGATACCCCCGCCGATCTGAGCGcc < 2:93550/139‑1 (MQ=255) ttATCGCCTTTT‑CCCCCTGTATCGGCTGTAAAGCCTGTAAGGGGGCGTGTGGGGGGGGGAACGACATCCGGGATGAAGTGGGGCACTGCGGGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTg < 2:150231/140‑1 (MQ=255) ttATCGACTTTT‑CCACCTGTATCGGCTGTAAAGCCTGTCGGGTGGCGTGTTGGGAGTGGAACGACATCGGTGATGAGGGGGGGCACTGCGTGGGGTTTTACAATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTg < 2:943939/140‑1 (MQ=25) aTCGACGTTT‑CCCCCTGTATCGGCTGTAAAGCCTGTCGGGGGGCGTGTTCGGAGTGGAACGACATCGGTGATGAGGTGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGAt < 2:452150/140‑1 (MQ=37) ttt‑CCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCttt > 1:1224216/1‑140 (MQ=255) tt‑CCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGTGCTTTa > 1:335428/1‑140 (MQ=255) tATGGGCTGTAAAGCCTGTCGGGTGGCGTGTTGGGAGTGGAACGACATCCGTGATGAGGGGGGGCACTGGGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGTGCTTTAGCGAAACCg < 2:335428/140‑1 (MQ=255) aTCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGa > 1:372940/1‑140 (MQ=255) aTCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGa > 1:707203/1‑140 (MQ=255) tCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGaa > 2:57803/1‑140 (MQ=255) tGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCt > 2:1208093/1‑140 (MQ=255) tGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCt > 1:993815/1‑140 (MQ=255) cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGa > 1:960017/1‑140 (MQ=255) cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGa > 1:1184611/1‑140 (MQ=255) cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGa > 2:881767/1‑140 (MQ=255) ggCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGc > 1:16865/1‑140 (MQ=255) gCGTGTTGGGAGTGGAACGACATCCGTGATGAAGTGGGGCCCTGGGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCt < 2:1184611/140‑1 (MQ=255) ttCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc > 1:975723/1‑140 (MQ=255) ttCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc > 1:1388628/1‑140 (MQ=255) acaggTGGAACGACATCGGTGATGAGGGGGGGCCCGGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGt < 2:1028218/136‑1 (MQ=255) aGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTaaa > 2:770044/1‑140 (MQ=255) gggAACGACACCCGTGATGAAGGGGGCCACTGCGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGa < 2:372940/139‑1 (MQ=255) ggAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGAc > 1:692124/1‑140 (MQ=255) gAACGACATCCGTGATGAGGTGGGGCACTGGGGCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg < 2:974344/140‑1 (MQ=255) gAACGACATCCGTGATGAAGGGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg < 2:820629/140‑1 (MQ=255) cATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTa > 1:960051/1‑140 (MQ=255) cATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTa > 2:1259423/1‑140 (MQ=255) tCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg < 1:770044/140‑1 (MQ=255) tCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTgctgt > 1:1028218/1‑136 (MQ=255) ccggTGATGAAGTGGGCCACGGCGTGGGGTTTTACAATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg < 2:975723/137‑1 (MQ=255) cGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg < 1:928030/130‑1 (MQ=255) cGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg > 2:928030/1‑130 (MQ=255) gTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCAc > 1:402666/1‑140 (MQ=255) gTGATGAAGGGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCAc < 2:16865/140‑1 (MQ=255) gatgaAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACtg > 2:600643/1‑140 (MQ=255) cACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg < 1:57803/140‑1 (MQ=255) cACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg < 1:881767/140‑1 (MQ=255) cACTGCGGGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg < 2:216583/140‑1 (MQ=255) cTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCt > 2:138217/1‑140 (MQ=255) gTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCt < 1:694466/140‑1 (MQ=255) gTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGc > 2:1145334/1‑140 (MQ=255) tttACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCg > 2:1377694/1‑140 (MQ=255) tttACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCg > 1:204799/1‑140 (MQ=255) tACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTg > 1:1275573/1‑140 (MQ=255) aTAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCg < 1:600643/140‑1 (MQ=255) aaCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTc > 2:1353753/1‑140 (MQ=255) aaCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTc > 2:1003027/1‑140 (MQ=255) | CGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTT‑CCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTC > NC_000913/1550530‑1550801 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |