Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,775,223 | 0 | T | G | 60.0% | 13.0 / 21.2 | 25 | T234P (ACC→CCC) | cysE | serine O‑acetyltransferase |
Reads supporting (aligned to +/- strand): ref base T (4/6); new base G (15/0); total (19/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.19e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-02 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGTTGCAGCACCACGGAACCTGCGCCAATCTTCGCGCCGCGCCCAACTTCAATATTGCCGAGGATTTTCGCGCCC > NZ_CP009273/3775136‑3775311 | tACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGAt > 2:536926/1‑90 (MQ=255) cGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGGGGGATgcg > 1:218643/1‑90 (MQ=255) gTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATgcgg > 2:195876/1‑90 (MQ=255) gTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATgcgg > 2:107121/1‑90 (MQ=255) gTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTATCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATgcgg > 1:287075/1‑90 (MQ=255) aaTGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGGgggggcc > 1:728761/1‑89 (MQ=255) gCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCg < 1:172332/90‑1 (MQ=255) tGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGt > 1:667466/1‑90 (MQ=255) atatCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGTTgcagc > 2:49333/1‑90 (MQ=255) atatCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGCGGCGGCCCCGGGTgccgc > 1:621330/1‑90 (MQ=255) tatCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGTTgcagca < 2:667466/90‑1 (MQ=255) atCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGGGGGATGCGGGGGCCCCCGGTGCAGcac > 1:134527/1‑90 (MQ=255) cATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGTTGCAGcaccac < 1:243878/90‑1 (MQ=255) cATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGCGGCGGGACCGGGTTCAGcagcac > 2:853411/1‑90 (MQ=255) aTTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGGGGGATGCGGCGGCACCCGGTGCAGCACCACg > 1:53099/1‑90 (MQ=255) tGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGCGGCGGCACCCGGTGCAGCACCACGGAAcc > 2:189909/1‑90 (MQ=255) tATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGGGGGGTGGGGCGGCCCCGGGTGCAGCACCACGGAACCCgcgc > 1:633170/1‑90 (MQ=255) tGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGCGGCGGCACCCGGTGCAGCACCCCGGAACCCGCCCCCATCt > 1:19729/1‑90 (MQ=255) gTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGCGGCGGCACCGGTTGCAGCACCACGGAACCCGCGCCAATCtt > 2:463315/1‑90 (MQ=255) gTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGCGGCGGCACCGGGTGCAGCACCACGGAACCTGCGCCCATCtt > 2:417164/1‑90 (MQ=255) gTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGGGGGATGCGGCGGCACCGGGTGCAGCACCCCGGGAACCGCGCCCATCTTcgcgc > 1:829651/1‑90 (MQ=255) aCAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGTTGCAGCACCACGGAACCTGCGCCAATCTTCGCGCCGCGCCCa < 1:623112/90‑1 (MQ=255) gAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGtt < 1:662375/36‑1 (MQ=255) gAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGtt > 2:662375/1‑36 (MQ=255) gCCAGCGGCGGTGGTATGCGGCGGCACCGGTTGCAGCACCACGGAACCTGCGCCAATCTTCGCGCCGCGCCCAACTTCAATATTGCCGAg < 2:880496/90‑1 (MQ=255) ggtgggATGCGGCGGCACCGGGTGCAGCACCACCGAACCCGCGCCCATCTTCGCGCCGCGCCCCACTTCAATATTGCCGAGGGTTTTcgc > 1:856423/1‑90 (MQ=255) ggtgggATGCGGCGGCACCCGGTGCAGCACCACGGGACCTGCGCCCATCTTCTCGCCCCGCCCCACTTCAATATTTCCGAGGATTTTtgc > 1:880496/1‑90 (MQ=255) ggtATGCGGCGGCACCGGTTGCAGCACCACGGAACCTGCGCCAATCTTCGCGCCGCGCCCAACTTCAATATTGCCGAGGATTTTcgcgcc < 1:437415/90‑1 (MQ=255) gtATGCGGCGGGACCGGGTTCAAAACCACGGAACCTGCGCCCATCTTCGCGCCGCGCCCCACTTCAAAATTGCCCAGGATTTTcgcgccc > 1:46637/1‑90 (MQ=255) | TACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGTTGCAGCACCACGGAACCTGCGCCAATCTTCGCGCCGCGCCCAACTTCAATATTGCCGAGGATTTTCGCGCCC > NZ_CP009273/3775136‑3775311 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |