| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,857,143 | 1 | . | A | 50.0% | 35.6 / 39.0 | 24 | coding (51/2562 nt) | mutS | methyl‑directed mismatch repair protein |
| Reads supporting (aligned to +/- strand): ref base . (9/3); new base A (4/8); total (13/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.95e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑AT‑C‑TCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCA > NC_000913/2857034‑2857254 | tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGAT‑‑‑AGTCGCTGAAAGCCCAGCATCCCGAGAt > 1:61746/1‑139 (MQ=255)tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAAGTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGAt > 1:980735/1‑139 (MQ=255) gAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCctgct < 2:61746/139‑1 (MQ=255) tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCctg < 1:1058828/123‑1 (MQ=255) tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATgg < 1:5390/139‑1 (MQ=255) tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCctg > 2:1058828/1‑123 (MQ=255) cATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGAT‑AGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGAtttt > 1:508767/1‑139 (MQ=255) tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTa > 1:1208907/1‑139 (MQ=255) atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa < 2:206369/139‑1 (MQ=255) ccGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa > 2:17351/1‑128 (MQ=255) ccGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGAT‑AGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa < 1:17351/128‑1 (MQ=255) ccccATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGTAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCctgc < 2:893700/89‑1 (MQ=255) ccccATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCctgc > 1:893700/1‑89 (MQ=255) ccATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAAcg < 2:508767/139‑1 (MQ=255) aTGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATcc < 2:900597/74‑1 (MQ=255) aTGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATcc > 1:900597/1‑74 (MQ=255) aTGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAAcgcg < 1:253553/139‑1 (MQ=255) gTGCAATAGACAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCCTCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTc > 1:120306/1‑139 (MQ=255) gTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTc > 1:1315544/1‑139 (MQ=255) gCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGc < 2:1208907/139‑1 (MQ=255) cAATAGAAAATTTCGACGCCCATACGCCCAAGATGCAGAAGT‑GA‑T‑AGTCGCTGAAAGAAAAGAATCTCGAGATCCTGCTGTTTTACCGGAAGGGAGATTTTTATGAACTGTCTTATGACGACGCAAAACGCGCGTCGCa > 1:912089/1‑139 (MQ=255) cGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCctgct < 2:15963/66‑1 (MQ=255) cGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCctgct > 1:15963/1‑66 (MQ=255) cgcccatacgcccatGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCa > 1:104398/1‑139 (MQ=255) | TACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑AT‑C‑TCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCA > NC_000913/2857034‑2857254 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |