Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 234,876 | T→C | 10.7% | W21R (TGG→CGG) | yafS → | putative S‑adenosyl‑L‑methionine‑dependent methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 234,876 | 0 | T | C | 10.7% | 81.3 / 4.0 | 38 | W21R (TGG→CGG) | yafS | putative S‑adenosyl‑L‑methionine‑dependent methyltransferase |
Reads supporting (aligned to +/- strand): ref base T (10/23); new base C (2/2); total (12/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.82e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01 |
TCAAAACCCAGATGTAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTG > NC_000913/234731‑235023 | tCAAAACCCAGATGTAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTggg < 1:745696/149‑1 (MQ=255) tCAAAACCCAGATGGAAATGACAGCAAAGGCGGTAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTggg < 1:213116/149‑1 (MQ=255) aaCCCAGATGTAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGaaa < 2:236730/149‑1 (MQ=255) ccAGATGTAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCt > 2:774864/1‑149 (MQ=255) cGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGGCCCTCCGACTGACGTGGATAACGCTTGCAGAAGCGATTTCCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGGGCGcc > 1:149835/1‑149 (MQ=255) cGGGAATACTGTTAAGATTCATAGATTAAACACGAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTCTATCGCAAGGCGCGGGATCGcc > 2:363872/1‑149 (MQ=255) gAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGc < 1:774864/149‑1 (MQ=255) aCTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCaa < 1:110288/149‑1 (MQ=255) gTTAAGTTTCATAGATTACCTCTCGGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCGAACTGTCGAGGCTCCTGATTGGTGGGGCGATTTGCCCTGGGGAAAGCTTTATGGCAAGGCGCTGGAGCGCCAGCTCAAccc < 2:226235/149‑1 (MQ=255) gaGATTCAGATGTTACCTCTCAGGGGGGAACGGAAGGTTGTGATGAAACCGGCAAGAGGCCCTCAAACGGTCGTGGCTCCTGATTGCTGGGGCGTTGGTCGCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTg < 1:347392/148‑1 (MQ=255) aaGATTTAGGGATTACCTCTCACTGTGACACGGAAGGTTGTGACGAAACCGGCAAGCGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCTAGGCGCTGGAGCGCCAGCTCAACCCGTg < 1:67506/149‑1 (MQ=255) aGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTgg < 2:77086/149‑1 (MQ=255) tCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCa < 2:725760/149‑1 (MQ=255) tCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCCGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCa > 2:403489/1‑149 (MQ=255) tCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCCGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCa > 2:305305/1‑149 (MQ=255) gATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaa < 2:97739/149‑1 (MQ=255) ctcAGTGTGAAACGGCAGGTTGTGATGAAACCGGCAAGAGTCCCTAAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGGGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATgg > 2:191156/1‑149 (MQ=255) ctcAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATgg < 2:79677/149‑1 (MQ=255) cgcAGTGTGAAACGGAAGGTTGTGATGAAACGGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTGGCCCGGGGGAATGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATgg < 2:710985/147‑1 (MQ=255) gtgtGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGtttt < 2:663948/149‑1 (MQ=255) tgAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGCGCGCCAGCTCAACCCGTGGTTCCCTACAATGTATGGCTTTCAt > 1:547510/1‑149 (MQ=255) gAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGCGTTTCATc > 2:608836/1‑149 (MQ=255) cGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCt < 1:633801/149‑1 (MQ=255) cGGAACGTTGTCATGAAAGCCGCAAGCGCCCCTCGAACTGTCGTGGCTCCTGATTGCTGGGGCGATTCGCCCTGGGGAAAGCTTTATCTCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCt < 1:86485/149‑1 (MQ=255) aaGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTTTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTaa < 2:252558/149‑1 (MQ=255) aaGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAGGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTaa < 1:307188/149‑1 (MQ=255) ggatgtgATGAAACCGGCAAGAGTCCCTCAACCTGTCGTCGCTCCTGATTGCAGGGGCGATTTGCCCTGGGGAAAGCTGTATCGCAAGGCGCTGGTGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGa < 2:552210/146‑1 (MQ=255) gtgATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTgg < 1:97259/149‑1 (MQ=255) aaCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCGGGAGCGCCAGCTCAACCCGTGGTTCACTAACATGGATGGTTCTCATCTGCTTAAGCTTGGCACTTTa > 2:11927/1‑149 (MQ=255) cGGCGAGAGTCCCTCAAGCTGTCGTGTCTCCTGATTGCTGGGGCGGTTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgc < 2:254230/149‑1 (MQ=255) tCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAACTGTATGGTTTTCATCTGCTTAAGATTGGCGATTTAAGCGCAGAAATCAAtt > 2:586021/1‑149 (MQ=255) tCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAAtt > 2:698806/1‑149 (MQ=255) cTCCTGATTGCTGGGGCGATTTGCCCCGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGcgcg < 1:403489/149‑1 (MQ=255) cTCCTGATTGCTGGGGCGATTTGCCCCGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGcgcg < 1:305305/149‑1 (MQ=255) gggCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAg > 1:282029/1‑149 (MQ=255) gCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGttt < 2:190785/149‑1 (MQ=255) cTGGGGAAAGCTTTATCGGAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTGAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTg < 1:618497/149‑1 (MQ=255) cTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTg < 1:431909/149‑1 (MQ=255) | TCAAAACCCAGATGTAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTG > NC_000913/234731‑235023 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |