| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,741,305 | 0 | A | C | 50.0% | 26.8 / 46.7 | 25 | intergenic (+496/+256) | gltA/yeiW | citrate synthase/putative oxidoreductase |
| Reads supporting (aligned to +/- strand): ref base A (1/0); major base C (9/3); minor base . (7/5); total (17/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.67e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-02 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCT‑CC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCT > NC_002947/4741192‑4741439 | gcagcaCGGAATCACGCGGCGTCACCCGCCATGGCGTATTGAACCGGGTATAAGTCCAGCTGTGGTCACCTGCGTGCCTGAG‑CAGCTGCT‑GCGACT‑TGCATGCGTGAA‑TCCCGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGa < 1:1377115/150‑1 (MQ=255) caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc < 2:1011798/147‑1 (MQ=255) caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt > 1:1027167/1‑150 (MQ=255) caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt > 2:955627/1‑150 (MQ=255) caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc > 1:1011798/1‑147 (MQ=255) cGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg > 1:1049718/1‑150 (MQ=255) aaTCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CACAG‑GCTGCGACT‑TGAATTAGAAAC‑TGCCAAAACAGGTTCCGTAGACGTCCTCCTGCAGCGCCTCGACCTTGGCCa > 1:752183/1‑150 (MQ=255) aTCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑TCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAg > 2:76468/1‑150 (MQ=255) cggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc < 1:955627/150‑1 (MQ=255) aCCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CCGCTGCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACGTTGGCCAGCGGCGCCTTGATc > 2:1044140/1‑150 (MQ=255) ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa > 1:243903/1‑150 (MQ=255) gTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGcc > 2:720676/1‑123 (MQ=255) gTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑TC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGcc < 1:720676/123‑1 (MQ=255) gTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTg > 1:1026165/1‑150 (MQ=255) gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGcc < 2:114837/93‑1 (MQ=255) gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGcc > 1:114837/1‑93 (MQ=255) ctctGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg < 2:1027167/150‑1 (MQ=255) ctctGGTCACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg < 1:76468/150‑1 (MQ=255) gTCACCTTCGTGCTTGAG‑CAGCT‑TC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAc > 2:778222/1‑150 (MQ=255) cACCTTCGTGCTTGAGCAGCTT‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac > 1:882617/1‑117 (MQ=255) cACCTGCGTGCTTGAG‑CAGCT‑TC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac < 2:882617/117‑1 (MQ=255) acgttCGTGCTTGAG‑CAGCT‑TC‑TGCGACTTGCATTCGTGAA‑TCCGGGCACAGACCCCGGACACGTCCTGCTGCAACGCCGCGACCTTTGCCAGCAGCGCCTTGATCACGGTCACCCCGCGCTACACCTTGTACTTGGACCTCGCGACCTc > 2:409388/4‑150 (MQ=255) gCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtc > 2:1392375/1‑150 (MQ=255) gCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtc > 1:1283012/1‑150 (MQ=255) aTTCGTGAA‑TCCAGGCACAGGCCCTGGAAACGTCCTCCTGCAACGGCTCGACCTTGGCCAGCAGCGCCTTTATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGTCCACCTAACTGAGCGACACCTTGATCGcgtccgcgtcct > 1:611010/1‑150 (MQ=255) | GCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCT‑CC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCT > NC_002947/4741192‑4741439 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |