| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_bme_pgi | 3,368,982 | 0 | G | C | 35.7% | 10.1 / 13.7 | 14 | intergenic (+28/+32) | yhcG/yhcH | DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily/DUF386 family protein, cupin superfamily |
| Reads supporting (aligned to +/- strand): ref base G (3/6); new base C (2/3); total (5/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCC > NC_000913_3_bme_pgi/3368843‑3369128 | tAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGa > 2:39941/1‑150 (MQ=255) aaCACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa < 1:3415/150‑1 (MQ=255) cGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCttt < 4:55039/150‑1 (MQ=255) gCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGa < 1:160618/150‑1 (MQ=255) cAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa < 1:120941/150‑1 (MQ=255) aaGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa < 1:99605/139‑1 (MQ=255) aaGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa > 2:99605/1‑139 (MQ=255) gACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa > 3:83257/1‑101 (MQ=255) gACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa < 4:83257/101‑1 (MQ=255) gAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACgg > 3:66393/1‑150 (MQ=255) gaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAg > 1:66767/2‑62 (MQ=255) gaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAg < 2:66767/61‑1 (MQ=255) tCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCCTCCTAGTTGAGGCGGTTCACCTGGCATAAACACGGCGaa < 1:39941/150‑1 (MQ=255) aCAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTccc < 4:66393/150‑1 (MQ=255) | TAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCC > NC_000913_3_bme_pgi/3368843‑3369128 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |