Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,336,785 | 0 | A | C | 57.1% | 11.1 / 15.3 | 14 | G288G (GGT→GGG) | yfaL | AIDA‑I family autotransporter adhesin YfaL/EhaC |
Reads supporting (aligned to +/- strand): ref base A (5/1); new base C (0/8); total (5/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.00e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGT‑AACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTG > NZ_CP009273/2336699‑2336859 | ccGCCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGT‑AACTTCACTTAAGCCGAGATACATAAAGCCAccc > 1:596675/1‑90 (MQ=255) ttATTCCAAATAACCAGCGTTTTCCCGGCGGCAATATAAAAGGA‑AACTTCACTAAAGCCGAAATACAAAAAGCCCCCCGCCGCAGaggaa < 1:808417/89‑1 (MQ=255) gttttGCCAATAACCACCGTTTTCCCGTCGGCAATAAAAAAGGT‑ACCTTCACTAAACCCGAGATAAAAAAAGCCCCCCGCCGCAGaggaa < 2:615860/87‑1 (MQ=255) ccAATAACCGCCTTTTTCCCGCCGGAAATA‑AAAAGGTCAACTTCCCTTATGCCGATATACATAACCCCCCCCGCCGCAGAGGAAGGACCa < 2:596675/90‑1 (MQ=255) cAATAACCAGCGTTTTCCCGCCGGAAATATCAAAGGT‑AACTTCCCTTAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCAt < 2:10348/90‑1 (MQ=255) cAAAAACCAGCTTTTTCCCGTCGGAAATAAAAAAGGT‑AACTTCACTAAAGCCGAGATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCAt < 1:251598/90‑1 (MQ=255) aaaTAACCAGCTTTTTCCCGCCGGAAATATAAAAGGT‑ACCTTCACTTAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCAt < 2:378456/89‑1 (MQ=255) tttCCGTCGGCAATATCAAAGGT‑AACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCctgtc > 2:374234/1‑88 (MQ=255) cccGTCGGAAATATAAAAGGT‑AACTTCACTTAAGCCGAAATAAAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCtgctgc < 2:786300/89‑1 (MQ=255) gTCGGCAATATCAAAGGT‑AACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCt > 2:823530/1‑90 (MQ=255) ccGAGATACATAAACCCCACCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTc < 1:832374/90‑1 (MQ=255) cGAGATACATAAAGCCACCCGCCGCAGaggaagga > 1:556813/1‑35 (MQ=255) cGAGATACAAAAAGCCACCCGCCGCAGaggaagga < 2:556813/35‑1 (MQ=255) gagaTACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTg > 2:352956/1‑90 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc > 1:651914/1‑36 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 2:651914/36‑1 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATc > 1:198035/1‑54 (MQ=255) | CCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGT‑AACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTG > NZ_CP009273/2336699‑2336859 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |