Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,169,232 | 0 | A | T | 28.1% | 138.2 / 14.6 | 57 | intergenic (+48/‑165) | rpoC/htrC | DNA‑directed RNA polymerase subunit beta'/heat shock protein |
Reads supporting (aligned to +/- strand): ref base A (11/30); new base T (16/0); total (27/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGAT > REL606/4169160‑4169307 | tCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGGTTTTTTTTgg > 1:1941509/1‑76 (MQ=255) cTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATggg > 2:562425/1‑76 (MQ=255) ggCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAACCCCGCTTCGGC‑‑GGGTTTTTTTATgggggg < 1:483738/76‑1 (MQ=255) ggCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATgggggg < 2:1016308/76‑1 (MQ=255) gTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCG‑GGGTTTTTTTTTggggggggg > 1:528456/1‑73 (MQ=37) gTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCG‑GGGTTTTTTTTTGGGGGGGGt > 1:1716423/1‑76 (MQ=37) ttCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGttt < 1:569577/76‑1 (MQ=255) tCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTa < 2:1179824/76‑1 (MQ=255) tCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGGTTTTTTTTGGGGGGGGtttt > 2:615862/1‑75 (MQ=38) cTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCG‑GGGTTTTTTTTTGGGGGGGGtttg > 1:1164457/1‑75 (MQ=37) cTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCG‑GGGTTTTTTTTTGGGGGGGGtttg > 2:791796/1‑75 (MQ=37) cTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCG‑GGGTTTTTTTTTGGGGGGGGtttg > 1:322801/1‑75 (MQ=37) tGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGGTTTTTTTTGGGGGGGGTTTAgg > 2:211064/1‑76 (MQ=38) gATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCG‑GGGTTTTTTTTTGGGGGGGGTTTggg > 1:587582/1‑76 (MQ=37) gATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTTTGGGGGGGGGTTAggg > 2:497691/1‑76 (MQ=38) gATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAggg < 2:410634/76‑1 (MQ=255) tAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGaa < 2:1941509/76‑1 (MQ=255) tAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGGTTTTTTTTGGGGGGGGTTTAGGGaa > 2:237478/1‑76 (MQ=38) tAACGAGTAACCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCG‑GGGTTTTTTTTTGGGGGGGGGTTggggg > 1:596026/1‑75 (MQ=25) aaCGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCCCCTCGGCGGGGGTTTTTTTTTGGGGGGGGTTT‑GGGaa > 1:1490834/1‑76 (MQ=255) aCGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAag < 2:1741231/76‑1 (MQ=255) aCGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAag < 1:830972/76‑1 (MQ=255) gAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAagag < 2:994278/76‑1 (MQ=255) aaTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTTTTGGGGGGGGTTAGGGGAAGAGGAtt > 2:552519/1‑76 (MQ=25) aaTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTTTGGGGGGGGTTTAGGGAAAGAGCAtt > 2:94117/1‑76 (MQ=255) aaTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCAtt < 1:153778/76‑1 (MQ=255) cGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGt < 1:1964988/76‑1 (MQ=255) gTTAATCCGCAAATAACGTAAAACCCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTc < 2:1345929/76‑1 (MQ=255) gTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGGTTTTTTTTGGGGGGGGTTTAGGGAAAGAGCATTTGTc > 1:1900204/1‑76 (MQ=38) ttAATCCGCAAATAACGTAAACACCCGCTTCGGC‑‑GGGGTTTTTTATTGGTGGGGTTTAGGGGAACGGCATTTGtta > 2:967544/1‑74 (MQ=25) ttAATCCGCAAATAACGTAAAAACCCGCTTCGGCG‑GGGTTTTTTTTTGGGGGGGGTTT‑GGGGAAGAGGAATTGTCa > 2:1316895/1‑76 (MQ=21) ttAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCa < 1:1763418/76‑1 (MQ=255) ttAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGGTTTTTTTTGGGGGGGGTTTAGGGAAAGAGCATTTGTCa > 2:1900767/1‑76 (MQ=38) ttAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGGTTTTTTTTGGGGGGGGGTTTGGGAAAGGGGAATTGTcg > 2:1192721/1‑75 (MQ=21) tAATCCGCAAATAACGAAAAACCCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAg < 1:869922/76‑1 (MQ=255) aaTCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGa < 2:1503990/76‑1 (MQ=255) tCCGCAAATAACGTAAAACCCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAat < 1:556967/76‑1 (MQ=255) ccGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAata < 2:52867/76‑1 (MQ=255) cGCAAATAACGTAAAACCCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAata‑t < 2:904111/76‑1 (MQ=255) cGCAAATAACGAAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAata‑t < 1:606531/76‑1 (MQ=255) gCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGGGTTTAGGGAAAGAGCATTTGTCAGAATA‑tt > 2:111193/1‑76 (MQ=255) gCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑tt > 2:102902/1‑76 (MQ=255) gCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGGTTTTTTTTGGGGGGGGGTTTGGGGAAGAGGATTTGTTAGGataat > 2:897840/1‑74 (MQ=16) aaaTAACGTAAAAACCCGCTTCGGCG‑GGGTTTTTTTTTGGGGGGGGTTT‑GGGGAAGAGGAATTTTCAGAAAA‑tttt > 1:666748/1‑75 (MQ=16) aagaaCGAAAAACCCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTaa < 1:714602/73‑1 (MQ=38) aaTAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTaa < 1:1577617/76‑1 (MQ=255) aaTAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTaa < 1:345716/76‑1 (MQ=255) aTAACGAAAAACCCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAg < 1:1228392/76‑1 (MQ=255) tAACGTAAAAACCCGCTTCGGC‑‑GGGGTTTTTTATGGGGGGGGTTTTGGGAAAGAGCATTTGTCAGAATA‑TTTAAgg > 2:262185/1‑76 (MQ=38) aaCGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTTTGGGGGGGGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGa > 2:899761/1‑76 (MQ=255) aaCGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTTTGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGa > 2:1836947/1‑76 (MQ=255) aCGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGaa < 2:216564/76‑1 (MQ=255) gTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTTTGGGGGGGGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAAtt > 2:759105/1‑76 (MQ=255) gTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAAtt > 2:1973937/1‑76 (MQ=255) tAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGGGTTTAGGGAAAGAGCCTTTGTCAGAATA‑TTTAAGGAAttt > 2:270901/1‑76 (MQ=255) aaaaaCCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTc > 2:1682721/1‑76 (MQ=255) aaaaaCCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTc > 2:420438/1‑76 (MQ=255) aaaaCCCGCTTCGGC‑‑GGGGTTTTTTTTGGGGGGGGGTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCt > 2:1939270/1‑76 (MQ=38) aaaCCCGCTTCGGCG‑GGGTTTTTTTTTGGGGGGGGTTTGGGGAAA‑AGGAATTTTCAGGATATTTTAAGGAtttttt > 2:922829/1‑74 (MQ=12) aaaCCCGCTTCGGC‑‑GGGGTTTTTTTTGGGGGGGGTTTTGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTg > 2:1160353/1‑76 (MQ=38) aCCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGaa < 2:1792258/76‑1 (MQ=255) cccGCTTCGGC‑‑GGGGTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAAt > 2:593097/1‑76 (MQ=255) ttCGGC‑‑GGGTTTTTTTTTGGGGGGGGGTTTGGGGAAGAGCATTTGTTAGAATA‑TTTAAGGAATTTTTGAATACTCa > 2:1228671/1‑76 (MQ=255) tCGGCG‑GGGTTTTTTTTTGGGGGGGGTTT‑GGGGAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTTTTTATACTCAt > 1:763260/1‑76 (MQ=255) gggTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCa < 1:932858/76‑1 (MQ=255) gTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCAAt > 1:1149441/1‑76 (MQ=255) gTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCAAt < 2:1556744/76‑1 (MQ=255) gTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCAAt > 2:1551383/1‑76 (MQ=255) ttttttATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCAATGt < 1:138571/76‑1 (MQ=255) ttttttATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCAATGt > 2:1937412/1‑76 (MQ=255) ttttATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCAATGTag < 1:1663711/76‑1 (MQ=255) ttttATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCAATGTag < 1:247135/76‑1 (MQ=255) ttATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCAATGTagag > 1:131639/1‑76 (MQ=255) tATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCAATGTAgaga > 2:1007845/1‑76 (MQ=255) aTGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGAt < 1:1792669/76‑1 (MQ=255) | TCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGC‑‑GGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATA‑TTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGAT > REL606/4169160‑4169307 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |