| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,254,698 | 0 | A | C | 70.8% | 34.5 / 19.0 | 24 | L512R (CTG→CGG) | dauA | C4‑dicarboxylic acid transporter DauA |
| Reads supporting (aligned to +/- strand): ref base A (5/2); new base C (10/7); total (15/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.69e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCTTCAAGACGTGACTCCAGGTCCGTG > NZ_CP009273/1254572‑1254829 | tctcGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCt > 2:191478/1‑139 (MQ=255) ctcGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCTc > 1:78914/1‑139 (MQ=255) tACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCTCggg > 2:116858/1‑138 (MQ=255) ataAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACg > 1:119286/1‑139 (MQ=255) ataAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACg > 1:222460/1‑139 (MQ=255) ccAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATc < 1:191478/139‑1 (MQ=255) gTCCCGGGATCGTTTCAATGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTAc < 1:350549/139‑1 (MQ=255) aTCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACGCTTCACAAAACACTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg < 2:44458/139‑1 (MQ=255) tGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACtt < 1:116858/139‑1 (MQ=255) aGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATcac > 1:386029/1‑138 (MQ=255) cgcAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGAccacca < 1:231729/90‑1 (MQ=255) cgcAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGAccacca > 2:231729/1‑90 (MQ=255) cgcAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGc > 2:17498/1‑139 (MQ=255) cAGTGGCTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCt > 1:132313/1‑139 (MQ=255) tGGCTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCTTCa < 2:286347/139‑1 (MQ=255) gCTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCTTc > 1:168068/1‑136 (MQ=255) cTGGAATTCCACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCTTCAAGa > 2:127447/1‑139 (MQ=255) tCCACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCTTCAAGACGTGACt > 2:4622/1‑139 (MQ=255) cACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCTTCAAGAc > 2:374060/1‑131 (MQ=255) cACGTTGCACACGCGCCGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCTTCAAGAc < 1:374060/131‑1 (MQ=255) cacaCGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCTTCAAGACGTGACTCCAGGTCCg > 2:177567/1‑139 (MQ=255) cacGCGCCGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCTTCAAGACGTGACTCCAGGTCCGTg < 1:127447/139‑1 (MQ=255) cAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCTTCAAGACGTGACTCCAGGTc < 2:12431/129‑1 (MQ=255) cAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCTTCAAGACGTGACTCCAGGTc > 1:12431/1‑129 (MQ=255) | TCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCATCCCACTTCAGAATCACAATCCGTTTGCCTTCAAGACGTGACTCCAGGTCCGTG > NZ_CP009273/1254572‑1254829 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |