| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,501,014 | C→A | 14.8% | G31V (GGG→GTG) | PP_3989 ← | DNA‑cytosine methyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,501,014 | 0 | C | A | 14.8% | 67.5 / 3.8 | 27 | G31V (GGG→GTG) | PP_3989 | DNA‑cytosine methyltransferase |
| Reads supporting (aligned to +/- strand): ref base C (11/12); new base A (2/2); total (13/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01 | |||||||||||
CTCAGGCGACGTCTCTAGTAGATCGACCTTCGAAGCCCTTTGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAGGGTTACTCTCGGAGGGTGCCCATTGCAAGCACTCCGATCCATTTTTACTCATGATGGTTTTTGCCTGTAAAAAAATGCCACAGTTTTACCATAGGCGGTCATCTGA > NC_002947/4500866‑4501159 | cTCAGGCGACGTCTCTAGTAGATCGACCTTCGAAGCCCTTTGTCCTGGATGGTTATGATTATACGTCTCCACAACATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGcccc > 1:159452/1‑150 (MQ=255) tCAGGCGACGTCTCTAGTAGATCGACCTTCGACGCCCTTGGTCCTGGGTGGTGATGATTATACGTCTCCTCAACATCACGACCGTCAACAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCCGATCCCATCCACCATCCCca > 1:403974/1‑150 (MQ=255) gTCTCTAGTAGATCGACCTTCGAAGCCCTTTGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAgg > 2:118820/1‑150 (MQ=255) gTCTCTAGTAGAGCGACCTACGAAGCCCTTTGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAACCCCGCATCAATGACCCCCCTATCCAACCCACCAGCCCCACAAAACTgg < 2:276067/150‑1 (MQ=255) tctAGTAGATCGACCTTCGAAGCCCTTTGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTa > 1:23924/1‑150 (MQ=255) tctAGTAGATCGACCTTCGAAGCCCTTTGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCACACAAAACAGGCTa > 1:81369/1‑150 (MQ=255) tctAGTAGATCGACCTTCGAAGCCCTTTGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCACACAAAACAGGCTa > 1:81320/1‑150 (MQ=255) ctAGTAGATCGACCTTCGAAGCCCTTTGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTaa < 1:118820/150‑1 (MQ=255) gACCTTCGAAGCCCTTTGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTgg > 2:305366/1‑150 (MQ=255) ttCGAAGCCCTTTGTCCTGGATGGTTATGATTATACGTCCCCACACCATACCGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGc > 1:209366/1‑150 (MQ=255) aGCCCTTTGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAg < 2:23924/150‑1 (MQ=255) cTTTGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAg < 1:107687/122‑1 (MQ=255) cTTTGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAg > 2:107687/1‑122 (MQ=255) cTTTGTCCTGGATGGTTATGATTATACGTCTCCACAACATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCa < 2:403974/150‑1 (MQ=255) tGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCACACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAggg < 2:81320/150‑1 (MQ=255) tGTCCTGGATGGTTATGATTATAAGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCACACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAggg < 2:81369/150‑1 (MQ=255) tGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAGGGTTAct > 2:213053/1‑150 (MQ=255) tcCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAGGGTTACTCTCGGAGGGTGCCCATTGCAAg > 2:392131/1‑150 (MQ=255) tCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAGGGTTACTCTCGGAGGGTGCCCATTGCAAGCTCTCCGATCCAttt < 2:159452/150‑1 (MQ=255) gAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTGGt > 2:336062/1‑67 (MQ=255) gAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTGGt < 1:336062/67‑1 (MQ=255) cATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCACCTTTAGTTGCTCCACGTCAGGGTTACTCTCGGAGAGTGCCCGTTGCAAGCACTCCTATCCATTTTTACTCATGATGGTTTTTGCCTGTAAAAAAATGc < 1:211782/150‑1 (MQ=255) tCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAGGGTTACTCTCGGAGGGTGCCCATTGCAAGCACTCCGATCCATTTTTACTCATGATGGTTTTTGCCTGTAAAAAAATGCca < 1:22012/150‑1 (MQ=255) aTGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAGGGTTACTCTCGGAGGGTGCCCATTGCAAGCACTCCGATCCATTTTTACTCATGATGGTTTTTGCCTGTAAAAAAATGCCACAg < 1:213053/150‑1 (MQ=255) ccaatccaCCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAGGGTTACTCTCGGAGGGTGCCCATTGCAAGCACTCCGATCCATTTTTACTCATGATGGTTTTTGCCTGTAAAAAAATGCCACAGTTTTACCATAGGc > 1:377870/1‑150 (MQ=255) aatccaCCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAGGGTTActc > 2:342290/1‑60 (MQ=255) aatccaCCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAGGGTTActc < 1:342290/60‑1 (MQ=255) aGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAGGGTTACTCTCGGAGGGTGCCCATTGCAAGCACTCCGATCCATTTTTACTCATGATGGTTTTTGCCTGTAAAAAAATGCCACAGTTTTACCATAGGCGGTCATCTGa < 1:305366/150‑1 (MQ=255) | CTCAGGCGACGTCTCTAGTAGATCGACCTTCGAAGCCCTTTGTCCTGGATGGTTATGATTATACGTCTCCACAGCATAACGATCGTGATCAGCAGCGAAAACAACATCGAAACCGGCATCAATGAACCCCAGATCCAATCCACCAGCCCCACAAAACAGGCTAACCAGCTTTGGTTGCTCCAGGTCAGGGTTACTCTCGGAGGGTGCCCATTGCAAGCACTCCGATCCATTTTTACTCATGATGGTTTTTGCCTGTAAAAAAATGCCACAGTTTTACCATAGGCGGTCATCTGA > NC_002947/4500866‑4501159 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |