| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 447,275 | 0 | G | C | 27.6% | 50.2 / 21.8 | 29 | A144G (GCT→GGT) ‡ | cyoE | protoheme IX farnesyltransferase |
| Reads supporting (aligned to +/- strand): ref base G (10/11); new base C (3/5); total (13/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.97e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATA > NC_000913/447134‑447412 | gccgccTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGCTGCGCCGCCgg < 4:96135/148‑1 (MQ=255) ccTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGc < 6:51117/149‑1 (MQ=255) ccTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGCTGCGCCGCCGgtgatc < 2:74726/149‑3 (MQ=255) gTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGaa < 5:120269/148‑1 (MQ=255) cTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCTCCGGCGCTGCGCCGCCGgtgatc > 6:16668/1‑137 (MQ=255) cTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCTCCGGCGCTGCGCCGCCGgtgatc < 5:16668/139‑3 (MQ=255) ttAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAAc < 8:123169/142‑1 (MQ=255) ttAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAAc > 7:123169/1‑142 (MQ=255) aGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTg < 8:32187/149‑1 (MQ=255) tGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTa > 3:36819/1‑149 (MQ=255) gCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTAACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCtata > 2:22971/1‑148 (MQ=255) ccACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGCTGCGCCGCCGgtga < 8:37965/95‑1 (MQ=255) ccACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGCTGCGCCGCCGgtga > 7:37965/1‑95 (MQ=255) acaGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGCTCTCCGGCGCTGCGCCGCCGgtgatc > 3:92954/1‑93 (MQ=37) acaGGCTGAAGAGAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGCTCTCCGGCGCTGCGCCGCCGgtgatc < 4:92954/95‑3 (MQ=25) aGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATa > 3:40061/1‑148 (MQ=255) aTAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTCACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAc < 1:22971/149‑1 (MQ=255) ccAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCg > 4:132510/1‑98 (MQ=255) ccAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCg < 3:132510/98‑1 (MQ=255) cgcTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTc < 6:41168/98‑1 (MQ=255) aCTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAgg > 4:106436/1‑98 (MQ=255) aCTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAgg < 3:106436/98‑1 (MQ=255) ccGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGccc < 5:128569/123‑1 (MQ=255) ccGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGccc > 6:128569/1‑123 (MQ=255) cGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGa > 1:142927/1‑149 (MQ=255) cGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGa < 4:36819/149‑1 (MQ=255) tAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACca < 4:40061/149‑1 (MQ=255) cGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCAcag > 5:330/1‑148 (MQ=255) cggcggCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATa > 4:92607/1‑148 (MQ=255) | GCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATA > NC_000913/447134‑447412 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |