Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 154,497 | 0 | T | G | 26.2% | 143.2 / 10.8 | 61 | A191A (GCA→GCC) | yadL | predicted fimbrial‑like adhesin protein |
Reads supporting (aligned to +/- strand): ref base T (9/36); new base G (16/0); total (25/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.01e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TACGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA > REL606/154422‑154572 | tACGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGt < 1:68930/76‑1 (MQ=255) gTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCa > 2:873581/1‑76 (MQ=255) tGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGgaga > 2:662915/1‑76 (MQ=255) aGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaa < 2:1824048/76‑1 (MQ=255) tGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAAt < 2:959247/76‑1 (MQ=255) gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAANTGTGGCAGGTGCAGAAAATTTa < 2:258100/76‑1 (MQ=255) gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTa > 1:1639606/1‑76 (MQ=255) gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTa > 1:132656/1‑76 (MQ=255) tGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAcc < 1:469013/76‑1 (MQ=255) tGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTAcc < 1:1359726/76‑1 (MQ=255) gCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTAccc > 1:1222519/1‑76 (MQ=255) cATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCg < 1:115487/76‑1 (MQ=255) cATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCg < 2:1804872/76‑1 (MQ=255) atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGc < 2:1541118/76‑1 (MQ=255) atatatAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTAccccc > 1:214120/1‑74 (MQ=255) tataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCAgg > 1:135972/1‑75 (MQ=255) ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTg > 1:126156/1‑76 (MQ=255) ataGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTg < 2:1827006/76‑1 (MQ=255) aGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGGAGAAAATTTACCCGCCGCGAc > 2:1185001/1‑76 (MQ=255) gTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACa < 1:628212/76‑1 (MQ=255) tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt < 1:36720/76‑1 (MQ=255) ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTg > 1:239332/1‑76 (MQ=255) cTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGGGAAAATTTACCCCCAGCGACATTGc > 2:516904/1‑76 (MQ=255) tGGATCTGTCAGTCTGTTAGTTATAGGTTACTGAAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTgat < 1:1637098/76‑3 (MQ=255) gATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTa < 2:871469/76‑1 (MQ=255) tCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACg < 1:1529458/76‑1 (MQ=255) cTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCGGGGGCATTGCTTATgg > 2:1908261/1‑76 (MQ=255) gTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGGGGAAAAATTAACCCCGGAGACATTGCTTTCGGtt > 2:823513/1‑76 (MQ=255) gTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTGCTTTCgggt > 1:1371447/1‑74 (MQ=255) tCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGttt < 1:1243520/76‑1 (MQ=255) tCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGttt > 1:1767582/1‑76 (MQ=255) cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtttt < 2:964050/76‑1 (MQ=255) cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTAACCGCAGTGACATTGCTTACGGtttt > 1:15112/1‑76 (MQ=255) gTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAt < 2:1303978/76‑1 (MQ=255) tCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtt < 2:289679/76‑1 (MQ=255) tCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTGCTTACGGTtttttt > 1:1882869/1‑76 (MQ=255) cTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAttt < 1:29009/76‑1 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt < 2:1738865/76‑1 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt < 2:1240561/76‑1 (MQ=255) ngttagttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGTGACATTGCTTTCGGGtttttttt > 1:31254/2‑76 (MQ=255) gttagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAttttt < 1:1500483/76‑1 (MQ=255) gttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcg < 1:1313527/76‑1 (MQ=255) gttaTAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTTTTTTTTTcg > 1:1648681/1‑76 (MQ=255) taGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacga < 2:574209/76‑1 (MQ=255) taGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacga < 2:244788/76‑1 (MQ=255) taGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGGGACATTGCTTTCGGTTTTTTTTTTTTttacga > 1:163407/1‑76 (MQ=255) taGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCAGTGACATTGCTTTCGGGTTTTTTTTTTTtgacga > 1:1288129/1‑76 (MQ=255) gTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCa < 2:1130479/76‑1 (MQ=255) gTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCa < 1:980945/76‑1 (MQ=255) ttACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAg < 1:515364/76‑1 (MQ=255) aCTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTTTTTTTTTTGACGACCAGTc > 1:24374/1‑76 (MQ=255) cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg < 2:1732721/76‑1 (MQ=255) gTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAg < 2:1232545/76‑1 (MQ=255) tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGc < 1:829908/76‑1 (MQ=255) aaaTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCa < 2:917192/76‑1 (MQ=255) aaTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCAc < 1:1152032/76‑1 (MQ=255) aTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACt < 2:1551034/76‑1 (MQ=255) aTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACt < 2:1526129/76‑1 (MQ=255) tgtgGCAGGGGCAGAAAATTTACCCCCCGGGACATTTCTTTCGGGTTTTTTTTTTTTTACGACCCGTCGAGCACTc > 2:403298/1‑76 (MQ=255) tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCCCTCAt > 2:527737/1‑76 (MQ=255) tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt > 1:578643/1‑76 (MQ=255) ggCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATc > 1:181438/1‑76 (MQ=255) ggCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATc > 1:181455/1‑76 (MQ=255) gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg < 2:779232/76‑1 (MQ=255) gCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCg > 1:1829511/1‑76 (MQ=255) cAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCgg > 2:227638/1‑76 (MQ=255) tGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGataa < 2:1626476/76‑1 (MQ=255) tGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGataa < 2:1560182/76‑1 (MQ=255) | TACGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA > REL606/154422‑154572 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |