| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,019,082 | 0 | T | G | 27.4% | 138.5 / 21.4 | 62 | V62G (GTA→GGA) | ycbT | predicted fimbrial‑like adhesin protein |
| Reads supporting (aligned to +/- strand): ref base T (9/36); new base G (17/0); total (26/36) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.22e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTT > REL606/1019032‑1019132 | cGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGt < 1:2833784/51‑1 (MQ=255) gTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTa < 1:6012888/51‑1 (MQ=255) gTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTa < 1:3365705/51‑1 (MQ=255) cAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGt > 1:5731047/1‑51 (MQ=255) aTTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGACAGGGGTAGTTa < 1:4298363/51‑1 (MQ=255) ttttGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTaa < 1:1210125/51‑1 (MQ=255) ttttGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTaa < 1:140583/51‑1 (MQ=255) tttGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTaaa < 1:4945779/51‑1 (MQ=255) ttGGCGTTAAACCCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAg < 1:6056523/51‑1 (MQ=255) tGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGa < 1:6234629/51‑1 (MQ=255) ggCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGAc < 1:1835075/51‑1 (MQ=255) gCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACa < 1:1828005/51‑1 (MQ=255) gCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACa < 1:5523986/51‑1 (MQ=255) gCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACa < 1:6298264/51‑1 (MQ=255) gCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACa < 1:5300106/51‑1 (MQ=255) cGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACaa < 1:156394/51‑1 (MQ=255) cGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACaa < 1:1732679/51‑1 (MQ=255) cGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACaa > 1:4571113/1‑51 (MQ=255) gTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACaaa > 1:7041157/1‑51 (MQ=255) ttAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAAt < 1:3150795/51‑1 (MQ=255) tAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACaaa < 1:1333858/49‑1 (MQ=255) tAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAAtt > 1:2069168/1‑51 (MQ=255) tAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAtt > 1:428239/1‑51 (MQ=255) aaaaCCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACaa < 1:1919773/47‑1 (MQ=255) aaaaCCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAttt > 1:5606507/1‑51 (MQ=255) aaaaCCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAatt > 1:4955308/1‑51 (MQ=255) aaaaCCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGGTAAAGACaaa > 1:4981030/1‑48 (MQ=255) aaaCCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAAtttt < 1:2114480/51‑1 (MQ=255) aaaCCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAAtttt < 1:6181515/51‑1 (MQ=255) aaaCCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGGTAAAGACAAAtttt > 1:3999807/1‑51 (MQ=255) aCCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTAt < 1:2541343/51‑1 (MQ=255) aCCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTAt > 1:2587309/1‑51 (MQ=255) cccTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTAt > 1:6488395/1‑50 (MQ=255) ccTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCa < 1:6211823/51‑1 (MQ=255) ccTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCa < 1:2128306/51‑1 (MQ=255) ccTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCa < 1:4487514/51‑1 (MQ=255) ccTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAtttt > 1:5960160/1‑47 (MQ=255) ggaggaAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTg > 1:2928657/1‑51 (MQ=255) ggaAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGaa < 1:3581623/51‑1 (MQ=255) ggaAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGaa < 1:446497/51‑1 (MQ=255) aagaagCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAAcg > 1:3499270/1‑51 (MQ=255) aagaagCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAAcg > 1:5082152/1‑51 (MQ=255) agaagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcga < 1:1780741/51‑1 (MQ=255) agaagCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGaa > 1:1033162/1‑48 (MQ=255) gaagCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAAcgaa > 1:1497865/1‑51 (MQ=255) aagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaat < 1:1535867/51‑1 (MQ=255) aagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaat < 1:6891728/51‑1 (MQ=255) aaaTACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaa > 1:1905835/1‑51 (MQ=255) aaTACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat < 1:1557344/51‑1 (MQ=255) aTACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATg < 1:2178434/51‑1 (MQ=255) aCGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGat > 1:5625945/1‑51 (MQ=255) aCGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGat < 1:1877570/51‑1 (MQ=255) cGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGatt > 1:6435417/1‑51 (MQ=255) cAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGattatt > 1:2664669/1‑51 (MQ=255) cAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGattatt > 1:6985948/1‑51 (MQ=255) cAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGattatt < 1:4029261/51‑1 (MQ=255) cAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGattatt > 1:4522363/1‑51 (MQ=255) aGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGAttatta < 1:3419043/51‑1 (MQ=255) aGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGAttatta > 1:2749902/1‑51 (MQ=255) ggggTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATtattat < 1:2086515/51‑1 (MQ=255) ggggTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATtattat > 1:1403058/1‑51 (MQ=255) ggggTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATtattat > 1:6928192/1‑51 (MQ=255) ggTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGt < 1:5572140/51‑1 (MQ=255) tAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGttt > 1:6802871/1‑51 (MQ=255) tAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGttt > 1:6162456/1‑51 (MQ=255) | CGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTT > REL606/1019032‑1019132 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |