| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,857,148 | 0 | G | C | 57.9% | 29.6 / 24.9 | 19 | coding (56/2562 nt) | mutS | methyl‑directed mismatch repair protein |
| Reads supporting (aligned to +/- strand): ref base G (0/0); major base C (7/4); minor base . (5/3); total (12/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAA‑CCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑AT‑C‑TCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGC > NC_000913/2857034‑2857243 | tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAA‑CCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGAT‑AGTC‑‑GCTGAAAGCCCAGCATCCCGAGAt > 1:18839/1‑139 (MQ=255) aCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAA‑CCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCctg > 2:569540/1‑139 (MQ=255) cGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAA‑CCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCctgc > 1:478177/1‑139 (MQ=255) aaaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAA‑CCCCATGAGTGCAATAGAAAATTTCGACGACCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCACGAGATCCTGCTGt > 2:155446/1‑139 (MQ=255) aCCATCACACCCTAGTCAGTATCAGGAAACCGGA‑ATAAGCGCCAGGAATGCAATAGGAAGTTTCGACGCCCATAAGCCCAAGATGCAGAAGT‑GA‑T‑AATCGCTGACAGCCCAGCATCCCGAGATCCTGCTGTTTTACCgg < 2:358131/139‑1 (MQ=255) atatCAGGGAACCGGACATAA‑CCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa < 1:569540/139‑1 (MQ=255) cAGGGAACCGGACATAA‑CCCCATGAGTGCAATAGAAAATTTCGACGCCCATCCGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATACTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGt > 2:356502/1‑139 (MQ=255) gggAACCGGACATAA‑CCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGttt > 1:82493/1‑139 (MQ=255) gggAACCGGACATAA‑CCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTa < 1:165847/111‑1 (MQ=255) gggAACCGGACATAA‑CCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTa > 2:165847/1‑111 (MQ=255) gggAACCGGACATAA‑CCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGttt > 1:687977/1‑139 (MQ=255) gACATAA‑CCCCAAGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGCAGGGTGATTTTTATGAACTGTTTTATgacga > 2:364718/1‑139 (MQ=255) gCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGc < 2:478177/139‑1 (MQ=255) cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCgaga < 1:183226/72‑1 (MQ=255) cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCgaga > 2:183226/1‑72 (MQ=255) cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCctg > 2:58928/1‑77 (MQ=255) cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCctg < 1:58928/77‑1 (MQ=255) aaTAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCaa < 2:687977/139‑1 (MQ=255) gAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAActgc > 1:404883/1‑139 (MQ=255) | TACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAA‑CCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑AT‑C‑TCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGC > NC_000913/2857034‑2857243 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |