Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 159,654 | 0 | A | C | 35.5% | 537.8 / 46.0 | 262 | G27G (GGT→GGG) | yadN | predicted fimbrial‑like adhesin protein |
Reads supporting (aligned to +/- strand): ref base A (105/64); new base C (0/93); total (105/157) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.71e-28 | |||||||||||
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 340 total aligned reads displayed.
ACCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTA > REL606/159579‑159729 | aCCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGacca < 1:998829/76‑1 (MQ=255) ccATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGaccac < 2:7679289/76‑1 (MQ=255) aTCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccacca > 2:8053799/1‑76 (MQ=255) aTCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccacca > 2:4914910/1‑76 (MQ=255) tCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccat < 1:4101673/76‑1 (MQ=255) tCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccat < 2:8709483/76‑1 (MQ=255) cTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc > 2:1222301/1‑76 (MQ=255) cTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatc > 2:114982/1‑76 (MQ=255) ttattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatcca > 1:8938351/1‑76 (MQ=255) tattcCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatccat < 2:6394708/76‑1 (MQ=255) tattaCCGCCGTCAACGCGGGTTTCACAGGTTTTATCAACAACCAGACCACTGATATTTAACTGACCAccatccat < 1:5770448/76‑1 (MQ=255) tattaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCAccatccat > 1:4550874/1‑76 (MQ=255) tttaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCata < 2:2698669/75‑1 (MQ=255) attaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCata > 2:3094891/1‑76 (MQ=255) ttaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAGCACCATCCatat < 2:3012536/76‑1 (MQ=255) ttaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCatat < 2:8859940/76‑1 (MQ=255) ttaCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCatat > 1:1045414/1‑76 (MQ=255) taCCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATc > 2:4335427/1‑76 (MQ=255) cgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcagc < 2:5600110/76‑1 (MQ=255) gccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagca > 1:6258146/1‑76 (MQ=255) ccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATcagcag < 1:7150017/76‑1 (MQ=255) ccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATcagcag > 2:8597111/1‑76 (MQ=255) tCAACGCGGGTTTCACAGGTTTTATCAACAACCAGCCCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAg < 2:748993/76‑1 (MQ=255) tCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAg < 2:8438818/76‑1 (MQ=255) aCGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAg < 1:693444/76‑1 (MQ=255) cgcgGTTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGt < 2:8179613/76‑1 (MQ=255) cgcgGGTTTCCAAGGTGTTACCAACATCCAGACCGGTGATTTTTAGCTACCCCCCATCCATATCAGCAGAAGCAGt < 1:6472007/76‑1 (MQ=255) cgcgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGt < 1:1669872/76‑1 (MQ=255) gcgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTa < 2:7254661/76‑1 (MQ=255) gcgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTa > 2:3241212/1‑76 (MQ=255) cgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTAc > 2:1983220/1‑76 (MQ=255) cgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTAc > 1:5697215/1‑76 (MQ=255) gTTTCACAGGTTTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTg < 2:2640862/76‑1 (MQ=255) cacaGGTTTTAACAAGACCAAGCCCCCTGATATTTAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAAc < 2:1433823/76‑1 (MQ=255) acaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAAcc > 2:3630659/1‑76 (MQ=255) acaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAAcc > 1:2191187/1‑76 (MQ=255) caGGTGTTATCAAAAACCAGCCCACTGATATTTACCTGTCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCa < 2:7084707/76‑1 (MQ=255) gggTTTTATCAACAACCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCAt < 2:5249224/75‑1 (MQ=255) aGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCAt > 2:2277367/1‑76 (MQ=255) aGGAGTTATCAACAACCAGACCACGGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCAt < 1:6676275/76‑1 (MQ=255) ggTGTTATCAACACCCAGACCCCTAATTTTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCAtt < 2:412481/76‑1 (MQ=255) ttgtAATCAACAACCAGACCACTGATATTAAACTCACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTg < 2:246997/75‑1 (MQ=255) gtgtTATCAACAACCAGCCCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTg < 1:2579606/76‑1 (MQ=255) ttATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCa > 2:7432042/1‑76 (MQ=255) taacaacacCCAGCCCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg < 2:5456333/73‑1 (MQ=255) tATCAACACCCAGCCCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg < 2:422761/76‑1 (MQ=255) tATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg > 2:7952481/1‑76 (MQ=255) tATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg > 1:3954264/1‑76 (MQ=255) aaTCAAAACCCGGCCCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAg < 1:5904055/75‑1 (MQ=255) aTCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGc > 1:266157/1‑76 (MQ=255) aTCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGc > 2:4742422/1‑76 (MQ=255) caacacCCAGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCAt < 1:2103128/76‑1 (MQ=255) aacaacCAGACCACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATa < 1:8083115/76‑1 (MQ=255) aacaacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATa < 1:1506361/76‑1 (MQ=255) acaccCAGCCCACTGATATTTACCTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa < 1:3989665/76‑1 (MQ=255) acaacCAGCCCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATaa < 1:1288181/76‑1 (MQ=255) cacccAGCCCACTGATATTTAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg < 1:3328812/73‑1 (MQ=255) caacCAGACCACTGATTTTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg < 1:2331555/76‑1 (MQ=255) caacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAg > 2:3441858/1‑76 (MQ=255) aacCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGa > 1:1764311/1‑76 (MQ=255) acCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGAc > 1:7637432/1‑76 (MQ=255) ccAGCCCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGAcc < 1:5009498/76‑1 (MQ=255) cAGCCCACTAATTTTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCa < 2:4612453/76‑1 (MQ=255) cAGACCCCTTATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCa < 2:162778/76‑1 (MQ=255) agcccACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg < 1:1222122/73‑1 (MQ=255) agcccACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg < 1:8652524/73‑1 (MQ=255) aGACCACTGATATTTAACTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg < 1:7421837/76‑1 (MQ=255) aGACCACTGATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAg < 2:599379/76‑1 (MQ=255) cccACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGAc < 2:3750501/75‑1 (MQ=255) cACTGATATTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAg < 1:4956317/76‑1 (MQ=255) cttGATTTTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAgg < 1:6850821/74‑1 (MQ=255) aCTGATTTTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAgg < 2:447342/76‑1 (MQ=255) tGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGc < 2:1755247/76‑1 (MQ=255) tatatTAAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa < 1:8434557/75‑1 (MQ=255) gATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa < 1:8740718/76‑1 (MQ=255) gATATTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCa < 1:2858962/76‑1 (MQ=255) atatTTAACTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa < 1:3110775/76‑1 (MQ=255) atatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa < 2:1804899/76‑1 (MQ=255) atatTTAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa < 1:4754545/76‑1 (MQ=255) aattttAACTGCCCACCACCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaa < 1:6147138/73‑1 (MQ=255) tatTTAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa < 1:4194304/76‑1 (MQ=255) tatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa > 2:5342537/1‑76 (MQ=255) tatTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaa > 1:4313737/1‑76 (MQ=255) ttttAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa < 1:5161818/75‑1 (MQ=255) atTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa > 1:8044402/1‑76 (MQ=255) atTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa > 1:5849208/1‑76 (MQ=255) atTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCaaaa > 2:612548/1‑76 (MQ=255) ttgaaCTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAc < 2:5838360/73‑1 (MQ=255) ttAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc < 2:6881617/76‑1 (MQ=255) ttAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc < 1:783483/76‑1 (MQ=255) ttAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAAcc > 2:6665660/1‑76 (MQ=255) tAACTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCt < 1:5354101/76‑1 (MQ=255) tAACGGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCt < 1:378906/76‑1 (MQ=255) aCTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTaa < 1:3379485/76‑1 (MQ=255) cggACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt < 1:8872969/74‑1 (MQ=255) cTGCCCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt < 2:2892182/76‑1 (MQ=255) cTGCCCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt < 2:509139/76‑1 (MQ=255) cTGACCCCCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAAt < 1:3561524/76‑1 (MQ=255) accaTCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTa > 1:1988236/1‑76 (MQ=255) accaTCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTa > 2:4023879/1‑76 (MQ=255) | ACCATCTTTATTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCACCATCCATATCAGCAGAAGCAGTACCTGCAACCATTGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTA > REL606/159579‑159729 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |