| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,442,007 | T→C | E57G (GAG→GGG) | msrA ← | methionine sulfoxide reductase A |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,442,007 | 0 | T | C | 96.8% | 102.3 / ‑4.4 | 31 | E57G (GAG→GGG) | msrA | methionine sulfoxide reductase A |
| Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (12/18); total (13/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.19e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01 | |||||||||||
TCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCG > NC_000913/4441872‑4442141 | tCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCca < 2:753191/139‑1 (MQ=255) aTGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGa > 2:485238/1‑139 (MQ=255) tGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGaa > 2:915079/1‑139 (MQ=255) gCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGaaa > 2:792438/1‑139 (MQ=255) cccGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCcac > 2:636019/1‑110 (MQ=255) cccGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCcac < 1:636019/110‑1 (MQ=255) gTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCaa < 1:474664/139‑1 (MQ=255) ccGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCa < 2:582230/139‑1 (MQ=255) ccGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCa < 1:953885/139‑1 (MQ=255) acTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAg > 1:927716/1‑139 (MQ=255) acTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAg < 1:991336/139‑1 (MQ=255) cTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAgg < 1:915079/139‑1 (MQ=255) ttCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGt > 2:459656/1‑139 (MQ=255) cccGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAGTGGCAATCTCCATTCCGTCAGGTAc < 1:191592/139‑1 (MQ=255) cccGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTc > 1:766678/1‑133 (MQ=255) cccGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTc < 2:766678/133‑1 (MQ=255) cGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACAt < 2:1356/139‑1 (MQ=255) taGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGa > 2:916564/1‑139 (MQ=255) cgGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATg < 1:485238/139‑1 (MQ=255) ggTATAGCCTGTGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc < 2:462204/87‑1 (MQ=255) ggTATAGCCTGTGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGc > 1:462204/1‑87 (MQ=255) ataGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGc < 1:153402/139‑1 (MQ=255) taGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCg < 1:916564/139‑1 (MQ=255) ggcggTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCCCCACACCCCAGGAACAACCCATCGCAAAAATGGCAATCTCCAATCCGTCAGGTAc < 2:37246/96‑1 (MQ=255) ggcggTGCTGTAAACGCCGGGTAACCGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCAATCCGTCAGGTAc > 1:37246/1‑96 (MQ=255) tgggTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTcc > 1:789978/2‑65 (MQ=255) cGGGTAACTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTcc < 2:789978/65‑1 (MQ=255) cTGCCAGAACAGACGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCcaggc < 1:792438/139‑1 (MQ=255) aGAACAGACGCCCCACACCCCAGAAACAACCCATCGC‑AAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTa < 2:562638/107‑1 (MQ=255) aGAACAGACGCCCCACACCCCAGAAACAACCCATCGC‑AAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTa > 1:562638/1‑107 (MQ=255) aCGCCCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATcggcg > 1:55268/1‑139 (MQ=255) | TCGTAAACAATGCGTACCGCTTCGGCATGACCCGTATCACCGGAGCACACTTCCCGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCG > NC_000913/4441872‑4442141 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |