Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_pae_tpiA | 4,262,171 | 0 | G | C | 21.7% | 35.5 / 11.4 | 23 | S196T (AGC→ACC) | dusA | tRNA‑dihydrouridine synthase A |
Reads supporting (aligned to +/- strand): ref base G (11/7); new base C (3/2); total (14/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTG > NC_000913_3_pae_tpiA/4262023‑4262319 | atgATTCTGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg > 2:261171/3‑149 (MQ=255) agtgcggtgacggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGccc < 2:581347/145‑1 (MQ=255) ggtgaAAGCGCGTATTGGCATCGACGACCAGGACAGGTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGTTTTGTTTATCATCGACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCACGAAAGGAAAcc > 2:250597/1‑149 (MQ=255) ggtaagtGGTATCGACGGCAGGGACGGCGATGAAGTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATcc < 2:590003/143‑1 (MQ=255) gacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATcc > 1:540399/1‑149 (MQ=255) tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGACTATCCGCATGTCTATCa > 2:429520/1‑149 (MQ=255) tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGa > 2:507189/1‑149 (MQ=255) cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaa > 2:618344/1‑70 (MQ=39) cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaa < 1:618344/70‑1 (MQ=39) gCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGCCTTTCCGCATCTCACAATGTCGATTAACggt > 2:35548/1‑149 (MQ=255) gaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAACCGTGACTTTCCGCATCTGACAATGTCGATTCAAGTGGGTATCCAGt > 2:275439/1‑149 (MQ=255) ttCATCATCCACGCACGTAAAGCCTGGCTTTCGGTGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCt > 2:686447/1‑149 (MQ=255) ccACGCACGTAAAGCCGGGGTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTACCGGTGGTATCAAGTCGCTGGAAGAgg < 2:394821/149‑1 (MQ=255) cacgcaggTAAAGCCTGGGTTTCGGGGTTAAGCACGAAAGAAAACGGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCTACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGc < 2:616648/149‑1 (MQ=255) tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCACGTCGCTGGAAGAGGCCAAAGcac > 2:92344/1‑149 (MQ=255) aaaGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAGTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCCACTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGcaca < 2:42860/149‑1 (MQ=255) gCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCt > 2:533251/1‑149 (MQ=255) cTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGCGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGCCCAAAGCACACCTGc > 1:657343/1‑149 (MQ=255) ttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCaa > 2:386246/3‑149 (MQ=255) tcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAgg < 2:609343/130‑1 (MQ=255) tcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAgg > 1:609343/2‑131 (MQ=255) gggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGTAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatgga < 1:616075/149‑1 (MQ=255) gggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCt < 2:43234/83‑1 (MQ=255) gggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCt > 1:43234/1‑83 (MQ=255) gCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGATGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTg > 2:477891/1‑149 (MQ=255) | TCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTG > NC_000913_3_pae_tpiA/4262023‑4262319 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |