Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 159,654 | 0 | A | C | 36.1% | 350.9 / 27.8 | 169 | G27G (GGT→GGG) | yadN | predicted fimbrial‑like adhesin protein |
Reads supporting (aligned to +/- strand): ref base A (76/32); new base C (0/61); total (76/93) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.54e-22 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 215 total aligned reads displayed.
ACCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTA > REL606/159579‑159729 | aCCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGacca < 1:1086466/76‑1 (MQ=255) cATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGaccacc < 1:5162883/76‑1 (MQ=255) aTCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccacca > 2:3826050/1‑76 (MQ=255) cTTTATTCCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc < 1:4205223/76‑1 (MQ=255) cTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc > 2:2825630/1‑76 (MQ=255) cTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc > 2:4706149/1‑76 (MQ=255) cTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc < 1:1209256/76‑1 (MQ=255) ttattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatcca < 2:858750/76‑1 (MQ=255) ttattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACAAccatcca > 2:1265522/1‑76 (MQ=255) tattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatccat < 1:3820235/76‑1 (MQ=255) tattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatccat < 1:3684325/76‑1 (MQ=255) tattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatccat > 1:2172968/1‑76 (MQ=255) attaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCata > 2:3521858/1‑76 (MQ=255) ttcccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCatat < 2:333126/73‑1 (MQ=255) tcccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATc < 2:5189957/74‑1 (MQ=255) tcccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATc < 1:1306112/74‑1 (MQ=255) cccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATca < 2:4686262/75‑1 (MQ=255) aCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATca > 2:2106005/1‑76 (MQ=255) ccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcag < 2:4363002/76‑1 (MQ=255) ccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcag < 2:4055585/76‑1 (MQ=255) gccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagca > 1:4394030/1‑76 (MQ=255) ccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcagcag < 2:1100251/76‑1 (MQ=255) ccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcagaag < 2:4088938/76‑1 (MQ=255) cgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGa < 2:2055657/76‑1 (MQ=255) gTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGaa > 1:2535713/1‑76 (MQ=255) gTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGaa > 2:520305/1‑76 (MQ=255) tCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAg > 2:4452148/1‑76 (MQ=255) tCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAg > 1:2888236/1‑76 (MQ=255) cAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGc < 2:4891731/76‑1 (MQ=255) cAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGc < 1:876530/76‑1 (MQ=255) aaaCGCGGGTTTCACAGGTGTTATCAACAACCAACCCCCGGATATTAAACTCCCCACCATCCATATCAGCAGAAGc < 2:1883399/75‑1 (MQ=255) cgcgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGt < 1:3760399/76‑1 (MQ=255) cgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTAc > 2:3421543/1‑76 (MQ=255) gttttCACAGGTTTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCt < 2:677324/74‑1 (MQ=255) ggTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCt < 1:2031878/76‑1 (MQ=255) tttCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGc > 2:1508638/1‑76 (MQ=255) cacaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAAc > 1:2612091/1‑76 (MQ=255) caGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCa > 1:4549737/1‑76 (MQ=255) caGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCa > 1:841619/1‑76 (MQ=255) aGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCAt < 1:2106533/76‑1 (MQ=255) aGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCAt < 1:4813187/76‑1 (MQ=255) gtgtTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTg > 1:1063405/1‑76 (MQ=255) tttATCAACAACCAGCCCACTGATTTTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGcc < 2:4512484/75‑1 (MQ=255) gtTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGcc > 1:4704150/1‑76 (MQ=255) gtTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGcc < 1:3749539/76‑1 (MQ=255) tATCAACACCCGGACCACTGAATTTTAACGGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg < 2:1041769/76‑1 (MQ=255) tATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg > 2:5178717/1‑76 (MQ=255) aTCAACACCCAGACCCCTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGc < 1:3912799/76‑1 (MQ=255) tCAACACCCAGCCCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCa < 2:3022393/76‑1 (MQ=255) ccaacaaCCAGACCACTGATATTAAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCa < 2:5112725/75‑1 (MQ=255) caacaaCCAGCCCACTGTTATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCAt < 2:2324853/76‑1 (MQ=255) caacaaCCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCAt < 2:3058153/76‑1 (MQ=255) caacaaCCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCAt > 2:2768215/1‑76 (MQ=255) acaccCAGACCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa < 1:529589/76‑1 (MQ=255) acaacCAGACCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa < 1:3442413/76‑1 (MQ=255) acaacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa > 2:4292038/1‑76 (MQ=255) caacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg > 2:3108760/1‑76 (MQ=255) caacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg > 2:1451008/1‑76 (MQ=255) caacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg > 2:1370689/1‑76 (MQ=255) aacccAGCCCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg < 1:161604/73‑1 (MQ=255) aaacCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg < 1:4315060/75‑1 (MQ=255) aacCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGa < 1:1018951/76‑1 (MQ=255) aacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGa > 2:3858486/1‑76 (MQ=255) ccAGCCCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGAcc < 1:3588079/76‑1 (MQ=255) gcccACTGATATTAAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGa < 1:4432150/74‑1 (MQ=255) cccACTTATATTTAACCGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGAc < 1:4786612/75‑1 (MQ=255) aCCACTAATTTTTAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGAc < 2:476334/76‑1 (MQ=255) ccACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACa < 2:4930703/76‑1 (MQ=255) cACTGATTTTTAGCTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg < 1:4014140/76‑1 (MQ=255) cACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg < 2:3220516/76‑1 (MQ=255) cACTGATATTAAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg < 2:1854086/76‑1 (MQ=255) aCTGATATTTACCTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAgg < 1:1788276/76‑1 (MQ=255) aCTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAgg > 2:3052089/1‑76 (MQ=255) aCTGATATTTAACGACCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAgg < 2:2696459/76‑1 (MQ=255) cTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAggg < 1:4077442/76‑1 (MQ=255) tGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGc > 2:3197527/1‑76 (MQ=255) gATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa < 1:2478587/76‑1 (MQ=255) gATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa > 2:4563043/1‑76 (MQ=255) gATATTTAACGGACCGCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa < 2:1702643/76‑1 (MQ=255) gATATTAACCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa < 1:3003344/76‑1 (MQ=255) atatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa > 1:3510480/1‑76 (MQ=255) atatTTAACTCCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa < 1:22002/76‑1 (MQ=255) atatTAAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa < 2:250053/76‑1 (MQ=255) tttttAACTACCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa < 1:1324989/74‑1 (MQ=255) ttttAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa < 1:3197746/75‑1 (MQ=255) atTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa > 1:1976856/1‑76 (MQ=255) ttaaaCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAc < 1:590664/73‑1 (MQ=255) ttAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc > 2:2137906/1‑76 (MQ=255) ttAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc > 1:4611117/1‑76 (MQ=255) tAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCt < 2:2457916/76‑1 (MQ=255) accTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTa < 1:3137317/74‑1 (MQ=255) aaCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTa < 1:510173/76‑1 (MQ=255) ccTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTaa < 2:239422/75‑1 (MQ=255) cTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt < 2:4394030/76‑1 (MQ=255) cTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt > 1:5065596/1‑76 (MQ=255) cTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt > 2:561815/1‑76 (MQ=255) tGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAtt > 2:3671645/1‑76 (MQ=255) gACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAttt > 1:1893014/1‑76 (MQ=255) gACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAttt < 2:1941458/76‑1 (MQ=255) ccaccaTCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAttttt > 2:3137347/1‑76 (MQ=255) | ACCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTA > REL606/159579‑159729 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |