Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,086,499 | 0 | C | A | 47.4% | 6.4 / 14.0 | 19 | V100V (GTG→GTT) | rhaM | L‑rhamnose mutarotase |
Reads supporting (aligned to +/- strand): ref base C (4/6); new base A (9/0); total (13/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCA > NZ_CP009273/4086419‑4086582 | atGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGc < 1:602695/90‑1 (MQ=255) tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGc > 2:413677/1‑90 (MQ=255) ttttCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGGAAAAACTTCTTTGAGCTCGCTa > 1:4812/1‑90 (MQ=255) tttCAAATCGAAATTTAAAATATTGTGCCGGAGGCGTCTCTGGCACATTGGGCAATTACGGCAGGGAAAACACCTCTTGCAGCTCGCCAc > 2:290904/1‑90 (MQ=255) gCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCg < 2:631635/90‑1 (MQ=255) gCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCg < 1:413625/90‑1 (MQ=255) cATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTTTTTCAGCTCGCTCCTCCCCCGGCTGTTTTTCGGGGTTGCGGGCATaaa > 1:631635/1‑89 (MQ=255) ctctGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATc < 2:202344/90‑1 (MQ=255) ttGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCatatat > 2:255677/1‑90 (MQ=255) ttGGGCAATTACGGCAGGGAAAAAACTTCTTTCCGCTCCCCACTCACCCGGCTGTTTTTCGGGGTAGCGGGCATAACATCGGTCatatat > 1:323099/1‑90 (MQ=255) ttGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCGCTACTCCCCCGGCCGTTTTCCCGGGTAGCGGGGATAAAATCGGGCatatat > 1:66728/1‑90 (MQ=255) ttGGGCAATTACGGCAGGGAAAAAACCTCTTTTAGCTCGCCACTCACCCGGCTGTTTTTCCGGGTAGCGGGCATAAAAACCGTCatatat > 1:410029/1‑90 (MQ=255) ttGGGCAATTACGGCAGGGAAAAAAATTTTTTCAGCTCCCTACTCACCCGGCCGTTTTCCCGGGTAGGGGGGATAAAATCGGGCatatat > 1:332761/1‑90 (MQ=255) tGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATAtt < 1:334075/90‑1 (MQ=255) gggCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATAttt > 2:2505/1‑90 (MQ=255) gggCAATTACGGCAGGGAAAAAACTTCTTTCAGCCCGCTACTCACCCGGCTGTTATCCGGGGTAGCGGGCATAAAATCGGGCATATAttt > 1:61281/1‑90 (MQ=255) ttACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAAAATCGGTCATATATTTccccca > 2:145514/1‑90 (MQ=255) ggCAGGTAAAAAACTTCTTTTAGCTCGCTACTCACCCGGCGGTTATTCGGGGTAGCGGGGATAAAAACGGTCATATTTTTTCACCCACGt > 1:542416/1‑90 (MQ=255) gTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCa < 1:413677/90‑1 (MQ=255) | ATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCA > NZ_CP009273/4086419‑4086582 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |