| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,596,055 | 0 | T | G | 66.7% | 29.0 / 18.2 | 24 | pseudogene (1933/3861 nt) | yneO | pseudogene, AidA homolog |
| Reads supporting (aligned to +/- strand): ref base T (4/4); new base G (8/8); total (12/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCGTCATTGACCAGCGTTTTACCGCTATAAGTATTCCCTGCATTGAACACCGTATGGCCTGCATCATGGATCACTTTTCCGTCTTTATCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGTATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACAC > NC_000913/1595916‑1596189 | ccGTCATTGACCAGCGTTTTACCGCTATAAGTATTCCCTGCATTGAACACCGTATGGCCTGCATCATGGATCACTTTTCCGTCTTTATCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGGATGATTaaa < 2:47028/149‑1 (MQ=255) aCCAGCGTTTTACCGCTATAAGTATTCCCTGCATTGAACACCGTATGGCCTGCATCATGGATCACTTTTCCGTCTTTATCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGGATGATTAAAGACAAAAAc < 3:113583/149‑1 (MQ=255) aGTATTCCCTGCATTGAACACCGTATGGCCTGCATCATGGATCACTTTTCCGTCTTTATCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGGATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACt > 1:161787/1‑149 (MQ=255) ccTGCATTGACAACAGTATGGCCTGCATCATGTATCACTTATACGTCTTTATAGTCAACTGTTATGAGAATATAGACCTGGTAGCCGGCATCACTGTGATTCGGATGATTAAAGCCAAAAACGCATTCGCCCAGACCAAACTCCACtt > 2:122098/1‑148 (MQ=255) cATTGAACACCGTATGGCCTGCATCATGGATCACTTTTCCGTCTTTATCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGGATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTtcgtc < 2:161787/149‑1 (MQ=255) cacCGTATGGCCTGCATCATGGATCACTTTTCCGTCTTTATCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGTATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTgg > 6:156729/1‑149 (MQ=255) cGTATGGCCTGCATCATGGATCACTTTTCCGTCTTTATCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGTATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTGGTGa < 5:156729/149‑1 (MQ=255) tGCATCATGGATCACTTTTCCGTCTTTATCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGGATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGTATTGGTGATAAATCCgg > 7:8120/1‑149 (MQ=255) atGGATCACTTTTCCGTCTTTATCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGTATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCt < 4:121791/148‑1 (MQ=255) cTTTTCCGTCTTTATCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTAATGGGATGATTAAAGACAAAAAAGCCTGCGCCAATACCAAACTCCACTTTCTTCGCATGTGTGATAAATACGGCATCTgccgcagc > 4:20212/1‑148 (MQ=255) cGTCTTTATCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGTATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTgg > 2:17502/1‑117 (MQ=255) tCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGgatgat < 8:134206/58‑1 (MQ=255) tCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGgatgat > 7:134206/1‑58 (MQ=255) gCATATCGACCTGGTAGCCGGCATCACTGTTATTGGGATGATTAAAGa < 3:103863/48‑1 (MQ=255) gCATATCGACCTGGTAGCCGGCATCACTGTTATTGGGATGATTAAAGa > 4:103863/1‑48 (MQ=255) tatCGACCTGGTAGCCGGCATCACTGTTATTGGTATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGagca < 8:44234/128‑1 (MQ=255) tatCGACCTGGTAGCCGGCATCACTGTTATTGGTATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGAGagca > 7:44234/1‑128 (MQ=255) aCCTGGTAGCCGGCATCACTGTTATTGGTATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTAcag > 8:58782/1‑148 (MQ=255) aCCTGGTAGCCGGCATCACTGTTATTGGGATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTAcagc > 6:78073/1‑149 (MQ=255) gCCGGCATCACTGTTATTGGGATGATTAAAGACAAAAACGCCTTCGCCAAGACCACACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTAca < 1:136958/139‑1 (MQ=255) gCCGGCATCACTGTTATTGGGATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTAca > 2:136958/1‑139 (MQ=255) ggCATCACTGTTATTGGTATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCaaaaa < 7:58782/149‑1 (MQ=255) ggaTCACTGTTATTGGGATGATTAAAGACAAAAACGCCTTCTCCAAGACCAAACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAac < 1:122098/147‑1 (MQ=255) aTCACTGTTATTGGGATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGTATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAacac < 8:8120/149‑1 (MQ=255) ttGGGATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGagca < 2:74484/99‑1 (MQ=255) ttGGGATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGagca > 1:74484/1‑99 (MQ=255) | CCGTCATTGACCAGCGTTTTACCGCTATAAGTATTCCCTGCATTGAACACCGTATGGCCTGCATCATGGATCACTTTTCCGTCTTTATCGTCACCTGTAATCAGCATATCGACCTGGTAGCCGGCATCACTGTTATTGGTATGATTAAAGACAAAAACGCCTTCGCCAAGACCAAACTCCACTTTCGTCGCATTGGTGATAAATCCGGCATCTGCCGCCGCCTCACCGTGAGCAGCACCAATGTTTAGCGTTCCTACAGCAGGCTCAAAAACAC > NC_000913/1595916‑1596189 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |