Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,574,833 | 0 | T | G | 55.0% | 23.9 / 27.0 | 20 | G613G (GGT→GGG) | tkt | transketolase |
Reads supporting (aligned to +/- strand): ref base T (0/9); new base G (11/0); total (11/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.95e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTA > NZ_CP009273/2574746‑2574895 | gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGAt > 2:9954/1‑90 (MQ=255) gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAt < 2:93371/90‑1 (MQ=255) gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTGcc > 2:295506/1‑90 (MQ=255) gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGTGGAAGCAGGGATTGcc > 2:172873/1‑90 (MQ=255) gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGGGGAAGCAGGGATTGcc > 2:332588/1‑90 (MQ=255) tCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACt < 1:332588/90‑1 (MQ=255) tCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAATGCCGATTACt < 2:382328/90‑1 (MQ=255) gCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCC‑TCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACa < 2:379336/90‑1 (MQ=255) cccAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGGGGGGATTGCCGATTATTGGTACa > 1:323478/1‑90 (MQ=255) gAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGtt < 1:349814/90‑1 (MQ=255) aGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTGGGACAAGTATGTTg > 1:110204/1‑90 (MQ=255) cGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa < 1:372741/90‑1 (MQ=255) cGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTGCCGATTACTGGGACAAAGATGTTGGGCTGGAAGGGGGaa > 1:355668/1‑90 (MQ=255) ttGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTTCCGGTTACTGGGACAAGTATGTTGGGCTGGAAGGGGGAAATGTc > 1:448218/1‑90 (MQ=255) tGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCg < 2:366910/90‑1 (MQ=255) cTAACGTTGCGGCTCGCGTGGCGGTGGAAGCCGGGATTGCCGATTTCTGGTACAAGGATGTTGGGTTGTAAGGGGCAATTGTCGGGGTgg > 1:443488/1‑89 (MQ=255) cTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGa > 2:278225/1‑90 (MQ=255) tAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGAc < 1:361183/90‑1 (MQ=255) ttGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGtt < 2:311055/90‑1 (MQ=255) tGCGGCTCGCGTGGCGGTGGAAGGAAGGATTGCCGGTTACTTGGAAAAGTATGTTGGGCTGAAAGGGGGAATTGTCGGGGTGACGGGGTa > 1:349977/1‑90 (MQ=255) | GACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTA > NZ_CP009273/2574746‑2574895 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |