| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 3,998,787 | 0 | A | C | 31.2% | 125.4 / 11.8 | 65 | T31P (ACC→CCC) | ubiB | putative ubiquinone biosynthesis protein UbiB |
| Reads supporting (aligned to +/- strand): ref base A (27/17); new base C (0/20); total (27/38) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.74e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTACGGCGCCTATATTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACTTTTAGGTGAG > REL606/3998712‑3998861 | gTACGGCGCCTATATTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCa < 2:1391814/76‑1 (MQ=255) tACGGCGCCTATATTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCAc > 2:1315718/1‑76 (MQ=255) ggCGCCTATATTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCt > 1:281158/1‑76 (MQ=255) ggCGCCTATATTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCt < 1:1622578/76‑1 (MQ=255) cTATATTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGc > 2:916245/1‑76 (MQ=255) atatTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTa < 1:885968/76‑1 (MQ=255) tatTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCCCCCTGCCGCTAc < 2:420535/76‑1 (MQ=255) tatTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTAc < 2:949312/76‑1 (MQ=255) atTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACGGATCCCCAAAATGCGCATCACCCTGCCGCTACg < 2:1409349/76‑1 (MQ=255) ttCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGc > 2:1586812/1‑76 (MQ=255) tcatcaTTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCt > 1:1628000/1‑76 (MQ=255) atcatTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTAt < 1:1018200/76‑1 (MQ=255) tcatTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATg < 1:1196630/76‑1 (MQ=255) catTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATgg > 1:551598/1‑76 (MQ=255) ttCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCg > 2:1451291/1‑76 (MQ=255) cGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGAt < 1:492618/76‑1 (MQ=255) cGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGAt < 1:154555/76‑1 (MQ=255) cGCACTTTTTTAAGCTACGGACTTGATGAACTGACCCCCAAAATGCGTACCCCCCTGCCGCTACGGCTATGGCGAt < 1:1286123/76‑1 (MQ=255) gCACTTTTTTAAGCTACGGACTTGATGATCTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATa < 1:1438949/76‑1 (MQ=255) cACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATCCGTACCCCCCTGCCGCTACGGCTATGGCGATAc < 2:436499/76‑1 (MQ=255) aCTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACt > 1:1608601/1‑76 (MQ=255) cTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTc > 2:1334485/1‑76 (MQ=255) ttttttAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCa > 1:1764707/1‑76 (MQ=255) tttttAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCat > 1:1331388/1‑76 (MQ=255) ttAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCattat > 2:1814075/1‑76 (MQ=255) tAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCattatt > 1:1176401/1‑76 (MQ=255) gCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTg > 2:1884736/1‑76 (MQ=255) cTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTgg > 1:1338335/1‑76 (MQ=255) tACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGa > 2:253628/1‑76 (MQ=255) aCGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGAt > 1:771788/1‑76 (MQ=255) gACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGcc < 1:2073486/76‑1 (MQ=255) gcTTGATGAACTGACCCCAAAAATGCGTACCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCa < 2:1474021/75‑1 (MQ=255) ttGATGAACTGACCCCCAAAATGCGTATCGCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCaaa < 1:912055/76‑1 (MQ=255) gatgaACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATc < 1:1164686/76‑1 (MQ=255) gatgaACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATc < 2:1671882/76‑1 (MQ=255) atgaACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCg < 1:399150/76‑1 (MQ=255) atgaACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCg < 2:516990/76‑1 (MQ=255) tgaACTGATCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCgg < 2:91210/76‑1 (MQ=255) tgaACTGATCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCgg < 2:482936/76‑1 (MQ=255) gaACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGc < 2:378437/76‑1 (MQ=255) gaACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGc < 1:57720/76‑1 (MQ=255) gaACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGc > 1:1318788/1‑76 (MQ=255) aCTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCAt > 1:1247481/1‑76 (MQ=255) aCTGATCCCAAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCAt < 1:1444453/76‑1 (MQ=255) cTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATa > 1:271809/1‑76 (MQ=255) cTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATa > 2:1013892/1‑76 (MQ=255) tGATCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATaa < 1:1586812/76‑1 (MQ=255) taTCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATaaa < 2:1841602/75‑1 (MQ=255) gATCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATaaa < 2:1885579/76‑1 (MQ=255) gATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATaaa > 2:382003/1‑76 (MQ=255) gATCCCCAAAATGCGTACCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATaaa < 2:1841626/76‑1 (MQ=255) aTCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAg < 2:1256978/76‑1 (MQ=255) tCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGa < 2:1037580/76‑1 (MQ=255) tCCCCAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGa < 1:2035966/76‑1 (MQ=255) ccAAAATGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACaa < 1:1443178/76‑1 (MQ=255) ccAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACaa < 1:2021978/76‑1 (MQ=255) ccAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACaa < 1:2021984/76‑1 (MQ=255) aaaaTGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAAc < 1:1160782/76‑1 (MQ=255) aaaaTGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAAc > 1:1414203/1‑76 (MQ=255) aaaTGGGAATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACt < 2:1877627/76‑1 (MQ=255) aaaTGCGTATCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACt < 2:1877615/76‑1 (MQ=255) aaTGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACtt < 1:1971345/76‑1 (MQ=255) aaTGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACtt > 2:1063748/1‑76 (MQ=255) gCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACTTTTa > 1:918255/1‑76 (MQ=255) cGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACTTTTAg > 2:1453096/1‑76 (MQ=255) aTCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACTTTTAGGTg < 1:1542479/76‑1 (MQ=255) aTCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACTTTTAGGTg < 1:1841684/76‑1 (MQ=255) aTCCCCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACTTTTAGGTg < 1:645272/76‑1 (MQ=255) cACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACTTTTAGGTGAg > 2:812518/1‑76 (MQ=255) | GTACGGCGCCTATATTTCATCATTCGCACTTTTTTAAGCTACGGACTTGATGAACTGATCCCCAAAATGCGTATCACCCTGCCGCTACGGCTATGGCGATACTCATTATTCTGGATGCCAAATCGGCATAAAGACAAACTTTTAGGTGAG > REL606/3998712‑3998861 |
| 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: ‑ |