Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 159,654 | 0 | A | C | 35.3% | 517.7 / 44.4 | 258 | G27G (GGT→GGG) | yadN | predicted fimbrial‑like adhesin protein |
Reads supporting (aligned to +/- strand): ref base A (101/66); new base C (0/91); total (101/157) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.09e-27 | |||||||||||
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 335 total aligned reads displayed.
ACCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTA > REL606/159579‑159729 | aCCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGacca > 2:1929715/1‑76 (MQ=255) aCCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGacca < 2:3504095/76‑1 (MQ=255) aCCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGacca > 1:3467115/1‑76 (MQ=255) ccATCTTTATTCCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGaccac < 1:7028457/76‑1 (MQ=255) aTCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccacca > 1:4514414/1‑76 (MQ=255) aTCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccacca > 2:6658792/1‑76 (MQ=255) aTCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccacca < 2:2988422/76‑1 (MQ=255) cTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc > 2:7581142/1‑76 (MQ=255) tttATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatcc > 1:1541996/1‑76 (MQ=255) ttattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACAAcaatcca > 1:1429682/1‑76 (MQ=255) ttaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCatat > 1:1888467/1‑76 (MQ=255) ttaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCatat < 2:3900358/76‑1 (MQ=255) tcccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATc < 2:4893323/74‑1 (MQ=255) taCCGCCGTCAACGCGGGTTTCACAGGTGGTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATc > 2:6410994/1‑76 (MQ=255) aCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATca > 2:1536804/1‑76 (MQ=255) aCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATca > 2:6288909/1‑76 (MQ=255) ccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcag < 2:233295/76‑1 (MQ=255) cgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagc > 1:164271/1‑76 (MQ=255) cgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagc > 2:6936785/1‑76 (MQ=255) cgccgTCAACGCGGGTTTCACAGGTGTTATAAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagc < 2:7196466/76‑1 (MQ=255) ccgTCAAGGCGGGTTTCACAGGTGTTATCAACAACCAGCCCACTGATATTTAACTGACCCCCATCCATATcagcag < 2:877050/76‑1 (MQ=255) ccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcagcag < 2:5896887/76‑1 (MQ=255) ccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagcag > 1:2227095/1‑76 (MQ=255) gTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGCCCACTGATATTTAACTGACCCCCATCCATATCAGCAGaa < 2:1055413/76‑1 (MQ=255) gTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGaa > 1:868692/1‑76 (MQ=255) gTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGaa < 1:5349550/76‑1 (MQ=255) tCAACGCGGTTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAg < 2:4645960/76‑1 (MQ=255) tCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAg < 2:426691/76‑1 (MQ=255) tCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAg < 2:5471273/76‑1 (MQ=255) tCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAg > 2:5799759/1‑76 (MQ=255) cAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGc > 2:2615658/1‑76 (MQ=255) aaCGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCa > 1:7623928/1‑76 (MQ=255) aCGCGGTTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAg < 1:4398756/76‑1 (MQ=255) aCGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAg < 2:7766900/76‑1 (MQ=255) cgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTAc > 2:5846453/1‑76 (MQ=255) gttttCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCt < 1:2915785/74‑1 (MQ=255) ggTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCt < 1:5959966/76‑1 (MQ=255) ggTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCt < 1:2048579/76‑1 (MQ=255) gTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCTTATCAGCAGAAGCAGTACCTg < 1:6188221/76‑1 (MQ=255) tCACAGGTGTTATCAACAACCAGACCACTGATTTTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCaa < 1:5203440/76‑1 (MQ=255) cacaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAAc > 2:5740940/1‑76 (MQ=255) acaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAAcc < 2:5249852/76‑1 (MQ=255) aGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCAt > 2:1710789/1‑76 (MQ=255) ggTGTTATCAACAACCAGCCCACTGATATTTAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCAtt < 1:6587452/76‑1 (MQ=255) gtgtTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTg > 1:6930442/1‑76 (MQ=255) tttATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGcc > 2:7885220/2‑76 (MQ=255) gtTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGcc > 1:2603290/1‑76 (MQ=255) ttATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCa > 1:3867953/1‑76 (MQ=255) ttATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCa > 1:6017598/1‑76 (MQ=255) tataAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg < 2:4094051/76‑1 (MQ=255) aacaccCACACCACTGATTTTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATa < 2:5034023/76‑1 (MQ=255) acaacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa < 2:708269/76‑1 (MQ=255) caacCAGCCCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg < 1:1017441/76‑1 (MQ=255) acccAGCCCACTGATATTTACCGGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGa < 1:5899696/74‑1 (MQ=255) aacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGa > 2:1804616/1‑76 (MQ=255) ccAGACCCCTAATATTAACCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGAcc < 1:3323172/76‑1 (MQ=255) ccAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGAcc > 2:3579443/1‑76 (MQ=255) cAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCa > 1:6818687/1‑76 (MQ=255) aGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg < 2:7426655/76‑1 (MQ=255) aGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg > 1:4132325/1‑76 (MQ=255) aGACCACTGATATTTAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg < 1:4847341/76‑1 (MQ=255) ccACTGATTTTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACa < 2:657083/76‑1 (MQ=255) ccACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACa < 1:562145/76‑1 (MQ=255) ccACTGATATTAAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACa < 2:7528705/76‑1 (MQ=255) cACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg < 2:1348651/76‑1 (MQ=255) cACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg < 1:7499809/76‑1 (MQ=255) cACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg < 1:7374297/76‑1 (MQ=255) cACTGATATTAAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg < 1:4018636/76‑1 (MQ=255) cACTAATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg < 1:72130/76‑1 (MQ=255) aCTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAgg < 2:3696097/76‑1 (MQ=255) cTGATTTTTAACTAACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAggg < 2:2425168/76‑1 (MQ=255) cTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAggg < 1:1149842/76‑1 (MQ=255) cTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAggg < 1:3599516/76‑1 (MQ=255) cTGATATTTAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAggg < 1:5627502/76‑1 (MQ=255) cTGATATTAACCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAggg < 2:2080232/76‑1 (MQ=255) gATTTTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa < 2:4815715/76‑1 (MQ=255) atatTTACCCGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa < 1:4994943/76‑1 (MQ=255) atatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa > 1:1057349/1‑76 (MQ=255) atatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa > 1:5764932/1‑76 (MQ=255) atatTAACCTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa < 2:4797064/76‑1 (MQ=255) tatTTACCGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa < 1:2141799/76‑1 (MQ=255) tatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa > 1:7950235/1‑76 (MQ=255) tatTAAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa < 2:4422220/76‑1 (MQ=255) tatTAAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa < 2:5625106/76‑1 (MQ=255) tatTAAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa < 1:624989/76‑1 (MQ=255) atTTAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa < 2:3287199/76‑1 (MQ=255) tttAACTCCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAc < 2:1093125/76‑1 (MQ=255) ttAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc < 1:654091/76‑1 (MQ=255) ttAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc < 1:2757826/76‑1 (MQ=255) ttAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa < 1:2430981/57‑1 (MQ=255) taccTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCt < 1:6281837/73‑1 (MQ=255) tAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCt < 1:5432164/76‑1 (MQ=255) tAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCt < 2:5983378/76‑1 (MQ=255) aaCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTa < 1:963172/76‑1 (MQ=255) aaCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTa < 1:7739363/76‑1 (MQ=255) aCTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTaa < 1:3119386/76‑1 (MQ=255) aCTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTaa > 1:3908110/1‑76 (MQ=255) cacccccATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAttt < 2:214565/75‑1 (MQ=255) accaccATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAtttt < 2:7710699/76‑1 (MQ=255) accaTCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTa > 2:805390/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: ‑ |