Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,560,991 | T→A | 18.1% | intergenic (+93/‑377) | yffL → / → yffM | CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,560,991 | 0 | T | A | 18.1% | 100.2 / 5.1 | 52 | intergenic (+93/‑377) | yffL/yffM | CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein |
Reads supporting (aligned to +/- strand): ref base T (24/18); new base A (3/6); total (27/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.76e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAG > NC_000913/2560845‑2561129 | gAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGG‑TCAAAATTTg > 2:195879/1‑149 (MQ=255) gAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGGGGG‑TCAAAAGTTg < 6:206373/149‑1 (MQ=255) aGGGAGTGCTACGAATACCTGTTATTATTCAATCCGTACAACATGAATTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCGTTGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAATGTTGcc < 5:85351/148‑1 (MQ=255) ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGTTGGG‑TCAAAAGTTGCAGt > 5:87320/1‑148 (MQ=255) tGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑Agg > 6:362613/1‑149 (MQ=255) cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGt < 5:362613/148‑1 (MQ=255) cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑GCAAAAGTTGCCGTTA‑AGGt < 1:399641/148‑1 (MQ=255) aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTgg > 1:356036/1‑149 (MQ=255) cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGt < 8:18518/149‑1 (MQ=255) cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGt < 7:147270/149‑1 (MQ=255) aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGTGCTATGTTCGCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATTGTGATGGG‑GCAAAAGATGGCGTTA‑AGGGGTGTGaaat > 4:282561/1‑148 (MQ=255) aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTACATGGATTCAGGTGTGTCAAAAGTTGACGTTATTTCTGTG‑TAAAAAGTTGAAGTTA‑AGgtgtgtaatat > 2:86000/1‑144 (MQ=255) cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCaaaatt > 5:288752/1‑147 (MQ=255) ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGTGCAATGTTCTAGCCTTTTCTTTTCCTTAATTACATGAATTCAGGTGGGTCAAAAGTTGACGTTAGTGGTGTT‑ACAAAAGTTGCCGTCA‑AGGTGGGTCTACAGTTg > 7:350182/1‑148 (MQ=255) ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGc < 7:11988/149‑1 (MQ=255) ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTAAGGTGGGTCAAAAGTTGCCGTTTATGGTTTG‑TAATATTTTGACTTTA‑AGGTGGGTGaatttttt > 3:308487/1‑142 (MQ=255) ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTAAGGTGGGTCAAAAGTTGCCGTTAGTGATGTG‑TAAAAATTTGCAGTTA‑AGGAGTGTCAAGAGTTg > 5:336177/1‑148 (MQ=255) attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCaaaa < 1:11061/142‑1 (MQ=255) attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAACAGTTGcgg > 5:211267/1‑147 (MQ=255) attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTACCCTTAATTCCATGAATTCAGGTGGCTCAAAAGTTTCCGTAAGTTGTGGG‑TCAAAAGATGCCGTTA‑AGGTGGGTCaaaa > 2:11061/1‑142 (MQ=255) ttattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAATGTTGCCTTAA‑AGGTGGGTCAAAATTTg > 1:205995/1‑145 (MQ=255) ttattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCATAAATTACATGAATTCAGGTCGGTCAAAAGTTGCCGTTAGTGGTGTG‑TCAAATTATGCCGTAA‑TGGGGGGTCAATAGTTgactg > 5:334095/1‑145 (MQ=255) aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TAAAAATTTGCAGTTA‑AGGTGGGTCAAAAGTTGCAGTTTAggg > 6:107979/1‑148 (MQ=255) aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGAGTGGTGGG‑TCAAAAGTTGCCGTTT‑AGGTGGGTCAAAAGTTGCCgagtaggg > 5:327987/1‑142 (MQ=255) gcagTACACCATGAAATAAAGGAATCAGTAGTGCAATTTTCACGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTgg < 4:328964/146‑1 (MQ=37) ggAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAATTTGCCGTTA‑AGGTGGGTCAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGCCGAGTAAg > 8:45963/1‑149 (MQ=255) ggcggcCACTTTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGTGGCGTCAAAAGTTGCCGTTAGTGGTGGG‑GCAAAAGTTGCCTTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta < 5:108021/149‑1 (MQ=37) ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TAAAAAGTTGCGGTGA‑AGGTGGCGCAAAATTTGCCTTTAATGTGGGGGAAAATTTTCCTAGTTAGAATCctc > 3:23727/1‑148 (MQ=255) ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTACATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGGGGG‑TAAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta < 7:377550/149‑1 (MQ=255) tttGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGCTGAAAATTTGACTAGTAAGAATTCTACAATTAAGAATTCAAGGCGACgg > 1:120079/1‑148 (MQ=37) tGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAATTTTGCAGTTA‑AGGTGTGTCAAAAGTTGCCGTTTAGGTGGGTGTAAAGTTGCCGAGTAAGGATTCTACTATTAATAATTCTAGGCGACGGGat > 1:344178/1‑148 (MQ=37) gCCTTAATTCCATGAATTCAGGTGAGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGt > 7:227476/1‑149 (MQ=255) aTTCCATGAATTCAGCTGGGTAAAAATGTGCCTTTATTGGTGGG‑TCAAAAGTTGCAGTTA‑ATGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa < 5:176133/149‑1 (MQ=37) cATGAATTCAGGGGGCTCAAAATAAGGCCTTAGTGGGGGG‑TCAAAAGATGCCGTTA‑CGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCt < 1:195879/149‑1 (MQ=25) aTGAATTCAGTTTGGTAAAAATAAGCCTTTATGGGTGGG‑TCATATCTTGCCGTTT‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTc < 8:353517/149‑1 (MQ=25) ggTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCGGAGTAAGAATTGTACTATTAAGAATTGTAGGCGACGGGAGTGTTTTAGCGTCCcggtcgg > 1:124521/1‑148 (MQ=37) ctaaaaGATGCCGTTATTGTTGCC‑TAAAAAGCAGCCGTTA‑ATGTGGGGCAAAAGTAGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATaga < 1:86000/147‑1 (MQ=255) ccAAAAGGAGCCTTTAGTGGTGGG‑TAAAAAGTTGCCTTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATaga < 5:385477/148‑1 (MQ=38) cacgTTGCCGTTAGTGGGGGG‑TCAAAACTTGCCGTCA‑CGGGGGGTCAAAAGTTGCCTTTAAGGTGCGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagaga < 7:340837/145‑1 (MQ=255) aaaGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagagag > 1:316488/1‑149 (MQ=255) aGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGGCGTGAAGGGGGGTCAAAAGGTGCCGAGATAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTc > 2:371509/1‑149 (MQ=38) aggccttACTGGGGGG‑TAAAAACAGGCCTTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGcc < 3:368017/144‑1 (MQ=38) gcGTTACGGGTGGG‑TCACAACACGCCGTTA‑CGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa < 7:364276/148‑1 (MQ=38) cccttATTGTTGGG‑TAAATACATGCCTTTA‑AGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa < 2:356036/146‑1 (MQ=38) tttctggtgcG‑GCAAAAGTGTCCGTTA‑AGTTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTagag < 7:323/145‑1 (MQ=38) ttAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAAGTGGGTCAAAAGTTGACGAGGAAGAATTCTAATATTAAGAATTCTAGCCGACTGGACTTTATAAGCATCAAGGTCCGTCATATAGAGTCGCATagat > 6:178068/1‑148 (MQ=21) taaGTGGTGCG‑TCAAAAGTTCACGATA‑ACGTGGGTCAAAAGATGCCGTAAACGAGGCTCAAAAGTTCCCCACTAAGAATACTACAATAAAGAATTCTCCACGACAGCACTTTATTAGCCTCCAGTACAGTCATAGAGAGTCGCCTagag > 6:183829/3‑149 (MQ=255) cttttggtggGCTCATATGATGCCTTTT‑ATGTGGGTCAAAAGTTGCACTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga < 2:46073/145‑1 (MQ=38) attcaggtggG‑TCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga < 7:40183/144‑1 (MQ=255) tggG‑TCAAAAGTTTCAGTTA‑AGGTGGGTCAATAGTTGCCGTTAAGGTTGGTCAAATGTTGCCGAGTAATTTTTCTACTATTAAGAATTCTAGGCGACGGGAGGTTTTTAGCCTCCCGGTCGGTCAtt < 1:389939/127‑2 (MQ=255) gggg‑TCATAAGATGCCTTTA‑AGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACa < 3:282561/148‑1 (MQ=255) tCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCACGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag > 2:103357/1‑146 (MQ=255) | GAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGG‑TCAAAAGTTGCCGTTA‑AGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAG > NC_000913/2560845‑2561129 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |