Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 635,361 | G→A | 14.8% | Q6* (CAA→TAA) C403C (TGC→TGT) |
ybdM ← ybdN ← |
ParB‑like nuclease domain‑containing protein YbdM putative PAPS reductase/DUF3440 domain‑containing protein YbdN |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 635,361 | 0 | G | A | 14.8% | 54.9 / 3.9 | 27 | Q6* (CAA→TAA) C403C (TGC→TGT) | ybdM ybdN | ParB‑like nuclease domain‑containing protein YbdM putative PAPS reductase/DUF3440 domain‑containing protein YbdN |
Reads supporting (aligned to +/- strand): ref base G (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.08e-01 |
CACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCA > NC_000913/635222‑635498 | cacaGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCAt > 2:262669/1‑143 (MQ=255) acaGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATa > 1:204863/1‑143 (MQ=255) cGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGTCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATcc < 1:445971/143‑1 (MQ=255) cGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAGATTGCGTTACATCCTGCGTTACTCGTTGTTGCATAGTacgc > 2:106409/1‑140 (MQ=255) cAATCCACGGTCTGCTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATcccc < 1:208272/143‑1 (MQ=255) cACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTc < 1:390048/143‑1 (MQ=255) cACGGGCTCTTCACGAAAGGGGCGAACCTGATGGATCGCCAGGCGGATTTCATTAAGCGCTTTGATGCGATCATCATCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTc < 2:8562/143‑1 (MQ=255) ttCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTAACGttt > 1:381797/1‑143 (MQ=255) ttCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCAGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATTCCCCATTCTTGACGttt > 1:381977/1‑143 (MQ=255) cGAGAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCt < 1:437736/143‑1 (MQ=255) tAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgtt > 1:289288/1‑143 (MQ=255) aTCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGc > 1:265625/1‑143 (MQ=255) gCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCGGGCGCGCATGAGACCAAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGtt < 1:262669/143‑1 (MQ=255) tCATGACTCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGAt < 2:246265/143‑1 (MQ=255) cATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGAtt > 1:559586/1‑143 (MQ=255) cTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGGTGTTGCATAGCATCCCCCATTCCTGACGCTTCCCTTTTATCCGTTCGTTATAACGCTGACAGTTCTCCGATCTTGCgggcc > 2:146885/1‑141 (MQ=255) cGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGAt < 2:433962/143‑1 (MQ=255) tCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTCTCCGTTCGTTATAACCCCGATAGTTCTTCGATTTTGTGGGGCTAAGTGATACTCCCCGa > 2:525327/1‑143 (MQ=255) gCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTACATGATAATGCCCGACACCAGTACTCg > 2:269061/1‑143 (MQ=255) aaTTGCGTTAAGTCCTGCGTCAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTg < 2:188098/143‑1 (MQ=255) cGTTAAATCCTGCGTTAATCGTTGTTACATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGcc < 1:318382/117‑1 (MQ=255) cGTTAAATCCTGCGTTAATCGTTGTTACATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGcc > 2:318382/1‑117 (MQ=255) cGTTAAATCCTGCGTTAATCGTTGTTACATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGcc > 2:314233/1‑117 (MQ=255) cGTTAAATCCTGCGTTAATCGTTGTTACATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGcc < 1:314233/117‑1 (MQ=255) ttAAATCCTGCGTTAATCGTTGTTGCATAGTATCCACCATAACTAACGTTTGTCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTaaaa < 1:106409/143‑1 (MQ=255) agcccgCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACttt < 1:225213/138‑1 (MQ=255) aTCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATccgcc < 2:265625/143‑1 (MQ=255) ttgttgCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAACTGATAATGCCCGACACCAGTCATCGTTATTGAGTAAAACTTTGCATATCCGCCCCCa > 1:65566/1‑143 (MQ=255) | CACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCA > NC_000913/635222‑635498 |
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: ‑ |