Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,741,301 | 0 | A | C | 63.3% | 24.2 / 24.8 | 30 | intergenic (+492/+260) | gltA/yeiW | citrate synthase/putative oxidoreductase |
Reads supporting (aligned to +/- strand): ref base A (5/6); new base C (7/12); total (12/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.12e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.30e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCC‑GCCAAC‑T‑GCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTT > NC_002947/4741174‑4741440 | cctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTACAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGt < 2:753443/150‑1 (MQ=255) ctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGtc < 1:37024/150‑1 (MQ=255) tGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTcctcct > 2:764361/1‑150 (MQ=255) aTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACC‑T‑GCATTCGTGAA‑CCCAGGCACACGCCCCGGACACGTCCTCCTGCAGCGcc < 2:426932/150‑1 (MQ=255) gtgCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGac < 2:97836/130‑1 (MQ=255) gtgCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTT‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGac > 1:97836/1‑130 (MQ=255) cagcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTT‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAc < 1:533042/150‑1 (MQ=255) aCGGAATCCCGCGGCGTCACCGGCAATGGCGTATTGAACATGGTATAGGGCCAGCTCTGGTTACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg < 2:258083/150‑1 (MQ=255) aCGGAAGCACGCGGCGTCACCGCCCATGGCGTATTGAACTCGGTATACGTCCAGCTCTGGTCACCTTCGTGCTTGCG‑CCGCTTCT‑GCCACT‑T‑GCATTCGTGAATCCAGGC‑ACAGCCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg < 1:427322/150‑1 (MQ=255) cggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGATCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc < 1:1124291/150‑1 (MQ=255) cACCGGCCATGGCGTATGGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt < 2:639702/150‑1 (MQ=255) aCCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTT‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc > 2:650717/1‑150 (MQ=255) ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCAGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg < 2:1204341/150‑1 (MQ=255) ccATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt < 1:764361/150‑1 (MQ=255) ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa > 2:379665/1‑150 (MQ=255) ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTT‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa > 2:547876/1‑150 (MQ=255) tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc < 2:605001/150‑1 (MQ=255) ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TAGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCtt < 1:482412/150‑1 (MQ=255) gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCACCTCGAcc > 1:481109/1‑99 (MQ=255) gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCACCTCGAcc < 2:481109/99‑1 (MQ=255) tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCTAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTgg > 1:703982/1‑150 (MQ=255) tCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGACACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTgg > 2:492817/1‑150 (MQ=255) gCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg < 2:1034092/118‑1 (MQ=255) gCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg > 1:1034092/1‑118 (MQ=255) ctctGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg < 1:1253778/150‑1 (MQ=255) tt‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcg > 2:257014/2‑150 (MQ=255) ttgtaTTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcg > 2:1423788/5‑150 (MQ=255) gCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtc < 2:1120166/150‑1 (MQ=255) cATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtcc < 2:345468/150‑1 (MQ=255) ttCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCtt > 1:58922/1‑150 (MQ=255) | CCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCC‑GCCAAC‑T‑GCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTT > NC_002947/4741174‑4741440 |
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: ‑ |