Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 578,130 | A→C | 40.2% | I8I (ATA→ATC) | rzpD → | DLP12 prophage; putative murein endopeptidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 578,130 | 0 | A | C | 40.2% | 67.5 / 63.7 | 65 | I8I (ATA→ATC) | rzpD | DLP12 prophage; putative murein endopeptidase |
Reads supporting (aligned to +/- strand): ref base A (21/18); new base C (13/13); total (34/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.04e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01 |
AGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATG > NC_000913/577986‑578268 | aGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGAGGACTCGTCGTGAGATTGCGCGTGCAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCg > 1:278439/1‑149 (MQ=255) cGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGctct < 3:28332/149‑1 (MQ=255) cAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTAt > 3:311577/1‑149 (MQ=255) aaaGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATc < 3:247801/149‑1 (MQ=255) aaGGCGCATGTGATCAGCTGGGTCGATGGACATTCGCTGGCTGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTATCCGCGATTATATCCGCTCTGATTATCt < 6:198033/149‑1 (MQ=255) ggCGCATGTGTTCAGCTGCGTCGATGGAGTTATGCTGGCGGTAAGCATTGGAAATGCCTGATGACTCGTCGTGAGATTTTGCGTGAATTCTGTTTGTGGGGGCAACAGTGTGCAGAGTAACGGCGATTATATCCGCTCTGATTATCTGc < 3:65863/149‑1 (MQ=255) aTCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCa < 4:311577/149‑1 (MQ=255) aGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTAt > 8:105390/1‑132 (MQ=255) aGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTAt < 7:105390/132‑1 (MQ=255) gCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGt > 7:274063/1‑148 (MQ=255) gcgggACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATgg < 6:35758/145‑1 (MQ=255) tGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCg > 8:352744/1‑149 (MQ=255) aCATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTc < 2:9859/149‑1 (MQ=255) aTACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGCCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCaa > 7:103696/1‑149 (MQ=255) tACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTTTCATGGGCGGTCAAt > 1:199497/1‑149 (MQ=255) gTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGt > 5:123472/1‑148 (MQ=255) tAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGc > 7:117125/1‑131 (MQ=255) tAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGc < 8:117125/131‑1 (MQ=255) cAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAAc > 5:133371/1‑149 (MQ=255) gAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAAt < 8:274063/149‑1 (MQ=255) aaGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCAGTCAATCATTACCGTGATAACGCAATCg > 7:220323/1‑149 (MQ=255) gCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATAATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGcc > 4:317541/1‑148 (MQ=255) gatgaCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGAATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACa > 6:187840/1‑148 (MQ=255) atgaCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACaaa < 6:123472/149‑1 (MQ=255) gtcgtGAGATTGAGCGTGAAGTCTCTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGATCAATAATTAACTTGATAAATCAATAGCCTACAAAGTCAAgc > 7:353498/1‑149 (MQ=255) ttGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCa < 5:319191/105‑1 (MQ=255) ttGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCa > 6:319191/1‑105 (MQ=255) ttGAGCGTGAAGTCTGTTTGTGCCGCCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTCATTATCTGCATCATCGTCAGCCTGTCATGGGCGGGCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCACCGCGACaaa > 5:4801/1‑148 (MQ=37) tGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGc > 3:255292/1‑146 (MQ=255) aaGTCTGTTTGTTGGGGCAACAGGGAGCAGAGTCACCGCGATTATATACGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCaga < 6:105064/148‑1 (MQ=255) gtattgttGTGGGGGCAACAGTGAGCAGAGTAACCGAGAGTATATCCGCTCTGATTATCTGCATCATCTTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGaa < 2:199497/143‑1 (MQ=25) ggggCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc < 4:255292/149‑1 (MQ=255) gCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGaa < 5:148594/149‑1 (MQ=255) gatgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGCTTATCTGCATCATCGTCTGCCTGTCATGGTCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAAc > 2:272482/7‑149 (MQ=11) atgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGc < 3:203957/121‑1 (MQ=11) atgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa > 4:315251/6‑128 (MQ=11) atgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa < 3:315251/123‑1 (MQ=11) tgtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCCTCGTCTGCCTGTCATGGGCGGTTTATCATTCCCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAAcgc > 3:300436/5‑149 (MQ=11) gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCaga > 1:183910/4‑128 (MQ=11) gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCaga < 2:183910/125‑1 (MQ=11) gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAAc < 1:217212/129‑1 (MQ=11) gtatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAAc > 2:217212/4‑132 (MQ=11) cAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAAcgcg > 7:344002/1‑149 (MQ=255) tatGAGCAGAGTCACCGCGATTATCTCCGCTCTGTTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATACCGGCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAAcgcg < 6:191777/146‑1 (MQ=11) tatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa < 5:105395/123‑1 (MQ=11) tatGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa > 6:105395/3‑125 (MQ=11) gTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCgg > 6:174608/1‑148 (MQ=255) tGAGCAGAGTCACCGCGATTATCTCCGCTCTGTTTATCTGCATCATCGTCTGCCTGTCATTGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGc < 6:109236/148‑1 (MQ=11) tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag > 3:311728/1‑126 (MQ=11) tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag < 4:311728/126‑1 (MQ=11) tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCtg > 3:230302/1‑136 (MQ=11) tGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCtg < 4:230302/137‑2 (MQ=11) gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGa < 7:210068/131‑1 (MQ=11) gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGa > 8:210068/1‑131 (MQ=11) gAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCtg < 2:218203/136‑2 (MQ=11) aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGc > 3:222005/1‑86 (MQ=11) aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGc < 4:222005/86‑1 (MQ=11) aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag < 6:6182/124‑1 (MQ=11) aGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCagag > 5:6182/1‑124 (MQ=11) aGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCAGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAAt < 8:220323/149‑1 (MQ=255) gCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAgcgc > 8:65908/1‑107 (MQ=11) agagTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCTGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc < 6:91822/134‑1 (MQ=255) agagTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc > 5:91822/1‑134 (MQ=255) ccGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATg > 2:96416/1‑149 (MQ=11) cgcgATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATg > 8:169899/1‑148 (MQ=255) | AGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATG > NC_000913/577986‑578268 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |