| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_pae_tpiA | 4,262,166 | 0 | G | C | 22.7% | 28.7 / 12.2 | 22 | G194G (GGG→GGC) | dusA | tRNA‑dihydrouridine synthase A |
| Reads supporting (aligned to +/- strand): ref base G (12/5); new base C (3/2); total (15/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-02 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGA > NC_000913_3_pae_tpiA/4262018‑4262312 | tggtgTCGAGTCCGGTGACGGTGCAAACGCGTATGGTCCTCGACGACCAGGACTGCTATGAATTTCTCTGCAATTTCATCAATTCCGGTTCTGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCgggg < 2:78970/149‑1 (MQ=255) gtCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGCGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCTGGGTTag > 2:351468/1‑148 (MQ=255) 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) | TGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGA > NC_000913_3_pae_tpiA/4262018‑4262312 |
| 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: ‑ |