Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,086,499 | 0 | C | A | 54.5% | ‑0.9 / 11.6 | 11 | V100V (GTG→GTT) | rhaM | L‑rhamnose mutarotase |
Reads supporting (aligned to +/- strand): ref base C (1/4); new base A (6/0); total (7/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.52e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCA > NZ_CP009273/4086422‑4086573 | tttAATTTTTCAGATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACATCTTTCAGc > 2:117107/1‑90 (MQ=255) tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAAATTTTTTCAGc > 2:3134/1‑90 (MQ=255) ttttCAGATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTa < 1:117107/90‑1 (MQ=255) ttttCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTa < 1:3134/90‑1 (MQ=255) cATCTCTGGCACATTGGGCAATTACGGCAGGGGAAAAAATTCTTTCAGCTCCCTACTCACCCGGGCGTTTTTCGGGGTTGCGGGGATaaa > 2:170095/1‑89 (MQ=255) ttGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCGCTACTCCCCCGGCTGTTTTTCGGGGTAGGGGGGATAACATCGGGCatatat > 2:250803/1‑90 (MQ=255) ttGGGCAATTACGGCAGGGAAAAAACTTCTTGCAGCTCGCTACTCACCCGGCTGTTTTCCGGGGTTGCGGGGATAACAACGGTCatatat > 2:96358/1‑90 (MQ=255) tGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATAtt < 2:280800/90‑1 (MQ=255) tGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCGCTACTCACCCGGCTGTTTTCCGGGGTAGGGGGGATAACAACGGGCATAAAtt > 1:296301/1‑90 (MQ=255) aTTACGGCAGGGAAAACACTTCTTTCAGCTCGCTACTCACCGGGCTGTTTTCCGGGGTAGCGGGGATAAAATCGGGCATATATTTccccc > 1:282579/1‑90 (MQ=255) ttACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTccacca < 1:250803/90‑1 (MQ=255) | TTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCA > NZ_CP009273/4086422‑4086573 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |