Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,927,661 | 0 | A | T | 50.0% | 23.4 / 21.9 | 32 | S401T (TCA→ACA) | edd | phosphogluconate dehydratase |
Reads supporting (aligned to +/- strand): ref base A (14/2); new base T (0/16); total (14/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.09e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGC > NZ_CP009273/1927573‑1927742 | gCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGAt > 2:61480/1‑90 (MQ=255) gCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGAt > 1:177986/1‑90 (MQ=255) gCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTc < 1:228168/90‑1 (MQ=255) ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc > 2:384060/1‑90 (MQ=255) ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc > 1:258357/1‑90 (MQ=255) ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc > 2:22582/1‑90 (MQ=255) ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc > 1:347402/1‑90 (MQ=255) ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc > 1:373257/1‑90 (MQ=255) ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc > 2:127344/1‑90 (MQ=255) ttAACTCTTTTGTCCCACAATGATGAGAAAAAGGTTTTTCGAAGGAAGCGATCACATTGCTGTCGAGTGTTTTTCCCGCCCCTTCCCGcc < 2:331786/90‑1 (MQ=255) ctttttGTCCCCCCATGATGAGAAAAAGGTTTTTCGAAGGAAGCGATCACTTTGCTGTCGAGTTTTTTTTCCCCCCCTTCCCGccagtcc < 2:258357/88‑1 (MQ=255) ctttttGCCCCCCCATGATGAGAAAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGTTTTTCCCCCCCCTTCCCGccagtcc < 2:373257/88‑1 (MQ=255) ccccccATGAGGAGAGAAGGGTTGTTCGAAGGAAGCGATCACTTTGCTGTCGAGTTTTTTTCCCGCCCCTCCCCGCCAGTCCAGTTCAcc < 1:130519/90‑1 (MQ=255) caTGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCattat > 1:346147/1‑90 (MQ=255) gatgaGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCa > 2:208716/1‑90 (MQ=255) gagagaAAGTTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGcc < 1:13486/90‑1 (MQ=255) gaaaaaGTTTGTTCGAAGGAAGCGACCACTTTGCTGTGGAGTTTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt < 1:22582/87‑1 (MQ=255) gaaaaaGGTTTTTCGAAGGAAGCGATCACATTGCTGTCGAGTGTTTTTTCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt < 2:177986/87‑1 (MQ=255) gaaaaaGGTTGTCGAAAGGAAGCGACCACATTGCTGTGGAGTGTTTTTTCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt < 2:235596/87‑1 (MQ=255) aGGTTGTTCGAAGGAAGCGTTCACATTGCTGTGGAGTTATTTTTCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTc < 2:347402/90‑1 (MQ=255) aGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTTTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTc < 1:167828/90‑1 (MQ=255) gtttttCGAAGGAGCCGATCACTTTGCTGTGGAGTTTTTTTCCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCaa < 1:145091/90‑1 (MQ=255) gttgttCGAAGGAAGCGACCACATTGCTGTCGAGTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCaa < 1:57255/90‑1 (MQ=255) gttCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAggg > 2:407149/1‑90 (MQ=255) aggaagCGATCACTTTGCTGTCGAGTTTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAAc < 1:251445/90‑1 (MQ=255) gCGATCACTTTGCTGTCGAGTGTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac < 1:73135/90‑1 (MQ=255) gCGATCACTTTGCTGGGGAGTTTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac < 2:352244/90‑1 (MQ=255) gCGACCACTTTGCTGTCGAGTTTTTTTTCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac < 1:4264/90‑1 (MQ=255) caTTGCTGTCGAGTTTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCaa < 1:61480/90‑1 (MQ=255) tGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTg > 1:26335/1‑90 (MQ=255) gTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGc > 1:88821/1‑90 (MQ=255) gTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGc > 1:13743/1‑90 (MQ=255) | GCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGC > NZ_CP009273/1927573‑1927742 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |