| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_002947 | 697,097 | A→G | intergenic (+21/‑730) | mmsA‑I → / → PP_16SD | methylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 697,097 | 0 | A | G | 83.3% | 16.8 / ‑2.6 | 12 | intergenic (+21/‑730) | mmsA‑I/PP_16SD | methylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA |
| Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (10/0); total (12/0) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGAGGCCGTCAGGCC‑TCCCCGCTTTGT‑TTTTGGGGGTTGGT‑TGTTGCTGTTGGGGAGGAACTGT > NC_002947/696964‑697158 | aTGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑Ggggggg > 1:4262556/1‑136 (MQ=255) ggTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGGCGGCACGGC‑Gcccccc > 2:4600133/1‑134 (MQ=255) gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCGGGCGGGCG‑CCCCCCCGTGGT‑TTTTGGGggggggt‑tt > 2:3952088/1‑135 (MQ=255) gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCCGGC‑GGCC‑CGCCCGCGGTGT‑TTTTGGGggtgggggtg > 2:3165536/1‑132 (MQ=255) gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGGCGGCCGGGC‑GCCCACCGTTGTGTTTTTGgggtgggtgt > 1:1906351/1‑134 (MQ=255) gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGCCCGGAGGCCGCCCCCCCGTTTT‑TTTTTGGggttggt‑tt > 2:2796116/1‑135 (MQ=255) aaGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGGCGGGAGGCCGCCCCCCCGTTTT‑TTTTTGGGGtt‑gt‑tgttg > 1:1935611/1‑136 (MQ=255) aaGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGGCCGCGGAGC‑GCCCCCCCGGGT‑TTTTTGGGGGTGGt‑tgttg > 1:1446709/1‑136 (MQ=255) aaGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGCAGGCCGTCAGGCC‑TCCCCGCGGTGG‑TTTTTGGGGGTGGt‑tgttg > 1:1759436/1‑136 (MQ=255) ttGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGGCGGGAGGCC‑GCCCCCCGGTGT‑TTTTTGGGGGTGGT‑GGTTGGGGTTGggggggggctg > 1:3125679/1‑131 (MQ=255) ttGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGCCGTGAACCC‑AACCCCCGGTGG‑TTTTGGGGGGTGGTGTGTGGTTGTGGGGGGGGATCt > 2:1202996/1‑136 (MQ=255) tGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCGGGCGGGGG‑GCCCCCGGTGGT‑TTTTTGGGGTGGGT‑TGTTG‑TGTTGgggggggggggt > 1:461888/1‑130 (MQ=255) tGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGAGGCCGTCAGGCC‑TGCCCGCTTTGT‑TTTTGGGGGTTGGT‑TGTTGCTGTTGGGGAGGAACTgg > 1:345486/1‑135 (MQ=255) | ATGCCTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGAGGCCGTCAGGCC‑TCCCCGCTTTGT‑TTTTGGGGGTTGGT‑TGTTGCTGTTGGGGAGGAACTGT > NC_002947/696964‑697158 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |