| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,741,305 | 0 | A | . | 66.7% | 43.1 / 44.0 | 33 | intergenic (+496/+256) | gltA/yeiW | citrate synthase/putative oxidoreductase |
| Reads supporting (aligned to +/- strand): ref base A (0/0); major base . (13/9); minor base C (5/6); total (18/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.88e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCC‑GCCAAC‑T‑GCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTT > NC_002947/4741185‑4741440 | gTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAgc > 1:1197399/1‑150 (MQ=255) aTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTT‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGcc < 2:1208866/150‑1 (MQ=255) acagAATCACGCGGCGTCACCGGCCATGGCGTATTGGACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg < 2:527236/147‑1 (MQ=255) aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg > 2:349433/1‑150 (MQ=255) aTCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctc < 2:1110846/126‑1 (MQ=255) aTCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctc > 1:1110846/1‑126 (MQ=255) tGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac < 1:1196120/150‑1 (MQ=255) gCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctcc < 2:540890/104‑1 (MQ=255) gCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctcc > 1:540890/1‑104 (MQ=255) aTTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCa < 1:904427/104‑1 (MQ=255) aTTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCa > 2:904427/1‑104 (MQ=255) gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATa < 2:472567/150‑1 (MQ=255) aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg > 1:629786/1‑150 (MQ=255) aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg > 1:372600/1‑150 (MQ=255) cTGGGTATAGGCCCCGCCATCGACCCCTTCGGGCTTGAG‑CA‑GGTGT‑TGCGAC‑TGGCATTCGTGAA‑TCCAGGCACAGGCCCCGGCCACGTCCTCCTCCAGCGCGGTGATCTCGGCTTGGTTCTCCTTGATCACGGTCACACCGCGATAGGc > 2:896096/1‑150 (MQ=255) tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg < 2:360361/128‑1 (MQ=255) tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg > 1:360361/1‑128 (MQ=255) ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGaca > 2:1455286/1‑80 (MQ=255) ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTa < 2:1668107/150‑1 (MQ=255) ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGaca < 1:1455286/80‑1 (MQ=255) taGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑T‑GCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTAc < 1:841227/150‑1 (MQ=255) ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt < 1:1540769/150‑1 (MQ=255) ccAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTTAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGa > 2:1559147/1‑150 (MQ=255) ctctGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg > 1:1750156/1‑150 (MQ=255) tCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTg < 2:1616792/74‑1 (MQ=255) tCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTg > 1:1616792/1‑74 (MQ=255) tCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑T‑ACATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAcc > 1:812914/1‑150 (MQ=255) gactt‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcc > 1:573084/5‑150 (MQ=255) actt‑GCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcct > 1:887332/4‑150 (MQ=255) gCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtc < 2:372600/150‑1 (MQ=255) aTTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc > 2:333837/1‑87 (MQ=255) aTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc < 1:333837/87‑1 (MQ=255) ttCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCtt > 2:56156/1‑150 (MQ=255) | GTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCC‑GCCAAC‑T‑GCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTT > NC_002947/4741185‑4741440 |
| 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: ‑ |