Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 154,497 | 0 | T | G | 28.8% | 114.5 / 10.9 | 52 | A191A (GCA→GCC) | yadL | predicted fimbrial‑like adhesin protein |
Reads supporting (aligned to +/- strand): ref base T (4/33); new base G (15/0); total (19/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.65e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA > REL606/154423‑154572 | aCGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATANGTTACTGTAAATGTGGCAGGTg < 2:232918/76‑1 (MQ=255) cGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGc < 1:762392/76‑1 (MQ=255) gTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCa < 1:1249727/76‑1 (MQ=255) gTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGcg > 1:249264/1‑75 (MQ=255) ttGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAg < 1:51872/76‑1 (MQ=255) ttGAGTGGGATGCATATATAGTTCCTGGACCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAg > 1:968138/1‑76 (MQ=255) gAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGaa > 1:54823/1‑76 (MQ=255) aGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGaaa > 1:221237/1‑76 (MQ=255) gTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGaaaa > 1:79376/1‑76 (MQ=255) ggATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAatt > 1:1466424/1‑76 (MQ=255) gATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTa < 1:143961/76‑1 (MQ=255) taGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGa < 2:64462/76‑1 (MQ=255) taGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCAGGGa > 1:383240/1‑76 (MQ=255) aGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGAc < 1:1325011/76‑1 (MQ=255) ttCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAAATTACCCGCCGTGAcct > 1:350834/1‑74 (MQ=255) ttCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGTGGCAt > 1:1407710/1‑76 (MQ=255) tCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACAtt < 1:1093185/76‑1 (MQ=255) ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTg < 2:507006/76‑1 (MQ=255) ccTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTg < 1:16053/76‑1 (MQ=255) cTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGTGACATTGc > 2:1032910/1‑76 (MQ=255) tGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCCCAGTGACATTTCt > 1:510676/1‑76 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCtt < 2:1666831/76‑1 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCtt < 1:1068744/76‑1 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCtt < 1:15693/76‑1 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCCGTGACATTGttt > 1:480818/1‑76 (MQ=255) ggATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTGCtt > 2:326913/1‑76 (MQ=255) tCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTTCg > 2:1408056/1‑76 (MQ=255) tCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAAATTACCCGCAGTGACATTGCTTACg > 1:1026348/1‑76 (MQ=255) tGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTTCGGt > 1:1342803/1‑76 (MQ=255) tCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGttt > 2:412129/1‑76 (MQ=255) cAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGtttt < 2:1621626/76‑1 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTa < 2:90625/76‑1 (MQ=255) gTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAt < 2:1667802/76‑1 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttt < 1:1226131/76‑1 (MQ=255) tGTTAGTTATAGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTAACCCCCGGGAAATTGCTTTCGGTtttttttt > 1:181416/1‑76 (MQ=255) ttagttaTAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTAtttttt < 2:1562724/76‑1 (MQ=255) tataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgac < 1:1507127/76‑1 (MQ=255) ataGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTcgacg < 2:21759/76‑1 (MQ=255) ataGGTTACTGTAAATGTGGCAGGGGCAGAAAATTTAACCGCAGTGACATTGCTTTCGGTTTTTTTTTTTTtgacg > 1:1350414/1‑76 (MQ=255) aGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCgacgac < 1:507155/76‑1 (MQ=255) ggTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGAcc < 2:245425/76‑1 (MQ=255) gTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCa < 2:1674048/76‑1 (MQ=255) cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg < 2:115027/76‑1 (MQ=255) cTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCg < 2:249264/76‑1 (MQ=255) cTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTTTTTTTTTCGACGACCCGTCg > 1:123169/1‑76 (MQ=255) cTGTAAATGTGGCAGGGGCAGAAAATTTACCCGCAGGGACATTGCTTACGGGTTTATTTTTTTCGACGACCAGTCg > 1:1325704/1‑76 (MQ=255) tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGc < 1:292315/76‑1 (MQ=255) tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGc < 1:1172371/76‑1 (MQ=255) tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGc < 1:1172378/76‑1 (MQ=255) tAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGc < 2:1161469/76‑1 (MQ=255) aaaTGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCa < 2:158893/76‑1 (MQ=255) tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt > 1:1497315/1‑76 (MQ=255) tgGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCAt < 1:975366/76‑1 (MQ=255) aGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGa < 1:555800/76‑1 (MQ=255) aGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGa > 1:929072/1‑76 (MQ=255) gTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGata < 1:933623/76‑1 (MQ=255) tGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGataa > 2:835872/1‑76 (MQ=255) | ACGTTGAGTGGGATGCATATATAGTTCCTGGATCTGTCAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCAGGTGCAGAAAATTTACCCGCAGTGACATTGCTTACGGTTTTATTTTTTTCGACGACCAGTCGAGCACTCATCGGATAA > REL606/154423‑154572 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |