Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,741,298 | 1 | . | T | 40.0% | 4.4 / 45.9 | 25 | intergenic (+489/+263) | gltA/yeiW | citrate synthase/putative oxidoreductase |
Reads supporting (aligned to +/- strand): ref base . (5/10); new base T (9/1); total (14/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.19e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTG > NC_002947/4741172‑4741441 | ctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGacac > 1:284426/1‑150 (MQ=255) tgCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGacac > 2:541590/1‑131 (MQ=255) tgCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGacac < 1:541590/131‑1 (MQ=255) caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt > 1:638884/1‑150 (MQ=255) caCGGAACCACGCGGCGTCACCCGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGTTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGACTTTCATTCCTGAA‑TCCAGGCACAGGCCCCGGACACGTTCTCCTGCAGCGCCTCGACCtt < 2:170468/150‑1 (MQ=255) aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg < 2:210785/150‑1 (MQ=255) cGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg < 1:126201/150‑1 (MQ=255) cACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcg < 2:284426/150‑1 (MQ=255) cGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc > 2:854269/1‑150 (MQ=255) tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg > 2:391723/1‑140 (MQ=255) tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg < 1:391723/140‑1 (MQ=255) tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑AGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc > 1:174317/1‑150 (MQ=255) ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCATCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa > 1:330160/1‑150 (MQ=255) gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATa > 2:253336/1‑150 (MQ=255) tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc > 2:479917/1‑150 (MQ=255) tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGac > 2:557734/1‑74 (MQ=255) tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGac < 1:557734/74‑1 (MQ=255) tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTgg < 1:854269/150‑1 (MQ=255) cAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAc > 1:35132/1‑150 (MQ=255) gTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGa > 2:727973/1‑141 (MQ=255) gTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGa < 1:727973/141‑1 (MQ=255) gTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAc < 1:204816/150‑1 (MQ=255) acttGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcct > 1:487598/4‑150 (MQ=255) gCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtc > 1:22840/1‑150 (MQ=255) tCGTCAA‑CCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTg < 2:330160/150‑1 (MQ=255) | CTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTG > NC_002947/4741172‑4741441 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |