Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,648,637 | 0 | A | C | 56.2% | 8.5 / 14.6 | 16 | V209G (GTG→GGG) | yhiD | MgtC/SapB family protein |
Reads supporting (aligned to +/- strand): ref base A (3/4); new base C (0/9); total (3/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.25e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTGATTATCTGGTGTGAATTTCAGGCTTACGGTGAGTCTGGCTACGCTGCCACACAGATTAGCTAATTGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATCAATGGCGACC > NZ_CP009273/3648556‑3648714 | atTGATTATCTGGTGTGAATTTCAGGCTTACGGTGAGTCTGGCTACGCTGCCACACAGATTAGCTAATTGAAACGCCTTTCACCCCTGcc > 1:268766/1‑90 (MQ=255) gAATTTCAGGCTTACGGTGAGTCTGGCTACGCTGCCCCACAGATTAGCTAATTGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATc < 2:341484/90‑1 (MQ=255) cAGGCTTACGGTGAGTCTGGCTCCGCTCCCACACAGATTAGCTAATTGAAACGCCTTCCCCCCCTGCCATACCTTTTAATAATCGCAACa < 1:331622/90‑1 (MQ=255) cAGGCTTACGGTGAGTCTGGCTACGCTGCCACACAGATTAGCTAATTGAAACGCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACa < 2:294854/90‑1 (MQ=255) cTTACGGTGAGTCTGGCTACGCTGCCACACAGATTAGCTAATTGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGtc < 2:76051/90‑1 (MQ=255) cTGCCACACAGATTAGCTAATTGAAACGCCTTTCCCCCCTGCCATACC‑TTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGa < 1:96248/90‑1 (MQ=255) cacacaGATTAGCTAATTGAAACGCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGc < 1:413847/90‑1 (MQ=255) acaGATTAGCTAATTGAACCCCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTgg < 1:213760/90‑1 (MQ=255) acaGATTAGCTAATTGAAACGCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTgg < 2:259185/90‑1 (MQ=255) ttAGCTAATTGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATc > 2:104112/1‑90 (MQ=255) aGCTAATTGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATCaa < 1:104112/90‑1 (MQ=255) ttGAAACGCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCAt < 1:436185/69‑1 (MQ=255) ttGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCAt > 2:436185/1‑69 (MQ=255) tGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATCAATGGCGAc < 1:355678/90‑1 (MQ=255) gAAACGCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATCAATGGCGAcc < 1:148634/90‑1 (MQ=255) gAAACGCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATCAATGGCGAcc < 2:369017/90‑1 (MQ=255) | ATTGATTATCTGGTGTGAATTTCAGGCTTACGGTGAGTCTGGCTACGCTGCCACACAGATTAGCTAATTGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATCAATGGCGACC > NZ_CP009273/3648556‑3648714 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |