| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 2,231,892 | 0 | A | C | 37.8% | 48.1 / 17.0 | 37 | T25P (ACC→CCC) | yejH | predicted ATP‑dependet helicase |
| Reads supporting (aligned to +/- strand): ref base A (20/3); new base C (0/14); total (20/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.11e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCATGATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTCGT > REL606/2231818‑2231966 | ccATGATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAc > 2:214827/1‑76 (MQ=255) aTGATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAccc < 1:1265687/76‑1 (MQ=255) gATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCt > 1:1304463/1‑76 (MQ=255) tttACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCg > 1:1804668/1‑76 (MQ=255) tACACTTCCCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGtt < 2:1907488/76‑1 (MQ=255) cacTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTAt < 1:1937579/76‑1 (MQ=255) acTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATc > 1:105942/1‑76 (MQ=255) ctccGCCCATATCAGAAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCg < 1:821901/73‑1 (MQ=255) cccATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGtgctg > 2:88830/1‑76 (MQ=255) ccATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGtgctgc > 1:38874/1‑76 (MQ=255) cATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGcc > 2:90364/1‑76 (MQ=255) cATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGcc > 1:2034639/1‑76 (MQ=255) atatCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGccg > 2:112328/1‑75 (MQ=255) atCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCAc < 2:270989/76‑1 (MQ=255) aGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCgg > 2:646028/1‑76 (MQ=255) aGCAAGAACCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCgg < 1:646749/76‑1 (MQ=255) aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGgcgc < 2:1729966/76‑1 (MQ=255) aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGgcgc > 1:49963/1‑76 (MQ=255) aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGgcgc < 2:644256/76‑1 (MQ=255) agaagCGGGGGATGCCCCGCTCATCCTTTTTCTTCGTCACAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCa < 2:1860839/76‑1 (MQ=255) aagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAgg > 2:1864921/1‑76 (MQ=255) gCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTa < 1:1284541/76‑1 (MQ=255) cGTGGATGCCCCGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTaaa < 1:2047950/76‑1 (MQ=255) cGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTaaa < 1:951498/76‑1 (MQ=255) gTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTaaaa < 2:41646/76‑1 (MQ=255) gTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTaaaa > 2:668793/1‑76 (MQ=255) tGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg < 2:1982627/76‑1 (MQ=255) aTGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGCAAAAGcat > 2:1070460/1‑74 (MQ=255) tGCCCCGCTCACCCATTTTCGTCCTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTg < 1:159626/76‑1 (MQ=255) tCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGc > 2:1208594/1‑76 (MQ=255) tCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGc > 1:340829/1‑76 (MQ=255) ccATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGaa > 1:342276/1‑76 (MQ=255) aTTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACt < 1:1384801/76‑1 (MQ=255) aTTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACt > 2:1092648/1‑76 (MQ=255) aTTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACt > 2:708575/1‑76 (MQ=255) aTTTTCGTCGTCAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAAcn < 1:846007/76‑2 (MQ=255) ttCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGc > 1:1767309/1‑76 (MQ=255) tcgtcAAAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACgg < 1:247079/76‑1 (MQ=255) gtcATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACggct < 2:1580557/76‑1 (MQ=255) tcATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACggctg < 2:832114/76‑1 (MQ=255) aaaaCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTc > 1:1571503/1‑76 (MQ=255) aaaaCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTc > 1:837901/1‑76 (MQ=255) aaaaCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTc > 1:350550/1‑76 (MQ=255) aaCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTCGt > 1:984602/1‑76 (MQ=255) | CCATGATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTCGT > REL606/2231818‑2231966 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |