| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 4,169,232 | 0 | A | T | 23.1% | 136.5 / 11.2 | 53 | intergenic (+48/‑165) | rpoC/htrC | DNA‑directed RNA polymerase subunit beta'/heat shock protein |
| Reads supporting (aligned to +/- strand): ref base A (10/30); new base T (12/0); total (23/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.13e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGAT > REL606/4169161‑4169307 | cTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATggg < 2:1019290/76‑1 (MQ=255)cTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTggg > 2:1236895/1‑76 (MQ=255)cTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTGTggg > 2:436989/1‑76 (MQ=255) ggCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATgggggg < 1:775579/76‑1 (MQ=255) gCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGa < 2:1125070/76‑1 (MQ=255) cGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTgtgggggg > 1:1734755/1‑74 (MQ=38) ggTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTTTGGGGGGGGt > 2:58936/1‑76 (MQ=255) gTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTGGGGGGGGtt > 2:571431/1‑76 (MQ=38) cTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAg < 2:556485/76‑1 (MQ=255) tGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTTGGGGGGGTTTTgg > 2:998253/1‑76 (MQ=37) tGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTGGGGGGGGTTTTgg > 1:183826/1‑76 (MQ=38) tGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTGGGGGGGGTTTAgg > 2:1523208/1‑76 (MQ=38) gATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAggg < 2:1782064/76‑1 (MQ=255) gATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTTGGGGGAGTTTTggg > 2:1363233/1‑76 (MQ=38) gATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTATTGGGGGGGGTTAggg > 2:873962/1‑76 (MQ=38) aTAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGGGTTTAGGGa > 2:1616532/1‑76 (MQ=255) aTAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGa < 2:534325/76‑1 (MQ=255) tAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGaa < 2:837918/76‑1 (MQ=255) aaCGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTGGGGGGGGGTTAGGGgaa > 2:987199/1‑76 (MQ=37) gAGTAATCGTTAATCCGCAAATAACGTAAAACCCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAagag < 1:1493386/76‑1 (MQ=255) gAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGGGTTTAGGGAAagag > 2:45196/1‑76 (MQ=255) aGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGGTTTTTTTTTGGGGGGGGTTTGGGGaaaagg > 1:1617891/1‑73 (MQ=25) aGTAATCGTTAATCCGAAAATACCGTAAAACCCCGCTTCGGG‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGc < 1:586117/76‑1 (MQ=255) tAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCAt < 1:1809111/76‑1 (MQ=255) tAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTGGGGGGGGGTTTGGGGAAGAGCAt > 2:1305539/1‑76 (MQ=25) tAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTATGGGGGGAGGTTAGGGAAagaggaa > 2:1305532/1‑73 (MQ=38) aaTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTTGGGGGGGTTTTGGGAAAGAGCAtt > 1:1858610/1‑76 (MQ=37) cGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTATGGGGGGGGTTTAGGGAAAGAGCATTTGt > 2:336554/1‑76 (MQ=255) gTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTc < 2:784718/76‑1 (MQ=255) ttAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCa < 1:394619/76‑1 (MQ=255) tAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGGTTTTTTTTTGGGGGGGGTTTGGGGAAA‑AGGAATTGTCAg > 1:381548/1‑76 (MQ=21) tAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTACGGGGGGAGTTTAGGGAAAGAGCATTTGTCAg > 2:1564212/1‑76 (MQ=255) tAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTGGGGGGGGGTTAGGGAAAGAGCATTTGTCAg > 2:1150339/1‑76 (MQ=38) aaTCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTTGGGGGGGTTTAGGGGAAGAGCATTTGTCAGa > 2:1445722/1‑76 (MQ=37) ccGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAata > 2:1540505/1‑76 (MQ=255) cGCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTGGGGGGGGGTTTGGGGAAGAGGATTTGTCAGAatat > 1:1733389/1‑76 (MQ=25) gCAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTGGGGGGGGGTTAGGGAAAGAGCATTTGTCAGAATAtt > 2:1499437/1‑76 (MQ=38) cAAATAACGTAAAAGCCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATAttt < 2:1728930/76‑1 (MQ=255) cAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTTTTGGGGGGGGTTTGGGAAAGAGCATTTGTCAGAATAttt > 2:536481/1‑76 (MQ=25) cAAATAACGTAAAAACCCGCTTCGGC‑GGGGTTTTTTTTGGGGGGGGGTTTGGGAAAGAGCATTTGTCAGAATAttt > 1:1751178/1‑76 (MQ=37) aaCGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTaanta < 2:1314325/76‑4 (MQ=255) aaCGAAAAACCCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGa < 1:1454277/76‑1 (MQ=255) cGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAAt < 2:179845/76‑1 (MQ=255) gTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAAtt < 2:757699/76‑1 (MQ=255) tAAAAACCCGCTTCGGC‑GGGTTTTTTTTTGGGGGGGGTTTTGGGGAAGAGCATTTGTCAGAATATTTAAGGAAttt > 2:1789789/1‑76 (MQ=37) tAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAAttt < 2:1094135/76‑1 (MQ=255) aaaaaCCCGCTTCGGC‑GGGTTTTTTTATGGGGGGGGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAAtttt > 2:1068636/1‑75 (MQ=255) aaaaaCCCGCTTCGGC‑GGGGTTTTTTTTGGGGGGGGTTTTTGGAAAAAGCAGTTGTTTGAATATTTTAGGGAtttt > 2:1101116/1‑75 (MQ=16) aaaaaCCCGCTTCGGC‑GGGGTTTTTTTTGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTc > 2:105797/1‑76 (MQ=255) aaccccGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTg < 1:1267768/73‑1 (MQ=255) aaaCCCGCTTCGGCGGGGTTTTTTTTTGGGGGGGGTTT‑GGGGAAGAGCATTTGTCCGAAAATTTAAGGAATTTCTg > 2:741987/1‑76 (MQ=25) aaaCCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTg < 1:1815712/76‑1 (MQ=255) aaCCCGCTTCGGC‑GGGTTTTTTTTTGGGGGGGGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGa > 2:31634/1‑76 (MQ=255) cccGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAAt < 1:1457285/76‑1 (MQ=255) ccGCTTCGGCGGGGTTTTTTTTTGGGGGGGGTTTGGGGAAA‑AGGATTTTTTAGAATATTTTAAGAAATTTTTAATa > 1:1242389/1‑76 (MQ=9) tCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCAt < 1:1793868/76‑1 (MQ=255) tCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCAt < 1:1037909/76‑1 (MQ=255) tCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCAt < 1:18443/76‑1 (MQ=255) tCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCAt < 2:1380052/76‑1 (MQ=255) tCGGC‑GGGGTTTTTTTTTGGGGGGGTTTTGGGGAAAAGCATTTGTCAGAATATTTAAAGAAATTCTGAATAATCAt > 2:1055527/1‑76 (MQ=255) c‑ggGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATc > 2:875287/1‑76 (MQ=255) gggTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCa > 1:1837202/1‑76 (MQ=255) gTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAAt < 1:1460392/76‑1 (MQ=255) tttttttATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATg > 2:791175/1‑76 (MQ=255) ttttATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTag < 1:1236895/76‑1 (MQ=255) ttttATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTag < 1:1095993/76‑1 (MQ=255) ttATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTagag < 1:1876397/76‑1 (MQ=255) ttATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTagag < 1:616015/76‑1 (MQ=255) tATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAgaga > 1:1219024/1‑76 (MQ=255) aTGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGAt > 2:1339015/1‑76 (MQ=255) | CTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGAT > REL606/4169161‑4169307 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |