| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 2,777,419 | 1 | . | A | 72.4% | 58.1 / 26.5 | 29 | coding (3964/4581 nt) | ypjA | adhesin‑like autotransporter |
| Reads supporting (aligned to +/- strand): ref base . (0/8); new base A (14/7); total (14/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.20e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑AAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCTTCAGGAATAT > W3110S.gb/2777363‑2777484 | ccAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑AAAACCGATAtgtg < 1:424075/71‑1 (MQ=255)ccAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑AAAACCGATAtgtg < 1:1329983/71‑1 (MQ=255)ccAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑AAAACCGATAtgtg < 1:2697229/71‑1 (MQ=255)ccAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑AAAACCGATAtgtg < 1:1708478/71‑1 (MQ=255)ccAACTGGCATAGCCGCCCAGAAAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑AAAACCGATAtgtg < 1:3968918/71‑1 (MQ=255) aaCTGGCATAGCCTCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑aaaa < 1:1173120/59‑1 (MQ=255) aaCTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑aaaa < 1:1226366/59‑1 (MQ=255) aaCTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑aaaa < 1:658811/59‑1 (MQ=255) aaCTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑aaaa < 1:3322309/59‑1 (MQ=255) aaCTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑aaaa < 1:1406651/59‑1 (MQ=255) aaCTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑aaaa < 1:3992697/59‑1 (MQ=255) cTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑aaaaac > 1:3167579/1‑57 (MQ=255) cTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑aaaaac > 1:3329012/1‑57 (MQ=255) cTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑aaaaac > 1:4099059/1‑57 (MQ=255) cTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑aaaaac > 1:1333402/1‑57 (MQ=255) cTGGCATAGCCGCCCAGAGAATAACTGCCCAAACTGCCATGTCCTCCGCGATC‑aaaaac > 1:1467145/1‑57 (MQ=255) tAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATgg > 1:1699752/1‑67 (MQ=255) tAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCAGCGATCAAAAACCGATATGTGAATGGGaa > 1:1335253/1‑70 (MQ=255) cgccCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGcc < 1:3806816/71‑1 (MQ=255) cgccCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGcc < 1:852153/71‑1 (MQ=255) cgccCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGcc < 1:1653241/71‑1 (MQ=255) cgccCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑AAAACCGATATGTGAATGGGAATAGcc < 1:309591/70‑1 (MQ=255) cgccCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑AAAACCGATATGTGAATGGGAATAGcc < 1:3578350/70‑1 (MQ=255) gagaATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGcc < 1:207489/65‑1 (MQ=255) gaATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAg > 1:2912763/1‑71 (MQ=255) gaATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAg > 1:825115/1‑71 (MQ=255) gaATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAg > 1:1294442/1‑71 (MQ=255) cacTGCCATGTCCTCCGCGATC‑AAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTa < 1:3338380/71‑1 (MQ=255) tcctccGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCtt > 1:1070055/1‑69 (MQ=255) tcctccGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCtt > 1:1315371/1‑69 (MQ=255) tcctccGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCtt > 1:3937228/1‑69 (MQ=255) tcctccGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCtt > 1:3932029/1‑69 (MQ=255) tcctccGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCtt > 1:217190/1‑69 (MQ=255) tcctccGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCtt > 1:3103562/1‑69 (MQ=255) ccGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCTTCAGGaa > 1:1696325/1‑71 (MQ=255) ccGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCTTCAGGaa > 1:3731286/1‑71 (MQ=255) ccGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCTTCAGGaa > 1:3242673/1‑71 (MQ=255) cgcgATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAAt < 1:201604/59‑1 (MQ=255) cgcgATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCTTCAGGaa < 1:2895512/70‑1 (MQ=255) cgATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCTTCAGGAatat < 1:1961447/71‑1 (MQ=255) | CCAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑AAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCTTCAGGAATAT > W3110S.gb/2777363‑2777484 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |