Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,472,583 | 0 | A | G | 57.1% | 7.3 / 16.3 | 14 | E34G (GAG→GGG) | ynbB | phosphatidate cytidylyltransferase |
Reads supporting (aligned to +/- strand): ref base A (0/6); new base G (8/0); total (8/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.33e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGC‑GG‑TTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCAC > NZ_CP009273/1472494‑1472651 | aTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGC‑GG‑TTA‑CCGGAGAGAAAAGGGGGTGa < 1:333995/90‑1 (MQ=255) cAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGC‑GG‑TTA‑CCGGAGAGAAAAGGGGGTGAGGTTACAt < 2:654908/90‑1 (MQ=255) acTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGC‑GG‑TTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACg < 2:57105/90‑1 (MQ=255) gttgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGC‑GG‑TTATCGGGGGGGAAAAGGGGGGGGGGTT‑CATTTCGTATTTGttc > 2:601276/1‑88 (MQ=255) tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGC‑GGTTTA‑CCGGGGGGAAAAGGGGGGGGGGGT‑CATTACTTATTTTTACCt > 2:490084/1‑90 (MQ=255) tAGCAACGCTGATAAACCGTTTTCTGCTATGGC‑GGTTTT‑CGGGGGGGAAAAGGGGGGGGGGTT‑CATTTTGTATTTTTTCCTTGTggggga > 2:486153/1‑88 (MQ=255) tAGCAACGCTGATAAACCGTTTTCTGCTATGGC‑GG‑TTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCa < 2:502211/90‑1 (MQ=255) tAGCAACGCTGATAAACCGTTTTCTGCTATGGC‑GG‑TTA‑CCGGAGAGAAAAGGGGGGGGGGGTACATTAAGTATTCTTACCTGggggggga > 1:276561/1‑88 (MQ=255) gATAAACCGTTTTCTGCTATGGC‑GG‑TTA‑CCGGGGAGAAAAAGGGGGGGGGGTTCATTAAGGATTTTTACCTGGTGGGGGATCGGCCTTTg > 2:46430/1‑90 (MQ=255) gATAAACCGTTTTCTGCTATGGC‑GG‑TTA‑CCGGAGAGAAAAGGGGGGGGGGGTATATTTCGTATTTGTATCTGCGGGGGGGTGGTCATTTg > 2:426994/1‑90 (MQ=255) gTTTTCTGCTATGGCGGG‑TTA‑CCGG‑GGGAAAAAGGGGGGGGGGTACATTTCGTTTTTGTACCTGGGGGGGGATCGTCGTtttttttttta > 2:331519/1‑87 (MQ=255) gTTTTCTGCTATGGC‑GG‑TTA‑CCGGGGGGAAAAGGGGGGGGGGGTACATTTCGTATTTTTACCCTGGGGGGGATCGTCATtttttttttta > 2:212631/1‑87 (MQ=255) cTGCTATGGC‑GG‑TTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTg < 1:586631/90‑1 (MQ=255) cTGCTATGGC‑GG‑TTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTg < 2:582634/90‑1 (MQ=255) tA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCAc < 1:684749/90‑1 (MQ=255) | ATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGC‑GG‑TTA‑CCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCAC > NZ_CP009273/1472494‑1472651 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |