Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,741,301 | 0 | A | C | 53.3% | ‑2.2 / 36.3 | 30 | intergenic (+492/+260) | gltA/yeiW | citrate synthase/putative oxidoreductase |
Reads supporting (aligned to +/- strand): ref base A (9/5); new base C (8/8); total (17/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.84e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACC‑T‑CC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTC > NC_002947/4741166‑4741437 | aTCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑T‑TC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGcccc > 1:789791/1‑150 (MQ=255) ccGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCgg < 2:789791/150‑1 (MQ=255) ctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGtc < 1:587973/150‑1 (MQ=255) aTGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctcc > 1:363425/1‑150 (MQ=255) tgtgCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGcc > 2:531571/1‑149 (MQ=255) tgtgCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGcc < 1:531571/149‑1 (MQ=255) gtgCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc < 1:1066839/150‑1 (MQ=255) gcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑T‑TC‑TTCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCt > 2:437597/1‑150 (MQ=255) caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt > 2:869800/1‑150 (MQ=255) aCGTAACCACGCGCCGTCACCGGCCATGGCGTATTGAACGGGGTATAGGTCCTGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑TGCT‑GCGACT‑TGCATTCGTGAC‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg < 1:1020835/150‑1 (MQ=255) gTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCCGGCACATGCCCCTGGCACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTg > 1:568424/1‑150 (MQ=255) ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg < 1:294939/150‑1 (MQ=255) aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTATCGAGCGGCGCCTTGATCACGGTCACACCGCGATAg < 1:869800/150‑1 (MQ=255) ataGGTCCAGCTCTAGTCACCTTCGTGCTTGAG‑CAGC‑T‑TCTGCGACT‑TGCATTCGTGAA‑TCCAGTCACAGGCCCCGGACACGTCCTCCTGCCGCGCCTCGACCTTTGCCAGCGACGCCTTGATCACGGTCACACCGCGATAGGCCTTGTa > 2:230192/1‑150 (MQ=255) gTCCAGCTTTGGTCACCTTCATGCTAGAG‑CATC‑T‑TC‑TGCGACTTGCATTCGTGAATCCAAAC‑ACAGGCCCCGGACACGTCCTCCTCCAGCGTCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg < 2:225982/150‑1 (MQ=255) gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg > 1:493459/1‑150 (MQ=255) ctctGGTCACCTTCGTGCTTGAG‑CAGCTT‑CT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg > 1:636104/1‑150 (MQ=255) tctGGTCACCTTCGTGCTTGAGCAGCT‑T‑CT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCgg > 2:270405/1‑150 (MQ=255) tctGGTCACCTTCGTGCTTGAG‑CAGC‑T‑TCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCgg > 1:585673/1‑150 (MQ=255) cTTCGTGCTTGAG‑CAGCTT‑CT‑GCGACT‑TGCATTCGTGAA‑TNCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCAc > 1:810017/1‑150 (MQ=255) tgcgacttGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATcc > 2:599220/8‑148 (MQ=255) tgcgacttGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATcc < 1:599220/141‑1 (MQ=255) ttGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc < 2:264700/98‑1 (MQ=255) ttGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc > 1:264700/2‑99 (MQ=255) tGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgt < 2:585673/150‑1 (MQ=255) gCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtc > 2:1126522/1‑150 (MQ=255) cATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt < 1:713143/91‑1 (MQ=255) cATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt > 2:713143/1‑91 (MQ=255) ttCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc < 2:400383/61‑1 (MQ=37) ttCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc > 1:400383/1‑61 (MQ=255) | ATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACC‑T‑CC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTC > NC_002947/4741166‑4741437 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |