| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,226,550 | 0 | C | A | 65.0% | 8.8 / 20.3 | 20 | A908E (GCG→GAG) | metH | homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent |
| Reads supporting (aligned to +/- strand): ref base C (3/4); new base A (7/6); total (10/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCGCGTACCGTTGGTGTGGTGGCGGCGCTGCTTTCCGATACCCAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGC > NC_000913/4226408‑4226693 | tCGCGTACCGTTGGTGTGGTGGCGGCGCTGCTTTCCGATACCCAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGAGgcgcg > 4:148735/1‑149 (MQ=255) tggCGGCGCTGCTTTCCGATACCCAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCtttt > 2:177745/1‑149 (MQ=255) ccGATACCCAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGAGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTa < 3:148735/148‑1 (MQ=255) cGATACCCAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTaca < 2:99557/149‑1 (MQ=255) cGATACCCAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTaca < 1:177745/149‑1 (MQ=255) gAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGAGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAgg < 6:123819/136‑1 (MQ=255) gAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGAGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAgg > 5:123819/1‑136 (MQ=255) aTTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGAGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCt > 8:128741/1‑149 (MQ=255) ttCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGAGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTg > 8:81445/1‑149 (MQ=255) cACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTaa < 4:162991/149‑1 (MQ=255) aaGAAACCGCGCACACCACCGGTCACGCTGGAAGAGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATc < 4:18413/147‑1 (MQ=255) aaGAAACCGCGCACACCACCGGTCACGCTGGAAGAGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATc > 3:18413/1‑147 (MQ=255) aaGAAACCGCGCACACCACCGGTCACGCTGGAAGAGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGa < 7:81445/149‑1 (MQ=255) aaGAAACCGCGCACACCACCGGTCACGCTGGAAGAGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGa > 3:150916/1‑149 (MQ=255) gAAACCGCGCACACCACCGGTCACGCTGGAAGAGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACt > 5:117099/1‑149 (MQ=255) aaCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTg > 7:119990/1‑148 (MQ=255) aCCGCGCACACCACCGGTCACGCTGGAAGAGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGa < 4:150916/149‑1 (MQ=255) ccaccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCTCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACg < 8:59648/117‑1 (MQ=255) ccaccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACg > 7:59648/1‑117 (MQ=255) gAAGAGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGc < 6:117099/148‑1 (MQ=255) | TCGCGTACCGTTGGTGTGGTGGCGGCGCTGCTTTCCGATACCCAGCGTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGC > NC_000913/4226408‑4226693 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |