| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,325,888 | 0 | G | T | 25.0% | 56.6 / 14.8 | 28 | intergenic (‑147/+511) | yjdN/yjdM | conserved protein YjdN/conserved protein YjdM |
| Reads supporting (aligned to +/- strand): ref base G (9/12); new base T (1/6); total (10/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.64e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-02 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCG > NC_000913/4325749‑4326020 | aCCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAg > 2:388897/1‑140 (MQ=255) gcggcgtgaacgccttatccggtctacacttcgcccataaaccgtaggcctgataagacgcgttagcgtcgcatcaggcatctgcgcacgacTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa < 1:341908‑M2/48‑1 (MQ=255) gcggcgtgaacgccttatccggtctacacttcgcccataaaccgtaggcctgataagacgcgttagcgtcgcatcaggcatctgcgcacgacTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa < 2:138950‑M2/48‑1 (MQ=255) aaGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa < 2:223266/140‑1 (MQ=255) cgtgaacgccttatccggtctacacttcgcccataaaccgtaggcctgataagacgcgttagcgtcgcatcaggcatctgcgcacgacTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAAttt > 2:420136‑M2/89‑140 (MQ=255) gtgaacgccttatccggtctacacttcgcccataaaccgtaggcctgataagacgcgttagcgtcgcatcaggcatctgcgcacgacTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTg < 1:254499‑M2/53‑1 (MQ=255) gtgaacgccttatccggtctacacttcgcccataaaccgtaggcctgataagacgcgttagcgtcgcatcaggcatctgcgcacgacTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTg < 1:304432‑M2/53‑1 (MQ=255) gtAGTTCGAAATTTATACAGATGAGAGG‑GGCCTGATAAGACGCGCGAGGGTCGCATCAGGCAGACGGTCCTGTTGCCGGATGCGGAGTAAACGCCTTATCCGACCTAAGTAGGGTGCGGGAATTTGTAGGACTGATAAGa < 2:489270/140‑1 (MQ=17) ccggtctacacttcgcccataaaccgtaggcctgataagacgcgttagcgtcgcatcaggcatctgcgcacgacTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAg < 1:420136‑M2/66‑1 (MQ=255) cGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACgcgc > 2:261538/1‑140 (MQ=255) cGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACgcgc > 1:450364/1‑140 (MQ=255) tttATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg < 2:43215/140‑1 (MQ=255) tttATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg < 2:201902/140‑1 (MQ=255) ataCAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCg > 1:472816/1‑140 (MQ=255) aTGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg < 2:238512/140‑1 (MQ=255) agagGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg > 2:457011/1‑128 (MQ=40) agagGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg < 1:457011/128‑1 (MQ=40) agGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGt < 1:388897/140‑1 (MQ=38) tGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCa < 1:454856/140‑1 (MQ=35) cgcgcgAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGa > 1:111197/1‑140 (MQ=37) cgAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcg > 2:468434/1‑140 (MQ=37) aTCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGc < 1:261538/140‑1 (MQ=35) tCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCgg < 2:111197/140‑1 (MQ=33) gCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt < 1:184305/140‑1 (MQ=32) tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAgg > 2:209276/1‑140 (MQ=32) cGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTg < 1:155482/140‑1 (MQ=17) gcggcgTGAACGCCTTATCCGTCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa < 2:126760/140‑1 (MQ=17) gCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg > 1:50902/1‑124 (MQ=32) gCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg < 2:50902/124‑1 (MQ=32) cTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgcg > 2:155482/1‑140 (MQ=32) | ACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCG > NC_000913/4325749‑4326020 |
| 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: ‑ |