| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,741,298 | 1 | . | T | 64.3% | 60.1 / 32.5 | 28 | intergenic (+489/+263) | gltA/yeiW | citrate synthase/putative oxidoreductase |
| Reads supporting (aligned to +/- strand): ref base . (4/6); new base T (11/7); total (15/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.33e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCCGCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTG > NC_002947/4741170‑4741441 | gCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGac > 1:248868/1‑150 (MQ=255) ctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtc > 1:595650/1‑150 (MQ=255) ttGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcc > 1:551409/1‑150 (MQ=255) gATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctc < 1:808021/150‑1 (MQ=255) aaTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAgcgc > 1:871795/1‑150 (MQ=255) aTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAA‑GCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGcc < 1:829735/150‑1 (MQ=255) cagcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc < 1:52017/150‑1 (MQ=255) aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTCTAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg > 2:109712/1‑150 (MQ=255) cGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg < 2:871795/150‑1 (MQ=255) cGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg > 1:985157/1‑150 (MQ=255) aCGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcgg > 1:37611/1‑150 (MQ=255) cGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCATGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc < 2:595650/150‑1 (MQ=255) ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg < 1:76145/150‑1 (MQ=255) ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTCGCCAGCGGCGCCTTGATCACg < 1:811845/150‑1 (MQ=255) cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc < 2:248868/150‑1 (MQ=255) aTGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca < 2:957473/150‑1 (MQ=255) ggCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcaca < 2:37611/150‑1 (MQ=255) gCGTCTTGCGCTGGGGATAGGTCCAGCTCGGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac > 2:991743/1‑150 (MQ=255) aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg > 2:760188/1‑150 (MQ=255) ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTa > 1:527382/1‑150 (MQ=255) tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTgg > 2:491523/1‑150 (MQ=255) gTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCtt > 2:888677/1‑132 (MQ=255) gTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAc > 1:592576/1‑150 (MQ=255) gTCACCTTCGTGCTTGAGCAG‑CTTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCtt < 1:888677/132‑1 (MQ=255) cgacttGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtc < 2:985157/145‑1 (MQ=255) cATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGCTCGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtcc > 1:780956/1‑150 (MQ=255) ttCGTGAA‑TCCAGGCACAGGCCCCGGACACGGCCTCCTGCAGCGCCTCGACCTTTGCCAGCGACGCCTTGCTCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCtt > 2:540057/1‑150 (MQ=255) tCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTg < 2:872814/150‑1 (MQ=255) | GCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCCGCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTG > NC_002947/4741170‑4741441 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |