Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 305,236 | 0 | A | C | 20.0% | 57.7 / 10.2 | 25 | F528V (TTT→GTT) | ecpD | fimbrial adhesin EcpD |
Reads supporting (aligned to +/- strand): ref base A (20/0); new base C (0/5); total (20/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.88e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAA > NC_000913/305176‑305291 | aGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAg > 1:807945/1‑65 (MQ=255) aGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAg > 1:2025352/1‑65 (MQ=255) ttAATGTTACGCCACGTCGCCTGAACATAAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGcc > 1:2097159/1‑65 (MQ=255) tGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTa > 1:81496/1‑65 (MQ=255) gTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTat > 1:12087/1‑65 (MQ=255) gTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTat > 1:3062853/1‑65 (MQ=255) gTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTat > 1:1858317/1‑65 (MQ=255) gTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTat > 1:1791295/1‑65 (MQ=255) gCCCCGTCGCCGGAACATGAGTCTCCCGTAAAGCGCTCATTTCGCCAACCCAGCCGTTATCATCt < 1:2937204/65‑1 (MQ=255) gCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCt > 1:1951787/1‑65 (MQ=255) gCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCt > 1:695831/1‑65 (MQ=255) gCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCt > 1:3053659/1‑65 (MQ=255) gCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCt > 1:273749/1‑65 (MQ=255) cACGTCGCCTGAACATAAACCTCTCCTGAAGCGCTGCCTCCGCCAACCCAGCCGTTATCATCTAc < 1:2945284/65‑1 (MQ=255) gTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGt > 1:2724553/1‑65 (MQ=255) aCATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGa > 1:326028/1‑65 (MQ=255) aCATAAATCTCTCTGGAAGCGCTGACTCCCCAAACCCAGCCGTTATCATCTACCGTACGCAGCGa < 1:2118973/65‑1 (MQ=255) tcctctcCTAAGCCGCTGTCGCCGACAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGc < 1:847016/63‑1 (MQ=255) tcctctcCGGAGGCGCTGCCTCCCCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGc < 1:2536389/63‑1 (MQ=255) aTCTCCCCCGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGc < 1:420241/65‑1 (MQ=255) tctctcCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCg > 1:1506205/1‑65 (MQ=255) ctctcCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGt > 1:556169/1‑65 (MQ=255) ccTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATc > 1:137814/1‑65 (MQ=255) aGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTg > 1:2968363/1‑65 (MQ=255) aGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTg > 1:2667311/1‑65 (MQ=255) aGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGCACGCAGCGCAATGGCGTTATCCATTg > 1:1962774/1‑65 (MQ=255) cgcTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGa > 1:1416947/1‑65 (MQ=255) cTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaa > 1:1689719/1‑65 (MQ=255) | AGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAA > NC_000913/305176‑305291 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |