Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,086,499 | 0 | C | A | 56.0% | 7.6 / 33.6 | 25 | V100V (GTG→GTT) | rhaM | L‑rhamnose mutarotase |
Reads supporting (aligned to +/- strand): ref base C (3/8); new base A (14/0); total (17/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.53e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGC > NZ_CP009273/4086410‑4086581 | agaTGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAac < 2:366113/90‑1 (MQ=255) gTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAg > 2:463435/1‑90 (MQ=255) tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGGAAAAACCTTTTTCAGc > 2:444358/1‑90 (MQ=255) tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGc > 2:265862/1‑90 (MQ=255) tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTGCAGc > 1:161714/1‑90 (MQ=255) tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTGCAGc > 2:274485/1‑90 (MQ=255) tAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGc > 2:137950/1‑85 (MQ=255) tAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGc < 1:137950/85‑1 (MQ=255) aTTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGc < 1:299350/90‑1 (MQ=255) aaaaTATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATc < 2:538861/90‑1 (MQ=255) aaaTATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATcc < 1:335934/90‑1 (MQ=255) tgtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCCCCACTCACCCGGCTGTTTTTCGGGGTa > 1:363633/1‑90 (MQ=255) gtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAg < 2:800646/90‑1 (MQ=255) gtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCGCTACTCACCCGGCTGTTTTTCGGGGTAg > 1:381334/1‑90 (MQ=255) gCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCg < 2:161714/90‑1 (MQ=255) tctcTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGGATAACAt > 2:47202/1‑90 (MQ=255) ctctGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATc < 1:579487/90‑1 (MQ=255) gCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCGCTACTCACCCGGCTGTTATCCGGGGTAGCGGGCATAAAATCGGTCa > 2:449386/1‑90 (MQ=255) acaTTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCCCTACTCACCCGGCTGTTTTCCGGGGTAGCGGGCATAAAATCCGGCata > 1:366778/1‑90 (MQ=255) aTTGGGCAATTACGGCAGGTAAAAAACTTCTTTCAGCTCGCCACTCCCCCGGCTGTTATCCCGGGTAGGGGGCATAAAATCGGGCatata > 1:165414/1‑90 (MQ=255) aTTGGGCAATTACGGCAGGTAAAAAACTTCTTGCAGCTCGCTACTCACCCGGCCGGTATCCGGGTTAGCGGGCAAAAAATCGGTCatata > 2:579487/1‑90 (MQ=255) ttGGGCAATTACGGCAGGTAAAAAACTTCTTGCAGCTCGCTACTCACCCGGCTGTTATCCGGGGTAGCGGGCATAAAATCGGTCatatat > 2:173159/1‑90 (MQ=255) gggCAATTACGGCAGGGAAAACACTTCTTTCAGCTCGCTACTCACCCGGCTGTTTTCCGGGTTAGCGGGGAAAAACTCGGGCATATAttt > 1:156916/1‑90 (MQ=255) gggCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCGCTACTCACCCGGCTGTTTTTCCGGGTAGGGGGCATAACATCGGTCCTATAttt > 1:734338/1‑90 (MQ=255) gggCAATTACGGCAGGGAAAAAACTTCTTGCAGCTCGCTACTCACCCGGCTGTTATCCGGGGTAGCGGGCATAAAATCGGGCATATAttt > 2:596048/1‑90 (MQ=255) aTTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTccacc < 2:700824/90‑1 (MQ=255) gggaaaaaaCTTCTTTCAGCTCGCTACTCCCCCGGCTGTTTTCCGGGGTAGCGGGCATAACATCCGGCATATATTTTCCCCCACGTTGGc > 1:769780/4‑90 (MQ=255) | AGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGC > NZ_CP009273/4086410‑4086581 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |